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Résumé. — À partir des considérations développées dans le travail pré 

les conditions initiales exactes auxquelles doit être soumise la densité es EE nn de 
sence w(x, y, z) des positions d’un microcorpuscule u. On montre que la fonction f(x, y, z) vérifiant 
la relation |}? — w, doit satisfaire à une équation différentielle du type de Sturm-Liouville. Les 
conditions de départ déterminent d’une manière unique ses coefficients et l’on trouve ainsi pour Î 
une équation de probabilité identique à l’équation correspondante de Schrôdinger. On généralise 
ces résultats pour les cas suivants : problème à deux et à trois dimensions, cas d'existence des 
potentiels électromagnétiques et de plusieurs corpuscules, cas où æ dépend du temps. Dans tous 
ces cas on aboutit à l’équation correspondante de Schrôdinger. 


… Abstract. — Taking into consideration his previous work [1] the author formulates the exact 
initial conditions to which is subjected the probability density of occurence of positions w{x, y, z) 
of a microcorpuscule u. It is shown that the function f(x, y, z) in which |f|? = w, should satisfy 
a differential equation of the Sturm-Liouville type. The initial conditions determine its coeffi- 
cients in a unique Way, thus getting a probability equation for f, identical to Schrüdinger’s equa- 
tion for this case. The following cases are summarized : à problem with two or three dimensions 
in case of presence of electromagnetic potentials, in case of dependence of time, and in case of 
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numerous corpuscules. 


In all those cases one comes to the corresponding Schrôdinger equation. 


Dans une publication précédente [1] nous avons 
exposé une nouvelle hypothèse physique sur la 


structure du champ et ses conséquences quali- 


tatives pour le mouvement d’un microscorpuscule , 
le mouvement ayant un caractère statistique. Sa 
description a été ramenée à une fonction w(x, y, 
z) [1] déterminant la probabilité de présence w dr 
de & dans un volume dr = dx dy dz. Pour avoir 
la description quantitative de la conduite de u on 
doit trouver la fonction w(x, y, z). Voilà pourquoi 
nous allons tout d’abord préciser nos affirmations 
de départ, sur lesquelles sera fondée la partie ma- 
thématique de la théorie. 


A. Conditions de départ déterminant ww. 

Considérons d’abord le mouvement d’un corpus- 
cule y sur l’axe OX. Admettons quele mouvement 
classique de & s’effectue sous l’action d’une force 


: à du 
attractive conservative 5 = — Sn provenant de O 
48 


et que le mouvement est limité, pour une énergie 
totale donnée, entre (— x9, æ), une oscillation 
durant un temps de l’ordre de T,. Admettons qu’il 
n'existe dans le champ qu’une seule formation ®; 
des corpuscules AS du subvac. L’interaction de tr 
avec ®, en tenant compte des fluctuations de ®, 
et des AS aura pour résultat, comme nous l’avons 
vu dans [1], que les positions de à sur OX auront 
une distribution aléatoire de densité de proba- 
bilité w,(x), donnée par [1, (2)] 


wa) = Jim à. (1) 
es ACER (1°) 
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où At est le temps de séjour de y dans l’intervalle x, 
æ + dx, et T, le temps d'observation. (Mais il faut 
admettre, comme dans tous les cas semblables 
d'application des probabilités, qu’en pratique w (1) 
conserve un sens aussi pour un T'fini, avec T > T). 
Nous dirons alors que u se trouve dans un état 
stationnaire À, de fonction de distribution w,(x) (1) 
avec condition de normalisation (1') et d’énergie 
totale fixe Æ, (énergie moyenne de pa). Dans cet 
état stationnaire le corpuscule pourrait rester un 
temps illimité. Donc la notion d’état stationnaire a 
un sens statistique. | 

Admettons que sous l’action, par exemple, d’un 
photon ®,, qui a créé avec ®, une nouvelle for- 
mation stable ®,, le corpuscule p arrive dans un 
nouvel état stationnaire A, d’énergie £, et de 
fonction de distribution w,(x) définie aussi par (1), 
(1°). & pourrait aussi rester dans l’état À, un temps 
illimité. 

Admettons maintenant que pour des raisons 
intérieures et par échange avec le milieu extérieur, 
le système (x, ®,) émette un photon ®,,, puis, 
après un certain temps, le reçoive de nouveau, etc. 
Il se créera de cette façon un équilibre statistique 
avec le milieu environnant grâce aux photons ®,, 
échangés. Ce phénomène pourrait être décrit par 
le schéma suivant. 

Soient {, un moment initial arbitraire et #;, 
Lo, ..., & — des moments successifs fixes, avec 


q 
All > tp le tit 1 00ient 
1 


T,, T, des intervalles de temps de l’ordre de gran- 
deur des périodes moyennes du corpuscule, respec- 
tivement dans les états A,, A,, en admettant que, 
chaque intervalle At, soit grand par rapport à T, 
T,. Donc, durant l'intervalle At, le corpuscule 
pourrait passer plusieurs fois par les états 4,, 4, 
de temps moyens de séjour correspondants 1, 
resp. {,, qui dépendent sûrement du système conte- 
nant le corpuscule et d’autre part, du milieu exté- 
rieur participant à l’échange mentionné. Nous 
admettons aussi que le passage de À, à 4,, et 
l’inverse, s’effectue vite, pendant un intervalle de 
temps petit par rapport à Af,. À cause de la dis- 
tribution statistique 1l y aura une probabilité cons- 
tante a? de trouver le corpuscule en l’état À,, et 
aussi une probabilité a? pour l’état A, avec 
di + ai — | évidemment. Mais pendant l’inter- 
valle Af, à un endroit fixe sur OX, ou mieux dans 
l'intervalle +, x + dx, le corpuscule séjourne seule- 
ment pendant un temps At, € Aty. D’après (1) 
la probabilité de présence de uw en dx sera 
Az At: [T, qui dépendra en général du moment 4,, 
donc de l'intervalle At, aussi. On pourrait done 
représenter la densité de la probabilité de présence 
par la fonction w(x, t) de x et de t, où t est une 
certaine valeur du temps dans l’intervalle A,. 
Notons ici une remarque importante, valable pour 
toutes les considérations suivantes contenant le 


temps. Quoique dans les formules qui suivent, le 
temps figure comme un paramètre continu, par 
rapport auquel on prend des dérivés etc., la défi- 
nition même de w(x, t) d’après (1) a un sens déter- 
miné non en un moment fixe {, mais en un inter- 
valle de temps contenant { Voilà pourquoi les 
différences At, introduites plus haut ne peuvent pas 
être rendues infiniment petites (A4, > Ti, T3). 
Elles ont le caractère de « différentielles physiques » 
semblables aux intervalles Af du mouvement 
Brownien ou des différentielles dn de la physique 
statistique (n, nombre des particules en un volume 
donné), où dn << n, mais di représente toujours un 
grand nombre de particules. 

La fonction w(x, t) pour chaque moment { (c’est- 
à-dire pour chaque intervalle At,) doit satisfaire à 
la condition 


fut, L) Ca ui + aè — À (2) 


exprimant le fait que 1 se trouve certainement 
dans l’un des deux états À, ou À,. (L'intégrale 
de (2), comme toutes les intégrales suivantes, est 
prise dans tout le domaine de variation de la 
variable d’espace.) Nous allons en outre exiger de 
la fonction w(x, t) la propriété suivante : sa valeur 
moyenne par rapport à { pendant un intervalle T 
contenant plusieurs intervalles At, ne doit dé- 
pendre que de x, et pour T — co elle doit repré- 
senter une somme des fonctions de distribution 
stationnaire af w,(x) et a? w,(x), correspondant aux 
états À,, A,. Cette exigence est fondée sur le mo- 
dèle physique dont nous nous servons, à savoir : 
quel que soit le mode de passage du corpuscule & 
au cours du temps entre les états À,, À, le carac- 
tère statistique du phénomène indique que l’effet 
des passages aléatoires doit s’effacer au cours du 
temps et qu'il ne doit rester à la limite que l’effet 
des distributions stationnaires. On doit done avoir, 
indépendamment du moment initial t,, 


DORSALE ; e 
ais Sa (er, 1) dt = w(x) = aï wilx) + a w,(x). (3) 


feat te = 1. fat TEHIS ai Pt (4) 


La généralisation s'impose d’une façon naturelle 
dans le cas où u peut apparaître dans un nombre 
arbitraire d’états À,, À... A, …… , de fonctions de 
distributions correspondantes w,(x)..…., wx(x), …, si 
chacun d’eux existait seul. Si à cause de l’échange de 
photons ®,,, .…., du, avec le milieu environnant, 
des formations ®,, …, ®,... se forment et dispa- 
raissent, alcrs le corpuscule 4 se trouvera plusieurs 
fois au cours des intervalles différents de temps A4, 
dans les états différents A, Az, …, de proba- 
bilités relatives correspondantes a?, ……., a... 
Dans le cas d'équilibre statistique, la densité de 
probabilité de présence de x sera une fonc- 


Nou 


tion w(x, t) de x et de { (dans le sens indiqué plus 
haut). La fonction w(x, t) possèdera les propriétés 
suivantes représentant des généralisations 1mmé- 
diates des propriétés (2), (3), (4). 


foie, ee UE C0 Mat | () 

; { to + T o 

Jim TJ w(æ, t) dt = w(x) — > ax aglæ) (6) 
fatal FN et EVE (2) 


Il est clair que la généralisation des considé- 
rations précédentes pour les problèmes à deux et à 
trois dimensions est immédiate. Le corpuseule x 
peut se trouver dans des états différents Az d’éner- 
gie Æx et de densité de probabilité mx, y, 2). 
Dans les formules correspondant à (2)-(7) l’inté- 
gration a lieu dans tout le domaine de variations 
des variables. 

2. Considérons de nouveau l’image donnée 
en [1, B]. Émettons d’un point M, du champ pen- 
dant des intervalles de temps successifs A,, 
At, des groupes de corpuscules x de même 
vitesse initiale &,, en négligeant l'interaction entre 
eux. Les corpuscules & (ou mieux leurs centres de 
masse) décriront des trajectoires différentes qui 
rempliront l’espace d’après une certaine loi de dis- 
tribution statistique stationnaire. Leur densité 
e(x, y, z) qui ne dépend pas du temps sera propor- 
tionnelle (nous écrirons égale) à la densité de la 
probabilité de présence du corpuscule individuel 


p(æ, y, 2) = w{x, y, 2). (8) 


Si le champ U varie avec le temps, mais suffi- 
samment lentement, de façon que le caractère sta- 
tionnaire du phénomène reste invariable, alors p 
et w seront des fonctions du temps et la relation (8) 
continuera à exister pour chaque {. On pourrait 
ainsi considérer les corpuscules comme des cons- 
tituants d’un fluide de densité p(x, y, z, t) variable 
en générale pour lequel dans chaque -point ou 
mieux dans chaque petit domaine contenant le 
point (x, y, z) du champ, où la vitesse moyenne du 
fluide est v, on a l’équation classique hydrodyna- 
mique de continuité 


ee + div wvù = 0. (9) 


Les conditions (9) et (5) sont compatibles 
puisque, en intégrant (9) par rapport à x, y, z, on 
arrive à (9). | 

Dans le cas stationnaire, l’équation (9) devient 
évidemment : 


div wv — 0. (10) 


On peut en outre écrire pour un corpuscule 1 
déterminé la formule suivante, valable aux envi- 
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rons d’un point d'espace déterminé M, et dans un 
intervalle de temps At 


Ur = — (1m) grad Sx 


où 5% est la fonction classique de Jacobi pour , 
avec les restrictions mentionnées, puisque dans ces 
conditions le mouvement de y peut être considéré 
comme classique. On ne s’occupera pas ici de savoir 
comment sont liées les différentes fonctions S, des 
différents corpuscules pour des positions et des 
moments différents. 

3. Dans le cas où l’on pourrait négliger la struc- 
ture discrète du champ c’est-à-dire les forma- 
tions ®; et les fluctuations des corpuscules AS 
(alors les états voisins A3, A4, et leurs énergies 
correspondantes Æ,, Ez,, ne se distinguent que 
très peu) le mouvement du corpuscule & doit 
devenir le mouvement classique. Nous appellerons 
principe de correspondance cet état de choses, 
puisqu'il lie d’une façon naturelle les mouvements 
quantique et classique. Son sens peut être exprimé 
plus précisément de la façon suivante : 

Considérons le mouvement classique périodique 
sur OX du corpuscule à de période T,. Quoique 
son mouvement soit précisément connu à chaque 
instant {, on peut tout de même parler de proba- 
bilité en présence en se demandant quel est son 
séjour relatif dp dans l'intervalle x, x + dx où sa 
vitesse est v(x). Si son séjour réel pendant une 
période T, y est dt, — dx/v, alors on aura 


GS = cer) (11) 


La grandeur dp peut ainsi être écrite de la 
façon suivante : 


dora dr Cie lv |) (12) 


où la constante C'est déterminée par la condition de 
normalisation 


1 —C Ÿ azote) = Clos OAI M) 


CD, 


Il est clair que dans le cas classique on peut 
appeler la grandeur dp — w(x) dx probabilité de 
présence « classique » de & dans l’intervalle dx. 

Cette considération peut être étendue à un mou- 
vement quelconque. Ainsi pour un mouvement non 
périodique de y sur OX si la distance x,, x, conte- 
nant dx est traversée une fois pendant le temps T,, 
alors la formule (11) de dp conserve le même sens 
de probabilité. De même si, dans le cas de mou- 
vement à deux ou à trois dimensions, on fait des 
observations pendant le temps T et si l’on trouve 
que le corpuscule est observé une fois durant le 
temps dé < T dans un domaine dr sur un élé- 
ment d{ de sa trajectoire où sa vitesse est v, on 
peut aussi parler de probabilité de présence de x 
dans dr définie par une formule, semblable à (11), 


wt) à = (deT) — (AT) (difo). (11”) 
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Ainsi dans le cas où l’on aura le droit de consi- 
dérer le mouvement quantique de comme un 
mouvement classique, on exigera la condition 
mathématique que la distribution de probabilité w 
du mouvement quantique devienne la distribution 
classique donnée par la formule correspon- 
dante (11') ou (11). 

4. Dans le cas particulier où ÜU — C, quand le 
mouvement classique est de vitesse constante pv, 
on aura dans le cas quantique aussi un mouvement 
uniforme, comme cela a été expliqué en [1]. Il sera 
caractérisé par une probabilité quantique w, qui 
est constante, 


æ — Cte, (14) 


Au fond ce résultat est contenu comme cas parti- 
culier dans le $ 3, mais nous le soulignons à part vu 
son importance pour les applications. 

5. Dans un cas particulier quelconque où la 
théorie quantique classique donne des valeurs 
quantifiées de l’énergie qui sont confirmées par 
l'expérience, la nouvelle théorie doit redonner pour 
le même cas les mêmes valeurs. On peut prendre 
comme exemple loscillateur harmonique auquel 
Planck a appliqué pour la première fois l’idée des 
quantas. L'énergie de l’oscillateur est donnée par la 
formule de Planck E, = nhy (n = 0, 1, 2...) bien 
vérifiée par l’expérience pour de grandes valeurs 
de n. Pour des valeurs très élevées de 7 la variation 
relative de Æ, devient très petite et l’énergie tend 
vers une distribution continue. Dans le cas de 
notre modèle de mouvement d’un corpuscule appli- 
qué à l’oscillateur harmonique, on a quelque chose 
de semblable. Pour des valeurs très basses de 
l’énergie, on aura des valeurs inconnues et soit £, 
leur ordre de grandeur. Pour des valeurs très 
élevées de l’énergie de l’ordre de £” (théoriquemeat 
si £" — co) on doit retrouver les valeurs continues 
de l’énergie de l’oscillateur harmonique classique. 
Nous admettrons que si l’énergie de l’oscillateur 
harmonique quantique est grande par rapport à Æ,, 
la formule de Planck est valable, £’est-à-dire 


En = niv pour EL Er SE DU) 

6. Si un corpuscule x (un électron par exemple) 
passe d’un état stationnaire d’énergie Æ, en un 
autre d'énergie E,, la différence Eym — En — En 
est liée à la fréquence v,, d’un photon ®,, émis ou 
absorbé (c’est-à-dire de la formation corres- 
pondante ® créée ou disparue) d’après la formule 
de Planck-Bohr 


19) = RVrm. (16) 

7. Admettons que pour un problème classique le 
potentiel U(x, y, z) est une somme de potentiels 
correspondants, fonctions d’une seule variable, 
c’est-à-dire 


U(x, y, 3) = Uilx) + Uily) + Ui(z) 


— 
> 
SI 

= 


N°4 


ou en d’autres termes, que les mouvements suivant 
les trois axes sont indépendants l’un de l’autre. 
Nous admettrons que pour le problème quantique 
correspondant la fonction w(x, y, z) qui détermine 
la probabilité de présence de y de coordonnées x, 
y,z,se décompose en un produit de fonctions w;(t), 
Wo(y), Wa(z), à savoir 


we, y, 2) = wi(x) w2(y) wa(a). (18) 

Autrement dit, les trois événements, à savoir que 
les trois coordonnées x, y, z de x ont des valeurs 
données, représentent trois événements indé- 
pendants. 


B. Détermination de la fonction *. 

En se basant sur l’exposé précédent, on peut trou- 
ver l’expression mathématique exacte de la fonction 
de distribution w pour chaque problème quantique 
correspondant à un problème déterminé de la 
mécanique classique. En même temps nous admet- 
tons que les vitesses v des corpuscules sont 
petites par rapport à la vitesse de la lumière c, donc 
les considérations faites ici sont fondées sur la 
mécanique classique, et non sur la mécanique rela- 
tiviste. 


1. PROBLÈME A UNE DIMENSION. — Soit un 
microcorpuscule 1 de masse »m en mouvement sur 
l'axe OX (— co < x < co) sous l’influence d’une 
force F provenant d’un potentiel U(x) qui déter- 
mine sur OX un mouvement classique limité. Il 
lui correspond, comme on l’a vu, un mouvement 
quantique. Alors le corpuscule peut se trouver 
successivement dans différents états Az(k = 1,2...) 
décrits par des fonctions w;(x) définies par (1), et 
au cours du temps dans un état plus général décrit 
par une fonction w(x, t). Admettons que pendant 
un certain temps le corpuscule u se trouve dans un 
état fixe Az d'énergie Æ, et de densité de proba- 
bilité de présence w,(x). Il est clair que par sa 
nature même la fonction w,(x) est non négative, 
wx(x) > 0, uniforme et continue, de valeurs à la 
limite nulles, w(+ co) — 0. Nous avons aussi 
admis en[1]qu’elle possède des dérivées. Mais au lieu 
d'étudier les propriétés des fonctions non néga- 
tives, il est plus facile et plus commode de ramener 
cette étude à des fonctions prenant des valeurs 
positives et négatives. Cela peut être fait de la 
façon suivante. À chaque fonction w;(x) on peut 
faire correspondre une fonction continue et uni- 
forme /;(x) qui s’annule en même temps que w,(x) 
et qui change de signe successivement pour les 
racines de wz(x). Cette correspondance se réalise 
d’une façon très naturelle par la relation suivante 


fklx) = Vwr(x) (19) 


ou 


(19°) 
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avec la convention que f(x) est successivement 
égale à Va(x) et —V/w(x) entre les racines 
de wx(x). (La question de l’unicité de cette rela- 
tion entre w, et fL sera discutée plus bas.) 

Il est clair que /,(x), ainsi que w,(x), est une 
fonction continue et uniforme, possédant des déri- 
vées, et aussi que f{(F co) — 0. Mais au lieu de la 
fonction réelle , (19) on peut considérer une fonc- 
ton complexe plus générale f(x) donnée par la 
représentation suivante (qui se montre aussi impo- 
sée par les considérations données plus loin, con- 
duisant à la détermination définitive de la fonc- 


tion fx(x)) 


rx) — gr(x) FU AT) 


ou 
welx) — |fe(x)|? — gK(x) 


la convention du signe alternatif de \/wz restant 
valable. On admet pour la fonction f, (20) comme 
généralisation de fx (19), qu’elle est aussi continue, 
uniforme et qu’elle possède des dérivées, avec évi- 
demment f;(+ co) = 0. Il est clair que si l’on 
connaît la fonction fx (19) ou (20), on aura immé- 
diatement la fonction cherchée w(x) (19’) ou (20). 

Donc, si le corpuscule 4 se trouve dans un état 
stationnaire fixe A7, on fait correspondre à la fonc- 
tion de distribution wz(x) la fonction f(x) (20). 
Admettons maintenant qu’au cours du temps le 
corpuscule passe successivement par des états diffé- 
rents Ax(k — 1, 2,...), comme cela a été expliqué 
en À,p.201. Alors la densité de probabilité de pré- 
sence totale sera évidemment une fonction w(x, t) 
de x et de t. On est naturellement conduit à la 
question : Comment la fonction w(x, t) peut-elle être 
représentée à l’aide des fonctions individuelles f(x) 
auxquelles correspondent les probabilités rela- 
tives a% ? Les considérations précédentes amènent 
à exiger que la fonction w(x, t) remplissse les condi- 
tions suivantes : (a) Si & se trouve dans un état 
fixe Az, c’est-à-dire tous les & = 0 (1 k) et 
a? — 1, alors (x, t) doit devenir la fonction w;(x) 
correspondante ; (b) w(x, t) doit satisfare à la 
condition (5) ; (c) w(x, t) doit satisfaire à la condi- 
tion (6). 

Dans l'esprit de ce qui a été dit, on doit ramener 
la fonction w(x, t) à une fonction correspon- 
dante f(x, t), la relation entre elles étant, naturel- 
lement, aussi de la forme (20”), 


(x, t) = |f(, t)[? (21) 


où f(x, t) dépend d’une certaine façon des fonc- 
tions f(x) (k — 1, 2,...). Les conditions (5)-(7) 
suggèrent immédiatement la relation la plus simple 
suivante entre f(x, t) et toutes les autres f(x) (la 
question d’unicité sera discutée plus tard) 


fl, à = 2 ar Cr(t) fe(x) (22) 
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où les Cz(t) sont des fonctions complexes de t avec 
la condition |C;(t)| = 1, c’est-à-dire (9, — fonction 
réelle) 


Clt) En errt) (22°) 


Il est facile de voir que la fonction w (21) ainsi 
construite satisfait aux conditions exigées (a)-(c). 
Ainsi si un nombre déterminé 4 — 1, et si tous les 
autres 4 — 0, on aura de (21), (22’) 

w(x, 1) = we(x) = |fr(x)|? (21) 
donc la condition (a) est remplie. 
. Remplaçons l'expression / (22) dans (21) en uti- 
Hisant (22'.) On aura 


mx, 0 = f*f = D aklfr/? 
+ à ak 4 Cri) Ci) fu) f(x). (23) 


D’après le choix (22') des fonctions Cz(t) la 
condition (c) est satisfaite aussi. 

Remplaçons maintenant w (23) dans la condition 
de normalisation (5). On trouve 


1 = foie, a) dhe = > a J inte)|* dx 


e a a Chlt) Cult) frite) filæ) dæ (24) 
(24) sera identique à (5) si l’on a 


ft) Fe de oe 0 


HEC dæ = fuite) Ci = A 


Les (25) expriment les conditions que toutes les 
fonctions f,(x) sont orthogonales entre elles, tandis 
que les (25') sont les conditions de normalisation 
acceptées (7). Ainsi la relation (24) devient 


IEC t) dx = D ag — 1. 


(25) 


(25°) 


(24°) 


Par conséquent la condition (b) est aussi remplie 
pour chaque £& (c’est-à-dire pour chaque inter- 
valle At, comme cela a été dit dans À, p. 201). 

On arrive donc à la conclusion que les fonc- 
tions f(x) déterminant la distribution des proba- 
bilités de présence dans les états stationnaires 4; 
remplissent les conditions suivantes : elles sont 
continues et uniformes, elles remplissent des condi- 
tions à la limite aux deux extrémités du domaine, 
«K( co) = 0, elles possèdent des dérivées succes- 
sives. Elles sont de plus orthogonales entre elles et 
peuvent être normées. D’un autre côté, chacune 
des fonctions f; correspond à une valeur fixée Æ, 
de l’énergie Æ. Ces propriétés suggèrent évidem- 
ment la pensée de considérer ces f; et Ex comme 
des solutions propres et des valeurs propres d’une 
équation différentielle d’ordre IT (à cause des 
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deux conditions à la limite aux extrémités de l’in- 
tervalle x — 00) pour la fonction f(x), de coef- 
ficients variables. Une équation de cette nature ne 
peut être que du type de Sturm-Liouville dont la 
forme la plus générale est la suivante 


ER dr enr 


26 
EC 072 (0 


f= 0. 


dx? 


Ici les coefficients b,, b, sont des fonctions incon- 
nues de x, dépendant sûrement aussi du poten- 
tiel U(x) qui détermine chaque problème concret, 
c'est-à-dire ba —0.(&, U(x)), bb: (x, UV)" 
est un paramètre réel introduit pour raison des 
dimensions alt= jm etre) 

Les valeurs propres attendues de équation (26) 
doivent évidemment être réelles. Mais la théorie 
des équations différentielles [2], [3] montre que 
cela nécessite que lopérateur linéaire appliqué sur / 
dans la partie gauche de (26) soit égale à son 
propre adjoint. Cela est possible dans les deux cas 
suivants : /A,) Les fonctions b, et b, sont réelles. 
(A,) b, est une fonction purement imaginaire et 
b,—réelle. Nous ne considérerons maintenant que le 
cas (A,) permettant de résoudre jusqu’au bout le 
problème posé, en laissant pour plus tard l’hypo- 
thèse (A,) qui mène à une autre variante, et plus 
générale, du problème. 

L’équation (26) écrite pour caractériser des états 
discrets doit aussi décrire le mouvement quantique 
du corpuscule & pour tous les cas, c’est-à-dire pour 
chaque potentiel U(x). Mais pour U(x) — U0 — Cie 
on a dans le cas classique un mouvement uniforme. 
Il doit lui correspondre un mouvement uniforme 
dans le cas quantique comme on l’a vu dans [1]. 
On doit avoir alors æ — Cf, d’où il suit que l’équa- 
tion (26) doit admettre une solution de la forme 

RENE) (27) 
(o(x) réelle). Cela n’est possible que si b, = 0, 
b, — 0} — CC, ce qui donne 
8 — À VAE CL 


a 


pla B, (E + 65 > 0). (28) 

Mais puisque b, dépend de U°, considéré main- 
tenant comme paramètre, et puisqu'on doit avoir 
b, = 0 pour chaque valeur de U9, on en conclut 
que b, = 0. (Le même résultat b, = 0 suit aussi 
de la théorie générale des équations différen- 
telles [3].) Ainsi l’équation (26) devient-elle 


(29) 


contenant dans son coefficient la fonction inconnue 
unique b, et le paramètre inconnu «. 

D’après la théorie générale [2], [3] l'équation (29) 
possède la propriété suivante. Si on l’applique 
pour une fonction b,(x) continue et pour un inter- 
valle borné x,, x,, avec des conditions à la limite 
linéaires aux extrémités x,, x, elle possède un 
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spectre discontinu Æ, par rapport au paramètre Æ, 
où |Æ,| — co si n — co. On sait également pour 
les valeurs asymptotiques des valeurs propres 
que (En — Ey_1)/En —0 avec n — co. Dans 
notre cas l'intervalle x,, x, est infini, æ, = — co, 
%, = 00 avec f(F co) = 0.Si — b,(x) > M > 0 pour 
%, —>— 00, æ, —> co (M donnée), l’équation (29) 
possède un spectre discontinu [2]. Nous avons 
montré dans un travail non encore publié (en consi- 
dérant cette propriété comme inconnue) que si 
b,(x) — — co pour x — Æ ©, donc le spectre 
de (28) est discontinu, on a encore, comme pour 
l'intervalle fini, la propriété (E, — Ey_;)/Er — 0. 
pour |Æ,] — co. Par conséquent, pour des E, 
grandes, c’est-à-dire pour des valeurs macros- 
copiques de l’énergie, les niveaux énergétiques de- 
viennent relativement très denses et pratiquement 
la situation n’est pas différente du cas classique où 
l’énergie prend des valeurs continues. La même 
propriété existe aussi pour des valeurs fimies de };, 
mais pour des valeurs suffisamment petites du pa- 
ramètre « dans (29). Par conséquent l’équation (29) 
pour des valeurs de « suffisamment petites doit 
donner des résultats qui tendent vers les résultats 
classiques, c’est-à-dire sa solution f(x) doit déter- 
miner une fonction w(x) — [f(x)|? qui tend vers la 
probabilité classique æ (11). On utilisera cette pro- 
priété pour déterminer la fonction inconnue b, 
dans (29). 
Posons en (29) pour simplifier 


BRL (30) 
et cherchons une solution de (29) de la forme 
pee (31) 
On trouve de (31) successivement f', f", et en 
les remplaçant dans (29), on aura 
y? + b + ay" = 0. (32) 


En considérant « comme petit paramètre on peut 
développer y(x) en série suivant les puissances de «, 


Y(x) = Yolx) +ioyalx) + ... + ot ynfx) + ..… (33) 


On tire de cette équation y, y”, et en les rem- 
plaçant dans (32), on aura 
(YO RU TE RE DEEE y ee la 0 4) 
(34) détermine les approximations successives. En 
se bornant à y,, y, (n — 1) on aura 


on RU ne À (35) 
d’où 
nn = SUV OS ie [Vo ax 
’ 1 d ’ CITE 
D = 2 de ya(e) = — In \/b(x). (36) 
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En remplaçant y, y, dans (33), (31), on aura 
approximation cherchée de la solution 


: 1 +aJ Vitide, 


, */B(x) FA 


D'un autre côté, dans le cas classique, on a 
d’après (11) 
PÉCIn=E Nm CV UE). C = 1/T,. | (38) 


En égalant w (38) à w — |f|? (37) et en prenant 
en considération (30), on peut écrire 


VmCIV2(E — U(r)) = 1[Vb{a) = 1JVh, TE 
d’où 
bic) + E = (2/mC?) (E — U(x)). (39) 


Cette expression remplacée dans l'équation (29) 
donne 


f + A(E— Ul(x))f—=0, A—92/mCto, (40) 
On à ainsi l’équation déterminant la fonction de 
distribution w — |f[? avec le seul paramètre in- 


connu À dans son coefficient. 

Afin de déterminer À appliquons l'équation (40) 
dans un cas particulier, à savoir à l’oscillateur har- 
monique (À, p.). On a dans ce cas 


Ulx) —=+2?, k = 4n2v2m (41) 
(v, — fréquence de l’oscillateur). En remplaçant 
U (41) dans (40) et en posant 
VAR) x =E€, VRAIRE = à (42) 
on trouve l’équation 
(MEL ANA ENT 0. (43) 


Comme il est connu cette équation possède des 
solutions continues, uniformes et s’annulant à 
l'infini (des polynômes d’Hermite) pour des valeurs 
de À formant un spectre discontinu, À = 2n +1 
(n — 0, 1, 2,..). En remplaçant ces dernières 
en (42), on trouve pour Æ À 


En = Vk]2A (2n + 1) = 2rv V2mfA (re Le J) (44) 


Comme nous avons dit en (A, p.), dans ce cas 
particulier nous exigeons que les valeurs élevées 
de E, (E KE < E, (15)) coïncident avec celles 
qui sont données par la formule de Planck E, = nhy, 
bien vérifiée par l’expérience. En comparant la 
dernière valeur de Æ, avec celle donnée par (44), 
dans laquelle on néglige 1/2 devant n (n grand), 
on trouve 


2rv VomfA — nhv (45) 


d’où, pour le paramètre inconnu À 


A —8Tr2 mr "2m? ñ — h[2r. (46) 
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En le remplaçant en (40) on aura finalement 
d? , 
“ (E — Ufx)) f = 0. (47) 


C’est l'équation cherchée déterminant la fonc- 
tion f(x) et par conséquent aussi la fonction de 
distribution w — |f|? dans le cas linéaire, voilà 
pourquoi nous l’appellerons ici, pour unicité, équa- 
tion de distribution de la probabilité ou simplement 
équation de probabilité. C’est évidemment l’équa- 
tion bien connue d'amplitude de Schrüdinger pour 
ce cas. Nous acceptons la validité de (47) pour 
chaque potentiel U(x) donné, quoique en l’établis- 
sant nous ayons admis certaines restrictions sur 
l'allure asymptotique de la fonction b,, done 
pour Ü aussi. On peut aussi voir que pour d’autres 
allures de la fonction U les propriétés de f, déter- 
minées par (47), sont en accord avec l’idée phy- 
sique qui est à la base de (47). 


Remarque. — Nous avons ainsi vu qu’en intro- 
duisant au lieu de la fonction w(x,t) la fonction 
f(x, t), resp. les fonctions f(x) à l’aide des rela- 
tions (21), (22), (20), l’étude de notre problème 
physique se ramène à la solution de l’équation 
différentielle (47). Mais peut-on dire que cela est 
la façon unique de faire correspondre à la fonction 
non négative w(x,t) la fonction de signe va- 
riable f(x,t) ou la fonction complexe /(x,t) et 
qu’une autre correspondance ne ramènera pas le 
problème à un autre appareil mathématique appro- 
prié ? Nous allons montrer que la solution donnée 
est unique dans le domaine d’une classe assez large 
de fonctions, à savoir des fonctions rationnelles 
entières, où, d’après nous, il est naturel de chercher 
cette solution. 

Supposons qu’il existe entre les fonctions w et f 
la relation quadratique (20') (resp. 21), mais qu’au 
lieu de l’égalité (22) exprimant f comme fonction 
des fonctions f; et des paramètres ay, il existe une 
autre relation. Puisque (22) satisfait aux conditions 
physiques préliminaires (5)-(7), cette nouvelle rela- 
tion, éventuellement plus générale que (22), doit la 
contenir comme cas particulier. Voilà pourquoi on 
peut admettre que f est une fonction des fx et ax 
telle que dans le développement de f suivant fr, 4x, 
l'expression linéaire (22) ne soit qu’un premier 
terme d’un développement plus général, à savoir 


f(x; 1) = 2 Ac Cxlt) fr(x) Le of, d:, Go, Cr fu 15e 0 °)) (22°) 


où le terme supplémentaire s’annulesi tous lesnom- 
bres 4, sauf l’un quelconque 4,, s’annulent. Aïnsi f 
(22/)satisfait-il à la première condition (a) imposée 
sur w. D'un autre côté en formant l’expression 
w — |ff?, parmi les autres membres on trouve 
aussi le terme |o|?. La condition (5) (ou 24) montre 
de sa part que l'intégrale de |@|? par rapport à x 
doit être nulle, et la condition (6) — que la valeur 
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moyenne de |o[? par rapport à t doit être nulle 
pour chaque t,, d’où il suit que o = 0. 

Admettons maintenant que la relation entre w 
et la fonction cherchée f soit donnée, au lieu de (21), 
par 


Mo) = eee 2e) (2/1) 
et la relation entre f et les f,— par (22) ou (22°). 
Alors les conditions (5), (6) mènent, comme il est. 
facile de le vérifier, à des conditions que des inté- 
grales par rapport à x des expressions de la forme 
D ATEN nr qi le 2;.) Fontiesales atzero, 
et que la valeur moyenne par rapport à {de fonc- 
tions positives de { dans l’intervalle (4,, co) pour 
chaque t, est aussi égale à zéro. [Il s’ensuit que la 
supposition faite n’est pas vraie. On tire aussi la 
même conclusion de l'hypothèse que w est une 
certaine fonction rationnelle entière de |f|?. Ainsi 
donc dans le domaine des fonctions continues, 
rationnelles et entières le problème indiqué de 
ramener la fonction w à une fonction de la nature 
de f n’a qu’une solution conduisant à l’équation de 
probabilité (47) pour f. 


2. CAS DES PROBLÈMES A DEUX ET A TROIS DIMEN- 
SIONS. — Considérons maintenant le mouvement 
d’un corpuscule & dans le plan (x, y) le champ 
U(x, y) étant donné. À cause des formations ® 
créées dans le champ et des fluctuations des cor- 
puscules AS, le corpuscule & se trouvera dans 
des états stationnaires différents Az d’énergies 
moyennes Z7. Nous allons encore chercher la pro- 
babilité de présence w(x, y) do de u dans le do- 
maine do — dx dy en sortant toujours de la défi- 
nition æ = lim Z Af;/T. Pour chaque état sta- 
tionnaire À, de probabilité totale constante aÿ < 1 
et de densité de probabilité w;(x, y) on introduit la 
fonction en générale complexe f, pareille à (20) 


fa(æ, y) = Vale, y) eixxtn (48) 


avec une convention analogue pour le signe de fz 
qui change si on traverse les lignes où s’annule 
wx(x, y). Pour les mêmes raisons, comme dans le 
cas de (20), on considérera les fonctions f, comme 
continues, uniformes, possédant des dérivés et 
nulles à l'infini. Au cours des créations ou des 
annulations des formations ®; le corpuscule peut 
passer par des états différents, par conséquent la 
densité de probabilité w(x, y, t) sera une fonction 
du temps aussi (dans le sens indiqué en (A, p.)). 
On aura de même des relations de la forme (5)-(7), 
à savoir 


EC y, t) dx dy = D a —1 (49) 


TT 
lim SJ (x, y, t) dé — > a% wr(æ, y) (50) 


[forte y) dr dy = 1 


(51) 
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La fonction w(x, y) peut être représentée, 
comme (21) sous la forme 


(x, y, t) = |f(x, y, 4)? (52) 

Îlæ, y, 0 = À ax Cu(t) flx, y) (53) 

où |Cxl = 1. En remplaçant w (52) dans (49) on voit 

que la condition (49) reste valable si l’on a les 

conditions suivantes d’orthogonalité et de norma- 
lisation des fx 


(54) 


[fit y) fil, 9) dx dy = 0 KZ 


fut, t)|? ae = 1, 


Par conséquent les fonctions f(x, y) peuvent 
encore être cherchées comme des fonctions propres 
correspondant aux valeurs propres Ex d’une équa- 
tion linéaire aux dérivées-partielles de second ordre. 
Elle peut être construite en suivant le schéma 
utilisé pour l’équation (47) dans le cas linéaire, 
mais elle peut aussi être directement écrite d’après 
les raisons suivantes : 1) Si la fonction U ne dépend 
pas d’une quelconque des variables x ou y, cette 
équation doit se transformer en l’équation (47) 
pour l’autre variable. 2) Si 


U(x,y) = Ui(x) + U,(y), 
on doit avoir f(x, y) — f(x) f(y), c’est-à-dire 
w(x, y) = w.(x) wA(y) d’après (A, p.203). On voit 
facilement que c’est l’équation 


(55) 


2 

Af + T5 (E — Ua, y)) f = 0. (56) 
La généralisation de ces résultats pour le cas à 
trois dimensions est immédiate. On aura pour la 
fonction f qui détermine la probabilité w = |f[? 

l’équation suivante correspondant à (56) 

2 

Af + (E— 


U(x, y; z)) Î A (57) 


3. DÉTERMINATION DE f EN FONCTION DU TEMPS. 
— Nous nous sommes occupés jusqu’à maintenant 
des états stationnaires du corpuscule y représentés 
par des fonctions f,; qui ne contiennent pas le 


temps et qui satisfont à l’équation (57). Le temps 


figure jusqu’à maintenant dans l’équation de conti- 
nuité (9) et dans la densité de probabilité w(r, t) (21) 
(on considère le cas général à trois dimensions, 
r = r(x,y,2)). Pour un état stationnaire fixe déter- 
miné par une fonction f, (20) (k donné), l’équa- 
tion (9) est satisfaite puisqu’alors dw/dt — 0 et l’on 
aura (10), div ww — 0. Mais elle sera également 
satisfaite si au lieu de la fonction f;(r) (20) on 
prend la fonction C;(t) f(:) de (22) (|Cx| = 1), 
puisque dans w — |f,|? l’effet du teinps disparaît. 
Donc la fonction C;(?) f4(r) contenant le temps peut 
décrire un état stationnaire 4; de densité de proba- 
bilité wy(r) = |fr(r)|2. Afin de déterminer la nature 
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des fonctions C;(t) en (22) on peut partir de l’indi- 
cation suivante, basée sur des raisons physiques. 
Dans la formation ®, introduite plus haut il existe 
un phénomène périodique d’une certaine nature, 
et nous avons admis que sa fréquence vw, est donnée 
par la formule de Planck Æx — hyw. Il est naturel 
d'essayer de lier le facteur C;(t) en (22) au phéno- 
mène périodique en ®; en considérant Cy}(t) comme 
une fonction périodique du temps de fréquence w», 
c’est-à-dire 


ST 
Crlt) = e TN Le À (Sgft) = —2rwt). (58) 
(On à pris devant Æ, le signe — pour conserver 


la forme des fonctions bien connues de la méca- 
nique quantique.) On aura alors pour la fonc- 
tion Cy(t) fx(r) que nous indiquerons par W, en 
l’écrivant dans les variables +, y, z, 


—#%, 


Pr, V; 3, ) TA f(x, y; z) € ê (59) 
HR) eue y che Hem, (59’) 

La fonction fx (59') satisfait à l'équation (57). La 
fonction x (59) satisfait à l’équation suivante 
(en omettant l’indice k) 


À J 
KR UY = __2imdT 
ñ ñ dt 


(60) 
que l’on vérifie immédiatement, ou bien que l’on 
trouve en éliminant £ de (57) et de la dérivée d#/à4 
tirée de (59). Nous admettons que (60) est l’équa- 
tion générale à laquelle est soumise la fonction Ÿ 
dans le cas aussi où le potentiel U contient le 
temps. Nous l’appellerons ici équation temporelle 
de distribution de la probabilité. Elle est, évidem- 
ment, identique à l’équation de Schrüdinger pour 
ce cas. 

En combinant l’équation (60) avec sa conjuguée 
on trouve, comme il est connu, l'équation de conti- 
nuité 


= + div pv = 0 (61) 
où p — wet 
_ À 1 
9 7 2im VE* 
et la condition de normalisation (en choisissant 
d’une façon convenable la constante multiplicative 


en T) 
Jirrar À. 


Nous avons admis au début l’équation (61) 
comme équation de départ (9) (où + avait le sens 
d’une vitesse moyenne) et nous avons alors écrit la 
condition de normalisation (63) sous la forme (5). 

Admettons maintenant que le potentiel U ne 
contient pas le temps. Alors la solution générale de 


(HREANEERTEENTTES) (62) 


(63) 


SUR L’INTERPRÉTATION DE LA MÉCANIQUE QUANTIQUE 209 


l'équation (60) est une combinaison linéaire de solu- 
tions de la forme (59) à savoir 


L#, Y, 3, t) = > Ale TE te Y, 2) e—"Tivrt 


(ax— des constantes réelles). La densité de la pro- 
babilité de présence sera 


(64) 


CAN AU AR ARE > aRlfelæ, y, z)[? 
2e à ar a eTivent fr(æ, Y, 2) filæ, y, z). (65) 


c'est-à-dire que w contient des membres pério- 
diques par rapport à { de fréquence vwy. (65) a 
exactement la même forme que la fonction w (23) 
que nous avons écrite pour des raisons générales, 
sans y préciser la nature des fonctions C;(1). 


4. CAS DE PRÉSENCE D'UN CHAMP ÉLECTRO- 
MAGNÉTIQUE. — Nous allons maintenant étudier la 
deuxième hypothèse possible (A,) pour le coef- 
ficient b, dans l’équation (26). 

Les résultats trouvés jusqu'ici se rapportaient à 
un corpuscule y sans charge dans un champ 
U(x, y, z). Admettons maintenant que le corpus- 
cule possède une charge e et qu’il existe outre le 
potentiel U un champ électromagnétique donné 
par le potentiel scalaire o(x, y, z) et le potentiel 
vecteur A(x, y, 2). Si À = 0, le cas ne sera pas 
essentiellement différent du précédent et pour 
englober le cas présent il suffit de remplacer sim- 
plement U par U, = U + eo dans l’équation (57), 
resp. (47). Nous allons maintenant faire la géné- 
ralisation pour le cas où o et À sont différents de 
zéro. On aura alors les relations bien connues de 
l’électrodynamique 

2 
PI 7 (p—°4) + eo + U, v— L(p—<a). (66) 

Nous allons d’abord considérer le mouvement 
à une dimension de x suivant OX en partant de 
l’équation générale de Sturm-Lliouville (26). En 
même temps nous considérerons le deuxième cas 
possible (A,) pour le coefficient b, indiqué après 
l'équation (20), à savoir b, = 14, (a, réel) et qui 
correspond tout à fait au cas du champ +, A. 

Puisque dans le cas où ® — 0, À — 0 on doit 
retomber sur l’équation trouvée (47), nous com- 
plèterons maintenant l’équation (26) d’un membre 
supplémentaire &,, de façon que l’équation cher- 
chée soit de la forme plus générale. 


RS 


L 14 
dx? 1 dx 


+ Ur Are (67) 


A2 

Nous chercherons les fonctions inconnues &4, & 
dépendant sûrement des ©, 4, à l’aide du principe 
de correspondance. Pour y arriver, effectuons 
en (67) la substitution 


(68) 
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qui transforme (67) en 


dr Dire De DES 
dx? EE (E U5) + &e ag iÙ dE ra] RE=N0! 


(69) 

Puisque les fonctions f et Æ d’après (68) 
s’annulent simultanément, les équations (67) et (69) 
auront les mêmes valeurs propres. Comme nous 
avons dit pour l'équation (47) (et de même 
pour (67), resp. (69)), la densité relative des valeurs 
propres tend vers zéro si |[Æ,] — co c’est-à-dire 
(En — En_) [En — 0 si |E,]— co, ce qui mène à une 
distribution uniforme des niveaux d'énergie. Cela 
est vrai, dans le cas de l’équation (47), aussi pour 
des valeurs finies |[Æ,| si en même temps # — 0. 
Il faut naturellement exiger que cette propriété 
existe aussi pour l'équation (67), resp. (69)), c’est-à- 
dire que pour À —> 0 on retombe sur le cas clas- 
sique. Pour cela il est nécessaire que, dans l’équa- 
tion (69), k? soit diviseur de tous les termes dans 
la parenthèse — dans le cas contraire on aurait eu 
la propriété mentionnée pour À 0 indépen- 
damment des valeurs des 4,, &,, c’est-à-dire des o, 
A. Il ne peut sûrement pas en être ainsi pour des 
raisons physiques. Voilà pourquoi nous rempla- 
cerons dans l’équation (69) les nombres &, 4ë, &, 
resp. par @&/h?, aÿ[h?, af (cette dernière pour 
raison de dimensions). En revenant à la fonction 
à l’aide de (68) on aura, au lieu de (69), 


GR VC; 1 
dz? ñ dx i ñ? 


= 2m(E — U, + a,]f = 0. (70) 


Faisons dans (70) la substitution f — exp (18 /h) 
semblable à (31), et dans le résultat remplaçons la 
fonction S par le développement S —Y(i})"S,, 


n 


en s’arrêtant aux deux premières approxi- 
mations 55, S$,. On arrive ainsi à deux égalités, 
d’où l’on tire facilement 


ls 
Rs Ge aÿ | 
So 3 + \ 2 + 2m(E — U;) + (71) 
LA 2 50 
S1 = 280 La (72) 


d’où l’on aura S,(x) et S,(x) par intégration. Alors 
la dernière substitution pour f donne 


f = exp (iS01% — S;), [fl = w = exp ( 


29,). (73) 


En utilisant les relations w — C,/v (C, est une 
constante) et (66) on aura 


C,Vm 


e—281 = ———  , 
VERT 


(74) 


L'égalité (74) sera satisfaite si d’un côté le numé- 
rateur S, en (72) est remplacé par S5 + a,/2 (cela 
est équivalent à l’addition du terme ca, f/2h à 
la partie gauche de (70)) et si d’un autre côté la 
fonction inconnue 4, dans S% (71) est remplacée 


N° #4 


par — a /4. En effectuant les opérations indiquées, 
on tirera de (73), (71), (72) (C, constante) 


Cy 


f(x) = UE — Us (75) 
éxe GE f (+ VImlE — Ui)) del]. 
Les formules (75), (66), (71) donnent 
et) 
Pen (76) 
mo = VomE — Us) = ps —< Ac See 


Si 4, — 0 on doit avoir a, — 0 (l’équation (70) 
doit devenir l'équation (47)), de façon que S4 = Ps 
c’est-à-dire a, — — 2eA,/c. En introduisant les 
valeurs trouvées de a,, a, en (70) et en effectuant 
des petites transformations on représente faci- 
lement (70) sous la forme 


d 1e 2 Im 
Ge 254) FF EUR ONE 
qui est l’équation probabiliste cherchée — c’est 


l’équation d'amplitude connue de Schrôdinger en 
cas de présence des potentiels électromagnétiques. 
Sa généralisation pour les cas de deux ou de trois 
dimensions s'effectue immédiatement de la façon 
dont on a trouvé l’équation (57) en partant de (47), 
et de même pour le cas où les potentiels dépendent 
du temps. Dans ce dernier cas, pour les mêmes 
raisons qui ont servi pour écrire l’équation (60), la 
partie droite de l’équation cherchée doit être iden- 
tique à la partie droite de (60). On arrive ainsi à 
l’équation 


æ,U,% ch 
2m 2im dY 
AE ce en 
5, CAS DE PLUSIEURS CORPUSCULES. — Consi- 


dérons plusieurs microcorpuscules M;(x;, y, &) — 
M(r;) (@& = 1, 2, …, p) en mouvement dans un 
champ donné si l’on connaît aussi leurs interactions, 
c’est-à-direqu’on a lechamp U(x;, …,r,)(parexemple 
des électrons dans un atome). D’après l’hypothèse 
d'existence dans le subvac des formations ® et des 
fluctuations, le problème général quantique ici aussi 
sera de trouver la probabilité de distribution dw 
et l’énergie correspondante d’abord pour un état 
stationnaire À; d'énergie Ey, 

GA à 


D 
ee Ole II dr; - 
i=1 


p 
L dr; dy; dz; 
i= j=1 


E (79) 


et puis pour des états de mouvement quelconque. 
En partant de w on trouvera, comme plus haut, la 


N° 4 


fonction f, w — |f[2. Puisque les mouvements des 
corpuscules en général ne sont pas indépendants on 
doit prévoir que la fonction f(r,,.….,r,)satisfera à une 
équation unique (au lieu d’avoir des fonctions indi- 
viduelles f;(r:) satisfaisant à des équations indivi- 
duelles). Si les corpuscules sont sans interactions, 


DATE 


vements sont indépendants, il est clair que la den- 
sité de probabilité w doit prendre la forme 


D 
= IT w(r). 
if 


? 
Tr») = À Ui(ri), c’est-à-dire si leurs mou- 
1 


Alors l'équation pour f doit se décomposer en p 
équations individuelles, et l’on doit avoir 


mo) 


où les /; sont des solutions des équations suivant 
lesquelles se décompose l’équation initiale. Bien 
entendu au cas où les corpuscules se réduisent à 
un seul, cette équation doit se transformer en 
l’équation (57). Alors il est facile de voir que pour 
une telle équation généralisée on doit prendre 
l’équation suivante 

Ja 2 


A; f + 


i=1 Pi ñ° 


(80) 


Si Æ, est une valeur propre de (80) à laquelle 
correspond la fonction propre f, on peut attacher 
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à fx un facteur représentant une fonction pério- 
dique par rapport au temps, et l’on aura ainsi une 
fonction W,, comme (59), 


Polris ce, re à) = falri, , rp) TPM Er = hwy, (81) 


RENE 


où Ÿ, satisfait à l'équation suivante correspondant 
à l'équation pour un corpuscule (60) (en omettant 
l'indice k) 


que l’on admet comme valable aussi dans le cas où 
le potentiel U dépend du temps. (C’est l’équation 
de Schrüdinger pour p corpuscules.) 


Remarque. — L’équation (82) ou (80) appliquée 
pour le cas de p électrons dans un atome donne en 
première approximation la solution connue menant 
aux termes « d'échange » que l’on interprète formel- 
lement par échange des positions des corpuscules. 
Puisque la définition de la probabilité æ (79) ici 
est toujours de la forme de (1), on peut accepter 
que le fait mentionné soit en réalité lié à l’échange 
successif au cours du temps des positions des cor- 
puscules et que la fonction w mesure l’effet moyen 
de cet échange. 


Manuscrit reçu le 9 décembre 1959. 
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ÉTUDE DE L'AMORTISSEMENT PAR RAYONNEMENT COHÉRENT 
EN RÉSONANCE MAGNÉTIQUE NUCLÉAIRE 


Par HEenri.BENOIT, 


Laboratoire d’Électronique et de. Radioélectricité de la Faculté des Sciences, B. P. n° 9, Fontenay-aux-Roses. 


Résumé. -— On calcule la forme de raie d'absorption ou d'émission en R. M. N. lorsque l’amor- 
tissement par rayonnement cohérent est notable. Mais contrairement à la situation habituelle, 
les passages par la résonance en balayant le champ magnétique directeur ou la fréquence d’exci- 
tation ne sont pas équivalents et sont étudiés séparément. En particulier, le second mode de 
balayage met en évidence la bande passante de l’amplification du type Maser lorsque l’on observe 
l'émission stimulée. Les prévisions théoriques sont comparées avec les résultats expérimentaux | 
obtenus dans des champs directeurs de l’ordre du gauss. Les mêmes équations ont aussi permis de | 
calculer la plage d’accrochage et l’entraînement de fréquence lorsque le système fonctionne en 
auto-oscillateur du type « maser ». 


Abstract. — First the author calculates the shape of absorption or emission lines in N. M. R. 
when the damping by coherent radiation is considerable. But here, contrary to the usual con- 
dition, passages by resonance sweeping the magnetic director field or frequency are not equiva- 
lent and are studied separately. Especially the second sweeping mode clearly points out the 
band with the Maser type amplification when we observe the stimulated emission of radiation. 
The theoretical predictions are compared withexperimental results obtained a in magnetic direc- 


tor field of about one gauss. 


With the same equations, it is possible to calculate the circuit frequency 


range of oscillation and pulling when the system works as an auto-oscillator of Maser type. 


I. Introduction. — La première étude du cou- 
plage entre les moments magnétiques nucléaires et 
le circuit oscillant est due à N. Bloembergen et 
R. V. Pound [1] qui ont montré que l’amortis- 
sement par rayonnement cohérent contribue à la 
disparition des composantes transversales du 
moment magnétique dans les expériences de pré- 
cession libre, contribution caractérisée par un 
temps Tr = (2ryn M,Q) 1. Les expériences ont 
été faites par A. Skoze et S. Meiboom [2]. Mais 
contrairement aux conclusions de Bloembergen et 
Pound cet amortissement contribue aussi à la lar- 
geur de raie dans les expériences de passage lent : 
la théorie en a été faite par Bruce [3] puis par 
Bloom [4] qui en donne un traitement général. 
Nous reprendrons ici les calculs de Bruce en les 
adaptant au mode de détection et de balayage que 
nous utilisons et décrirons les expériences de pas- 
sage lent correspondantes. 


IT. Dispositif expérimental [5]. — Nous étudions 
la détection de la R. M. N. par la méthode du 
€ Q mètre » [6] : un générateur de f. é. m. e cos wt 
attaque un circuit oscillant dont la bobine B est 
traversée par un courant de benzène où se produit 
la résonance magnétique des protons, et nous mesu- 
rons Ja tension Æ aux bornes de cette bobine en 
l’envoyant après amplification sur la déviation 
verticale d’un oscillographe (fig. 1). 

La f. 6. m. e cos ot est obtenue en couplant très 
faiblement la bobine B par mutuelle induction avec 
une autre bobine alimentée à courant constant par 
un générateur basse fréquence. Un dispositif multi- 


phcateur de Q permet de régler la surtension du 
circuit oscillant. 

La bobine B est située au centre d’un système de 
bobines de Helmoltz qui créent un champ magné- 
tique À de l’ordre du gauss. Le courant de benzène 


€ cos wt 


NS 


HEC 


qui la traverse séjourne préalablement dans l’en- 
trefer d’un électroaimant un temps suffisant pour 
y attemdre la polarisation d'équilibre. La durée de 
trajet entre l’électroaimant et la bobine B est suffi- 
samment courte pour que la dépolarisation par 
relaxation soit faible de sorte que nous bénéficions 
de la polarisation obtenue dans l’électroaimant. 
Une bobine B' placée dans le champ de fuite de 
l’électroaimant permet de retourner le moment 
magnétique nucléaire par passage adiabatique en 
champ inhomogène |5] permettant d’observer 
l'émission stimulée au lieu de l’absorption (fig. 2). 

L'observation de la raie de résonance peut se 
faire soit en variant le champ magnétique 7 : c’est 
la technique habituelle de balayage en champ, soit 


N°4 


en variant la fréquence w /2x du générateur B. F. : 
on étudie ainsi la bande passante du système. Ces 
différents balayages se font en plusieurs secondes 
de manière à réaliser les conditions de balayage 


Pompe 


He 20000 G 


B’ 


| Electro-aimant | 


F1G. 2.— Schéma de principe de la circulation de benzène. 
La bobine B est la bobine du circuit oscillant de la 
figure 1. La bobine B’ permet de retourner éventuel- 
lement les moments magnétiques par passage adia- 
batique rapide en champ inhomogène. 


lent, ils commandent aussi l’entrée horizontale de 
l’oscillographe grâce à un amplificateur continu. 


IIT. Mise en équation. — Calculons la tension Æ 
aux bornes de la bobine B (fig. 1). Le courant 1 est 
donné par l’équation : 


te il 
e + i(iLo + Cle 0 


la self Z valant : 
L = Lifi + &rn(x — jx°)] 


x! et x” étant respectivement les susceptibilités 
correspondant à la dispersion et à l’absorption 
introduites par la R. M. N. et n le coefficient de 
remplissage de la bobine. 

Appelons w,1a pulsation de résonance du cireuit 
électrique (L, Co 5 — 1), Q sa surtension. En sup- 
posant Q élevé et w voisin de &w,, d’après (1) et (2) 
la tension £ — Lui est donnée par l’équation : 

+ PTE arnx | Cp | ÉS arnx| = 0. (3) 


Les susceptibilités y’ et y”, lorsque la saturation 
est négligeable sont données par les formules de 
Bloch : 


roi Q — © ; 
/ oo: + (Q — o)° Ti (a) 
72 1 1 F5 
nono To) Ts (5) 


en appelant Q la vitesse angulaire de Larmoor des 
spins nucléaires, Q = y, y le rapport gyro- 
magnétique, 7’, le temps de relaxation transversal 
et M, le moment macroscopique par unité de 
volume. 
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En introduisant les variables réduites suivantes : 
Q = 1/27yn MT; 
x = (Q— o)T, (6) 
B — (© — «;)/o 


l’équation (3) devient : 


 . Fe nc | 2e] k Je rl | à a. 
Da TDR Ce Me NS Er Er | 
SOLUTION APPROCHÉE. — On donne habituel- 


lement une solution approchée de cette équation en 
supposant 8 — 0 et Q < Q,. La variation rela- 
tive AE JE de la tension au moment du passage 
par la résonance est alors : 


AËJE = — hrny"Q 


et l'appareil est sensible au signal d'absorption. 

Nous étudierons ici la solution de l’équation (7) 
en abandonnant les hypothèses simplificatrices pré- 
cédentes et comparerons les résultats aux expé- 
rlences correspondantes soit lorsqu'on balaye en 
champ magnétique faisant varier alors Q, soit 
lorsqu'on balaye en fréquence, c’est alors & qui 
varie dans les équations (6) et (7). 


IV. Forme de la raie lorsqu’on balaye en faisant 
varier le champ magnétique. — Nous laissons la 
fréquence w/27r du générateur constante et égale 


LU 
S1 
a 
[a 
Q 
&=" À HT, 
2 L > 
0 1 2 3 & 
F1G. 3. — Forme de la raie lorsqu'on fait varier le champ 


magnétique À, pour différentes valeurs de Q[Q; notées 
sur chaque courbe. © est la surtension du circuit oscil- 
lant et Q — 1/27yn M, T, la surtension limite d’accro- 
chage. La mi-largeur à mi-hauteur est indiquée sur 
chaque courbe par un trait vertical. 


à la fréquence résonance w,/2x du circuit oscillant. 
Nous avons donc © — «4 soit 8 — 0. La valeur 
de « est alors « (Oo TE = AA en 
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appelant AH le champ de balayage. L’équation (7) 
devient alors : 


mor) +48) 0 


e 


La figure 3 représente (Æ) en fonction de & pour 


différentes valeurs de Q/Q,, la valeur de |e[Q étant 
prise arbitrairement égale à 1. 

Nous avons tracé aussi les courbes correspon- 
dantes à Q, < 0 c’est-à-dire à M, négatif : on 
observe alors l'émission stimulée. 

On constate que la mi-largeur à mi-hauteur des 
courbes d'absorption est supérieure à { (en valeur 
relative), tandis qu’elle est inférieure dans le cas de 
l'émission stimulée. Pour 10/0, | l’augmen- 
tation relative e de la largeur de raie est donnée 
au deuxième ordre par la formule : 


3 0 3 (£ 2 
LAC <) (9) 
n O1 32 \O: 


L’amplitude de la tension Æ lors du passage par 


[a] 


la résonance est proportionnelle à =, 
| | | + O0 

résulte une dissymétrie notable entre les signaux 

d'absorption et d'émission respectivement propor- 


- | 1 1 
tionnels à — 


D HQIQ 1 TOI 


b 


F16. 4, — Signal d'absorption (a) et signal d'émission (b) 


pour Q — 260 et Q — 110. La largeur du balayage en 
champ magnétique est de 33 mG. 


No 4 


L'accord entre ces prévisions théoriques et les 
résultats expérimentaux représentés par les photo- 
graphies n°8 4 et 5 est satisfaisant. 


b 


Fic. 5. — Signal d'absorption (a) et signal d'émission (b) 


pour Q, — 260 et Q — 135. La largeur du balayage en 
champ magnétique est de 33 m G. 


V. Forme de la raie lorsqu’on balaye en faisant 
varier la fréquence du générateur. — Le champ 
magnétique À reste constant et réglé à une valeur 
telle que YH = Q — w,. La valeur de 8 


B = (o — %0)1© 
est alors : 
œ or 


o Ta : 2 Qn 


en posant Q, — &w7,/2, Q, est la surtension de la 
résonance magnétique nucléaire. L’équation (7) 
s’écrit alors : 


= 


LE CHER œ PA OR 1 Ce 

E | Qn | red G+arl-° 
(10) 

E en fonction de « ne peut plus être calculé en 


fonction du seul paramètre Q/Q, contrairement au 
paragraphe précédent. Nous n’avons pas entrepris 


No 4 


de tracer le réseau de courbes correspondant, mais 
nous étudierons cette équation dans les cas corres- 
pondants aux conditions effectivement réalisées. 
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trique ou au contraire une saillie dans le cas de 
l’émission stimulée : 


c’est l’amplification du type 


£ 
F1G. 6. — (a) Courbe de réponse du circuit oscillant. Nous 
en déduisons sa surtension Q — 110. 
(b) Signal d'absorption. 


(ce) Signal d'émission. La largeur du balayage en fré- 
quence est de 18,7 H, autour de 3 656 H,. 


Sur les oscillogrammes réprésentés figure 6, la 
condition Q € Q, est suffisamment réalisée pour 
que l’équation (8) se simplifie et devienne : 


E = Ecfi — 4rny'Q) 


où £, donné par : 
e w ] 


ON, 


est la tension aux bornes du circuit oscillant 
lorsqu'il n’y a pas d’absorption par R. M. N. 

On détecte donc le signal d’absorption dû à la 
résonance magnétique nucléaire qui crée une cre- 
vasse dans la courbe de réponse du circuit élec- 


—10 


b 


F1G. 7. — (a) Courbe de réponse du circuit oscillant, d’où 
Q = 760. 
(b) Signal d'absorption dans les mêmes conditions. La 
largeur du balayage en fréquence est de 18,7 Hz. 


Fic. 8. — Courbe calculée à partir de l’équation 10 et 
correspondant aux conditions expérimentales de la 
figure 7 : 


OÙ 60 0 On 1070 61 Qi — 260: 
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« Maser ». Nous pouvons alors déterminer sans 
correction la largeur à mi-hauteur, donc le temps de 
relaxation transversale 7, — 0,09 s. La valeur 
correspondante de Q, est : 


Oro 12 020: 

Pour une valeur Q/Q, — 2,9, condition effecti- 
vement réalisée dans les résultats représentés 
(fig. 7) nous avons calculé £ par l’équation 10. Le 
résultat de ce calcul est représenté par la figure 8, 
l'accord avec les résultats expérimentaux est bon 
et montre bien en particulier la diminution notable 
de la tension aux bornes du cireuit au moment dn 
passage par la résonance ; 1l y a par contre aug- 
mentation de la tension sur les «ailes » de la courbe. 
La largeur et la profondeur de cette crevasse vont 
d’ailleurs en augmentant avec la surtension du 
cireuit comme le montre les résultats décrits dans 
la référence 5. 


VI. Calcul de la plage d’acerochage et de l’en- 
traînement de fréquence d’un Maser. — L’équa- 
tion (7) permet de calculer la plage d'accrochage et 
l’entraînement de fréquence d’un Maser à R. M. N. 
Dans le cas d’un fonctionnement en Maser nous 
observons une tension Æ alors qu'il n’y a pas 
d’excitation (e — 0). II faut donc : 


äl œ 
De et 

(11) 
OR 
OPERA 


Ce système n’a de solution que pour Q, négatif 
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et Q € —- Qi, nous retrouvons donc la condition 
d'accrochage d’un Maser (réf. 8) : 
271 My l20P 2 
La solution des équations (11) est alors : 
D RSC 
(O7 2 0! Qn V Q 


: Ou (Q — 
On UC 


(12) 


a) on 6). (13) 

L’équation (13) permet de calculer la fréquence 
d’oscillation w/27 en fonction de la fréquence de 
résonance w,/2r du circuit oscillant et de la fré- 


quence Q7/27 de la R. M. N. Écrite sous la forme : 


ep 2 ol Q) 
elle représente l'entraînement de fréquence. 

L’équation (12) done la plaque d’accrochage, 
désaccord maximum entre Q et w,, tel que l’on 
observe quand même l’auto-oscillation. Nous re- 
trouvons effectivement les résultats déjà établis par 
I. Solomon [9]. La comparaison des résultats 
expérimentaux avec les formules précédentes a été 
faite dans la référence [5]. 

Les équations précédentes ne permettent pas de 
calculer l'amplitude effective des oscillations : pour 
y arriver il faudrait tenir compte des termes non 
linéaires correspondants à la saturation que nous 
avons explicitement négligés dans les équations (4) 
et (5) donnant y’ et y”. Ce calcul compte tenu du 
mouvement dans la bobine a été fait par une autre 
méthode dans la référence [5]. 


(14) 


Manuscrit reçu le 16 janvier 1960. 
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ACTION DE L'OPÉRATEUR DE COLLISION ÉLASTIQUE DE BOLTZMANN 
SUR UNE FONCTION ISOTROPE DES VITESSES, DANS UN GAZ DE LORENTZ IMPARFAIT 


Par E. MOREAU et J. SALMON, 


Commissariat à l'Énergie Atomique 


Résumé. — On explicite l'opérateur de collision élastique de Boltzmann pour la partie isotrope 
d’un gaz de Lorentz imparfait. On retrouve la formule de Chapman et Cowling et l’on en précise 
les domaines de validité, on propose d’autres expressions s'appliquant dans les domaines où cette 


formule s’avère en défaut. 


Abstract. — The elastic collision operator is explained for the isotropic part of an imperfect 


Lorentzian gas. 


The formula of Chapman and Cowling is reobtained and the fields in which it is valid are spe- 


cified. 


Other expressions are proposed for fields where this formula proves to be inapplicable, 


Introduction. — L'utilisation de l'opérateur de 
collision élastique de Boltzmann se révèle dans la 
plupart des applications extrêmement ardue. Or le 
développement de la théorie cinétique des gaz 
donne une importance croissante à cet opérateur. 
Nous chercherons, ici, à l’expliciter, pour une fonc- 
tion isotrope, et un gaz de Lorentz. Nous retrou- 
. verons la formule de Chapman et Cowling, nous en 
discuterons les conditions de validité et, nous pro- 
poserons d’autres expressions s’appliquant au do- 
maine, pour lequel cette formule n’est plus satis- 
faisante. 

Un gaz de particules lourdes de masse M est en 
équilibre thermique. Sa densité est NV, et sa fonc- 
ton de distribution est maxwellienne, à la tempé- 
rature du milieu. Désignant cette fonction par F(V) 
on peut écrire : 


3/9 MV? 
AU M es > 2T 
F(F) » (= ) € (1) 
Avec : K': constante de Boltzmann, 


V : vecteur vitesse d’une particule Y, 
V : module de cette vitesse. 

À ce gaz est mélangé un gaz de particules légères, 
de masse m, et de concentration nr. On suppose 
n € N. Dans ces conditions la fonction de distri- 
bution f(&) de ces particules ne dépend que des 
collisions avec les particules AZ. 

Désignons alors, dans une collision de ce type, 
les vitesses avant et après choc par ® et v’ pour la 
particule m et par V et V’ pour la particule #. 
g et g' représentent les vitesses relatives corres- 
pondantes : 

ENT EUE (2) 
g' = v — V!. 


Les lois de conservation des quantités de mou- 
vement et d'énergie cinétique s’écrivent : 

mo + MV = mv + MV’ (3) 

mu? + MV?— mv°? + MV'?. (4) 


D'où l’on tire aisément : 
BB (5) 
Supposant la loi d'interaction centrale, et con- 
naissant celle-ci, on peut en déduire la section effi- 
cace différentiellec(0, g) danslesystème du centrede 
masse. 0 désigne l’angle des vecteurs g et g’. Nous 
définirons également l’angle y comme l’angle de la 
projection de g’ sur un plan perpendiculaire à g 
avec une direction arbitraire de ce plan. 


Les paramètres 0 et n avec les équations (3) et (4) 
déterminent complètement v' et V’ en fonction dev 
et y. 

Ainsi l’opérateur de collision élastique de Boltz 
mann s'écrit : 


2TT T 
Te ane WE’ fF) g o(6, g) dnsin 0 d6 dV. 
(2) 
On a désigné par f, f', F, F',les quantités j(v), 
f(v'), F(v), F(v'), et l’on a étendu l’intégration sur 
la totalité de l’espace des vitesses P dont dV 
désigne l’élément différentiel. 


Action de / sur une fonction isotrope /(v) : 
Posons : 


sl) 
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Avec : 


NS Er e : : 
La fonction de distribution des molécules restant 
maxwellienne, nous avons : 
fo F° — fo F) = 0 (10) 
puisque : 
mv? + MV? = mv? + MV’? (11) 
donc : 
LE—CRF = jo @ 8" — fo of = fo Flo — +) (12) 
où l’on a désigné par +’ et o les quantités œ(v') 
et p(v). 
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= D.- -.(1 + 52 + 2s cos 6)12 


% 


N° & 


L'opérateur de collision, dans ces conditions, 
s'écrit : 


= marrante fe, A, LE LE, 


(@'— p)e— MP g 6(0, g) dn sin 0 dO sin 6 dé dy V? dy. 


Les angles &, y désignant la colatitude et la 
longitude de V. 


Caleul de la vitesse v’. — Ce calcul est effectué 
en appendice I. Il montre-que le module de v' ne 
dépend que de 0, n, B, V et v. Ce module s’exprime 
par : 


L y? 1e 
É ji 2 e [(s — 1) cos B (1 — cos 0) + sin B sin 0 sin n (s + 1)] + 2 2 (1 — cos | (4) 


1 & s2 + 9s cos Ô 


si l’on a posé : m]M =. 

Si l’on porte v’ tiré de (14) dans (13), on obtient 
une expression de /(f) correcte, mais imextricable. 
Nous chercherons des approximations ayant un 
sens physique, en divisant les vitesses v en diffé- 
rents groupes à l’intérieur desquels les approxi- 
mations seront d’un type particulier. 

Nous distinguerons le groupe des vitesses éle- 
vées, de loin le plus important, et le groupe des 
vitesses très faibles. 

L'influence du groupe des vitesses moyennes 
pourra le plus souvent être négligée. 


A) GROUPE DES VITESSES ÉLEVÉES. — On dé- 
signe par là les vitesses v pour lesquelles l’inégalité 


(V2) Jo € 1 (15) 
ou ce qui est identique : 
8kT ]Mv? & 1 (16) 


est satisfaite. V? représentant la vitesse quadra- 
tique moyenne des particules M. 

Ces inégalités signifient que pour la quasi tota- 
lité du spectre des vitesses des particules M l’iné- 
galité : 


(Vo) <1 (17) 
est vérifiée. On posera par conséquent : 
(18) 


WIDE 


+ 2s cos 0 + s2)12 


et, on considérera cette quantité comme petite 
devant l’unité dans l’expression de J. 
Nous avons déjà posé : 


SUN ne PHOITAEEQUE (19) 


Les quantités s!/? et x’ sont du même ordre de 
grandeur pour un mélange de deux gaz en équilibre 
thermique, en l’absence de forces extérieures. Dans 
ce cas, en effet, la distribution des particules m 
est isotrope et maxwellienne 


ere ( te le — mo /2KT 
— —= nn \- (9 Ë 
Q Toul 


Et celle des particules M étant : 
F = N(MI2r KT)? e—MY°PRT, 


On déduit : 


(20) 


(21) 


V?[0? = m/M. (22) 


La quantité V/v — «x! est donc, en moyenne de 
l’ordre de grandeur de st”, 

En présence d’une force extérieure accélératrice, 
le rapport des vitesses quadratiques moyennes 
sera plus petit, le gaz de particules m possédant une 
température cinétique plus élevée que T. 


Approximation de J(f) : 
Considérant &’ et s/? comme des quantités du 
même ordre, et nous limitant au développement 


du 2° ordre, nous obtenons en tenant compte 
de (14) 


[ + 2a’[(s — 1) cos B(1 — cos 0) + (s + 4) sin 6 sin n sin 0] + 2æ&/2{1 — cos I (23) 


1 + 2s cos 8 + s2 
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Et: 


v'œv [1 + &” (sin 6 sin 0 sin n — cos B(1 — cos 8)) 


2 
Z + &'2(1 — cos 0) — É {sin & sin 0 sin n — cos B{1 —- cos 0)]? — sft — cos o)|. (21) 


= 


D'où les expressions : 


Ce . : : 
— : [sin 8 sin 0 sin n» — cos B(1 — cos 0)]? — s(1 — cos o)] (25) 
(v' — v}? & x’? v? (sin 8 sin n sin 0 — cos B{1 — cos 0))? (26) 
La quantité v' — vest donc, dans ces conditions, Dans le même esprit, 


une quantité petite vis-à-vis de v et nous pouvons 
considérer o(v') comme peu différent de o(v). 
Ceci permet d’obtenir une approximation de 


g=e—_V (29) 


g = (0? + V?— 2vV cos B)"? 


@(v") — o(v) par un développement Taylorien, soit : ë 
/ PS ET re. de v'—v2ide — 01 — 2° cos 8 + a/?2)v2 (30) 
go gb= tp ne+( =") Se (7 : 
: CHI 
ou avec (25) e (26) ï op (1 — 7 COS B + 7 n° ) 


t 
ÉMIS Ê Penn con pe D) et comme g diffère peu de », 


! 0 2e 9 42 
= Cos oO) la s) se AE do(v, 0) | (g— v)? d?o(v, 0) 
x'2 de o(g, ) + o(v, ) Î (g v) do Ï 2 dv? 

—— (sin $ sin 9 sin 0 — cos F1 — cos ojj®] ie 
È (31) 
LE d? o 
—— (sin B sin n sin 6 — cos (1 — cos 0))? 0? —+. EE : : : ; 
IR SRRS pt TE Réunissant les approximations (28), (30, (31), 


(8) il vient pour J(f) : 


o(v)) eMP2ET 96(g, 0) dn sin 0 d0 sin B df V?dV (32) 


DEN M 3/2 CO T TT 2TT s 
Aer Cr) fo RE ete) 


l'intégration sur y étant immédiate. L'intégration sur 1 est élémentaire car seul œ(v') — o(v) est fonc- 
tion de n d’où : 


J(f) | " Her N 


27510 
12 


: sin? 6 sin? 0 , Sarl 
le 1 hs [Le — s) (4 — cos 0) — «(1 — cos 0) cos Ê — = JAN EE + cos? 6 (1 — cos o2|| DE 


Open ciné 2 «+ 2 : s r9 r 99 
e EE ES + cose p{1 = cos 0°] p? al MP golg, 0) sin 0 dO sin B dB V?dV. (33) 
- ne 


À pA 


En tenant compte de (31) et (30), nous obtiendrons, une approximation de o(6, £) : 


CCR don tu eD Re AA: x 
o(0, g) = 0(0, v) E «'v (£ sin? B — cos 4) . 7 COS” Be: (34) 
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Dans J(f), les produits effectués, en négligeant les termes supérieurs à «'? xs? et s,nous amènent à : 
M ee en [2 > 49 ni I | ,9 P \ , 
œ &4T? 2 — 5) (1 — cos 0) — (1 — cos Ô) cos Ê 
ere) HN Lee (x? — s) ( 
rene eine | 
= F _ . L cos? B(1 — cos 6) ] + x’2{1 — cos 0) cos? | 
do x?fsin?Bsin20. ; 
+ D vo? cos? B(1 — cos o)] L- vè Pr [ 2 L cos? B(1 — cos 0)?||. (35) 
e—MV#2kT sin 0 46 sin £ d£ V?d”. 
L'intégration sur 8 donne ainsi : 
M \3%/2 
= 
J( = 8% (= %) fo N 
: 3 3 42 A : | 
“ Ne [ea — cos Ô) (v° : 6 — = Fe É ee = me 5) 27 o(1 cos o] e—MV®#/2KkT sin 6 d8 V2? dV. (36) 
V-0 /6-0 
Rappelons que : donc : S 
= V Jo Ci: (18) do SM (: mm ] "orme (É … MU ] (44) 
Re do — 27 #1 do KT 
AT n(— Fr) eMVIERT QV — 1 (37 : ee. 
Doll ) d° 9 it mm Fe +2 /2kT (45) 
du? —_ (2x aT a 
M \#2 f 
72 —MVIXT AV — 3kTIM.. (: 2 
a kr V? _ ) e—MV?RkT AV = 3kT]M.o (38) É f . 2mv de (1 
N— ’ À AE do NN ETAT TETE : 
L'intégration sur V est immédiate : 
Et finalement : 
J(ÿ =} [SE (eer dé v?) d KT df 
TZ Jo [Oo ——, m ; k 16 
é © Aie Mo? a (° ss (+ 7 mv LL pi 


N o{v, 0) (1 — cos 0) sin 0 11] 


T 
Î 2 
0 


" Le et 
M ‘dv 


KP QËe 


dolv. 0) 
0) 


27 Nv (1 — cos 0) sin 0 del 


F 27 Nov, 0) (1 — cos 0) sin 0 do]. 
0 


(M av? 
(39) 
La quantité : 
— ne 2r volv, 0) (1— cos 0) sin 040 (40) 


appelée « fréquence de relaxation d'ordre 1 » part 
MM. Bayet, Delecroix et Denisse nous permettra, 
en remarquant que : 


VU N [lex (1 — cos 6) sin 0 dO (41) 

oi 7 = N for ue (1 — cos 0) sin 0 d0 (42) 
d'écrire (39) sous la forme : 

na = Gr (en 2) mn) sua 

(43) 


Mais nous avions posé : cf. (8) et (9) : 
= fo p = neMRET (m]27r kT)3/2.œ 


C’est la formule de Chapman et Cowling. 

À l’équilibre f —f,, ® —1 et J(f) s’annule: 
Cette formule, malheureusement, présente des 1n- 
convénients, puisqu'elle n’a de sens que pour # € V 
En particulier, lorsqu’on désire calculer une fonc- 
tion de distribution des particules m en résolvant 
par une méthode numérique l’équation de Boltz- 
mann, la présence d’un point singulier à l’origine 
introduit des difficultés. Physiquement la formule 
de Chapmann et Cowling est fausse pour les cas où, 
en moyenne «’ est supérieur ou égale à 1 c’est-à- 
dire pour v < (V?)%. Ceci, nous conduit à re- 
chercher une formule pour les vitesses faibles. 


B) GROUPE DES VITESSES FAIBLES. — Elles sont 
définies par : 


PETER (&7) 


Examinons, d’abord le cas v — 0. 

Vitessei=v— 0. 

Reprenant (14) et tenant compte du fait ques 
est petit devant l’unité, nous obtenons : 


2 TeSiniO® 
Ù = —— 


Dre 9 VEsin,0 f2 


(48) 
et de (30) : 
1 


o 
(e] 


N°4 


d’où une expression de J(f) : 


L'OPÉRATEUR DE COLLISION ÉLASTIQUE DE BOLTZMANN 221 


: M \3%2 fo fr _— n 
Je 2 N (= 5) si e—MPRKT|o(2V sin 0/2) — o(0)) VS o(V, 0) sin 0 dO dF. (49) 
Vitesse sv —\(. 
On serait tenté d’effectuer un développement en série du type : 
3 ue ef v L ù He 
DA NOP [1 + AV Sins 02 (sin 0 sin $ sin n(s + 1) + cos B(1 — cos 0) (s — 1))] 4 (50) 


Malheureusement ce développement n’est pas 
valable lorsque V sin 0/2 est très petit. Il convient 
donc d’effectuer des coupures. On choisira deux 
bornes 0,, et V,, telles que 

2Vnsin 0/2 > v. (51) 


Pour v' —2V sin 0/2 et V sin 0/2 très petit, 
v' tendra vers vet J(f) vers zéro, donc, en dessous 
des bornes ci-dessus nous supposerons J(f) —0 et 
nous aurons : 


M 


3/2 TT 
TJ =8h.N x) 44 oi e—MVPBET [e(2V sin 0/2) — ©(0)] V3 o(F, 6) sin 0 d dY 


en admettant que : 
DE 


Malgré l'inconvénient de l'opérateur inté- 
gral (52), cette approximation de J(f) aux vitesses 
faibles, devrait permettre de lever à la fois, les 
difficultés mathématiques rencontrées dans l’usage 
de la formule de Chapmann et Cowling à l’origine, 
et l’erreur physique commise en utilisant la même 
formule aux vitesses faibles. 


Appendice I : Caleul de la vitesse v’ après choc. 
— La fonction de distribution des particules m 
étant isotrope, nous prendrons comme système de 
coordonnées, un système de coordonnées carté- 
siennes où ® est porté par Oz d’où les notations : 


y, B angles polaires de V (avec Ox arbitraire) ; 
& colatitude de gi, équipollent à g passant par O ; 
n, Ô angles polaires de g’ dans un systeme 
d’axex Ox/y'z/ où gn et sur Oz' et, Ox’ dans le 
plan », V (rotation de y par rapport au 3e axe). 
_ Si nous représentons la matrice de rotation de 


. (52) 


angle y par rapport à l’axe ? par 6;, nous aurons 
pour les coordonnées de g’ dans Le système Oxyz: 


ge 
|| — 63,y : 02,a + O3,n - 02,0 | (3) 
82 8 
Sa “0 0 D 
|| = mr mal | = | 0 | + Le] (4) 
Bz _£ — V v 


La conservation de l’impulsion et d’énergie,nous 
amène à : 


(m +- 


M) v = mv + MV + Mg’. (5) 


Done, après avoir effectué (3) et (4) nous aurons : 


gx = — V sin Ê cos y = g sin & cos y 

gy = —- V sin B sin y = gsin «sin y (6) 
Dé — D — cos À — 12 CUS, 

(2èts - 
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D'où : 
ge — (v — V cos B) cos n’ sin 0 — V cos 0 sin 
gi — gsin n sin 0 : (7) 
ge — V sin B cos n’ sin 0 + (v — V cos «) cos 0. 
Et : 
("= =) Ve — (1 — V cos GB) cos n’ sin 0 cos y 
+ Fsin 8 (1 — cos 0) cos y — gsin r’ 10 Sin ÿ 
(5) vy — (v — V cos Ê) cos n° sin 0 sin y 


+ V sin 8 (1 — cos 0) sin y + g sin ny sin 0 cos y 


im + M\ : : ; ) 
(+) v; = V sin 6 cos »’ sin 0 + v(F - cos 1) 


Re 
+ V (1 — cos 0) cos £. (8) 


Le carré du module de &’ sera donc obtenu, en 
faisant la somme des carrés des 3 équations (8). 
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En remarquant que les 2 premières équations sont 
de la forme : 
A cos y — Bsin y et À sin y + B cos y 
nous voyons que la somme de leurs carrés est : 
A? + B? c'est-à-dire que le module de &' est indé- 
pendant de +. Nous prendrons 7 — 0 dans (8) et 
NOUS ArrIverons à : 
A 
(m + M)? 
+ 20 V((s + 1) sin 8 sin 0 cos n° + (s— 1) cos B(1 — cos 0)) ,| 
+ 2P2{1 — cos 0)) (9). 


(u2{1 + s? + 2s cos 0) 


OF 


(m + M)]M =1 +s s = m|M 


la formule (9) serait la formule (14) du texte si 
l’on remplaçait l’angle arbitraire 1! par son com- 
plément n. Cette substitution n’amène aucun chan- 
gement aux résultats, et n’est motivée que par des 
raisons mnémotechniques. 


avec 


Manuscrit reçu le 28 septembre 1959. 
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ÉTUDE EXPÉRIMENTALE 
DU SPECTRE DES ÉLECTRONS D'AUTOIONISATION EN RADIOACTIVITÉ B 


par F. SUZOR, 


Laboratoire de Physique Nucléaire, Orsay. 


Résumé. — Faisant suite à des expériences déjà publiées sur le %2P, 358, 147Pm, les résultats 
sont donnés pour le °°Y, 1#Pr et ?2Na. Contrairement aux intensités des raies X, Æ ou L, qui sont 
en bon accord avec la théorie de l’autoionisation, les spectres continus, mesurés entre 1 et 13 keV, 
des électrons d’autoionisation émis simultanément avec le rayonnement 8 sont en contradiction 
avec les prévisions théoriques (intensité plus grande, distribution spectrale différente et loi inversée 
en fonction de Z). Ceci s’expliquerait par une contribution des couches externes du cortège élec- 
tronique considérablement plus importante que celle prévue par la théorie de l’autoionisation. 


Abstract. — Following already published experiments on #?P, %5%$, 147Pm, results are given 
for °°Ÿ, 18Pr and Na. The intensities of the X, Æ or L lines, are in good agreement with the 
theory of autoionisation. On the other hand the continuous spectra, between 1 and 13 keV, of 
autoionisation electrons emitted simultaneously with 8 radiation disagree with theoretical predic- 
tions. (The intensity is greater, the spectral distribution is different, and the law of variation 
versus Z is inverted.) This could be explained by a contribution from external electronic shells 
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considerably more important than that predicted by the autoionisation theory. 


Introduction. — Lorsqu'un noyau radioactif se 
désintègre par émission B, le cortège électronique 
subit une perturbation le laissant parfois dans un 
état ionisé. La théorie de ce phénomène, appelé 
autoionisation, a été effectuée par Migdal [1] et 
Feinberg [2], puis par Levinger [3] et d’autres 
auteurs (voir pour une étude récente et une biblio- 
graphie détaillée la thèse de Grard [4]. Il existe 
trois moyens expérimentaux de vérifier ces théories : 
1) La mesure des raies X consécutives au réarran- 
gement de l’atome aprèsautoionisation; 2) la mesure 
de la charge de l’atome résiduel ; 3) la mesure du 
spectre continu des électrons émis lors de l’auto- 
ionisation et des raies Auger consécutives au 
réarrangement du cortège électronique. 

En dehors du cas de la radioactivité fB*, ce 
dernier type de mesure est impossible par les 
méthodes ordinaires de spectrographie magnétique 
et nécessite que soient enregistrés en coïncidence 
les électrons du spectre 8 et les électrons d’auto- 
ionisation. Les résultats expérimentaux ainsi obte- 
nus et exposés dans cet article, sont la suite de 
travaux commencés avec G. Charpak [5 et G]. 

Après une première expérience [5] obtenue avec 
un appareillage mal adapté au but recherché, nous 
avons construit [7] un spectromètre constitué de 
deux grands compteurs proportionnels accolés lun 
à l’autre. La source étant vue de chacun des comp- 
teurs sous un angle solide voisin de 27, l'excellent 
rendement ainsi obtenu rend possible, par la mé- 
thode des coïncidences, l’étude de phénomènes peu 
fréquents. Le fait que la source puisse, sans l’inter- 
médiaire d'aucun absorbant, être en contact direct 
avec le gaz de l’un des deux compteurs permet 
l'étude du spectre des électrons jusqu’à des éner- 


gies aussi basses que 1 keV environ. Enfin les 
grandes dimensions des compteurs entraînent pour 
les électrons énergiques du spectre 8 une perte 
d'énergie suffisante pour que la séparation entre 
eux et les électrons de plus faible énergie émis 
simultanément soit possible. L'étude de ce type de 
spectromètre, lés précautions à prendre et l’élimi- 
nation des coïncidences parasites ont déjà été expo- 
sées dans les publications précédentes [6 et 7]. 


Yttrium 90. — Ce radioélément se désintègre par 
émission f— et possède un spectre dont l’énergie 
limite est de 2,26 MeV. Les sources étaient obtenues 
par séparation sur résine échangeuse d’ions à partir 
d’une solution contenant en équilibre °°%Sr et °0Y ; 
une activité de l’ordre du centième de microcurie 
était déposée par évaporation thermique sous vide 
sur une membrane mince de formvar aluminisé 
d'environ 40 ug/em? et d’un diamètre de 5 mm. Les 
compteurs proportionnels étant remplis de pro- 
pane, une série de mesures a été faite sous une 
pression de 57 cm de mercure et une autre sous une 
pression de 8 em. La source était en contact direct 
avec le gaz du compteur 1 et les impulsions enre- 
gistrées dans celui-ci étaient, après amplification, 
comptées dans un sélecteur à 50 canaux ; dans le 
compteur 2, seules étaient enregistrées après ampli- 
fication et passage dans un sélecteur à seuil, les 
impulsions correspondant à une énergie perdue 
dans le compteur supérieure à 25 keV. Les coïnci- 
dences obtenues entre les deux compteurs four- 
nissent entre 4 et 13 keV environ le spectre des 
électrons émis simultanément avec le rayonne- 
ment 8. Les phénomènes parasites qui pourraient 
s'ajouter à ce spectre sont les suivants : 
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4) Les photons émis simultanément avec le 
rayonnement & ; leur importance relative est négli- 
geable, comme le montre d’ailleurs une expérience 
annexe obtenue en recouvrant la source avec un 
écran mince perméable aux photons et arrêtant les 
électrons 

2) Les coïncidences obtenues par un électron £ 
perdant une très petite partie de son énergie dans 
le compteur 1 et associé à un électron ou à un 
photon d’énergie supérieure à 25 keV et enregistré 
dans le compteur 2. Les dimensions des compteurs 
sont suffisantes pour que, avec un remplissage en 
propane sous une pression de 57 em de mercure, le 
spectre B ait un maximum vers 15 à 20 keV et 
fournisse très peu d’impulsions inférieures à 10 keV. 
Dans ces conditions ce type de coïncidences para- 
sites est absolument négligeable jusqu'à 10 ou 
12 keV ; par contre c’est lui qui interdit dans le 
cas présent l’exploitation des résultats pour des 
énergies supérieures 

3) Les coïncidences obtenues par un électron B 
ayant un parcours de quelques millimètres dans le 
compteur { et ensuite renvoyé dans le compteur 2 
par rétrodiffusion sur une molécule du gaz. Ce 
phénomène peut être étudié en faisant varier le 
diamètre de la membrane supportant la source qui 
constitue un diaphragme de séparation entre les 
deux compteurs et en modifiant la pression de rem- 
plissage des compteurs. Dans les conditions expéri- 
mentales utilisées, ce phénomène est négligeable. 

4) La projection par le rayonnement B d’élec- 
trons dans le support de la source ou dans la 


100 


Nombre d'électrons par keV pour 10 “rayons B 
on 
o 


Rice 


matière environnante. Ceci a déjà fait l’objet d’une 
étude détaillée [5] montrant que ce phénomène 
ne peut être responsable que d’une faible partie du 
spectre obtenu ; il est d’ailleurs d’autant moins 
important que l’énergie du rayonnement & est plus 
grande. Dans le cas présent le spectre de la figure 1 
a été obtenu après soustraction d’une correction 
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estimée [6]de ce phénomène dont l’importance rela- 
tive varie de 40 % à 30 % entre 5 et 13 keV et peut 
être pratiquement négligée en dessous de 5 keV. 
Le spectre des électrons représenté sur la figure 1 
et associé au rayonnement B de l’yttrium 90 se 
décompose en un fond continu décroissant avec 


l'énergie et un pic correspondant aux électrons 


Auger de la couche L. Grâce au sélecteur à 50 ca- 
naux utilisé et à la durée de chaque mesure excé- 
dant plusieurs heures, les fluctuations statistiques 
ne sont pas gênantes pour la précision des mesures, 
celles-ci étant surtout limitées par la stabilité de 
l’'amplification du compteur proportionnel et par le 
recouvrement des séries de mesures faites avec des 
valeurs différentes d'amplification. Cette précision 
est de l’ordre de 5 % au début du spectre à 1 keV et 
décroit jusqu’à environ 20 ou 30 % à la limite du 
spectre vers 13 keV. 

Une prévision du spectre continu, tirée de la 
théorie de autoionisation a été tracée en pointillé 
sur la figure 1. On voit que l’accord n’est pas satis- 
faisant comme le montrent les chiffres du tableau 
suivant. 


ONE de 1 à de 3 à de 8 à 
3 keV 8 keV 13 keV 
Spectre continu expé- 
rimental JA 38.104 C0 
Théorie de l’autoioni- 
sation 16 10 14.10—4 3 40 


La courbe théorique de 1 à 13 keV sur la figure 1 
a été obtenue en utilisant, outre la forme théorique 
du spectre pour les couches X et Z une probabilité 
d’autoionisation de 0,9/Z? pour la couche X et une 
probabilité de 8/72? pour la couche Z, tenant ainsi 
compte de l’effet d’écran suivant les estimations de 
Levinger ; Ja théorie prévoit dans ce domaine 
d'énergie une contribution négligeable pour les 
couches 4 et suivantes. 

Il n’a pas été tenu compte de la rétrodiffusion 
des électrons par la membrane de formvar suppor- 
tant la source. Même dans les conditions de géo- 
métrie utilisées, il ne semble pas que cet effet amène 
une surestimation du spectre supérieure à 30 % 
environ. Par contre tout effet d’autoabsorption 
entrainerait une sous-estimation du spectre ; étant 
donné la netteté du pic des électrons Auger de la 
couche Z, et la pureté du radioélément employé, il 
est légitime de penser que l’autoabsorption des 
électrons du spectre continu doit être faible. 


Praséodyme 143. — Ce radioélément se désin- 
tègre par émission f— et possède un spectre dont 
l’énergie limite est de 0,93 MeV. Les sources étaient 
obtenues comme dans le cas de l’yttrium par évapo- 
ration thermique sous vide sur une membrane 
mince à partir d’un radioélément pur en prove- 
nance d’Harwell, contenant, d’après la notice d’en- 
voi, 40 mg de matière par millicurie. Les compteurs 
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proportionnels étant remplis de propane sous une 
pression de 57 cm de mercure, les conditions expé- 
rimentales étaient par ailleurs semblables à celles 
décrites précédemment pour lyttrium — (les 
spectres de la figure 2 ont été obtenus après sous- 
traction d’une correction relativement plus faible 
que celle indiquée pour l’yttrium et correspondant 
à la projection par le rayonnement 8 d’électrons 
dans la matière environnant la source), — Sur la 
figure 2, le spectre 1 correspond à une source recou- 
vrant une surface de l’ordre du millimètre carré et 
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devant contenir environ 25 ug/em? de matière : le 
spectre 2 correspond à une source plus largement 
étalée et devant contenir environ 2 ug/em?; la 
quantité appréciable de matière contenue dans la 
source 1 est responsable d’une auto-absorption du 
spectre décroissant avec l’énergie et parfaitement 
compatible avec la loi d'absorption des électrons 
dans ce domaine d’énergie [8]. Il est même très pro- 
bable que, en dessous de 4 à 5 keV, le spectre 2 est 
encore inférieur au spectre réel par effet d’auto- 
absorption. L’intensité minimum du rayonnement 


Pr 143 
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nécessitée par l’expérience n’a cependant pas per- 
mis d’étaler plus la source ou de diminuer son 
intensité. 

Le spectre des électrons représenté sur la figure 2 
et associé au rayonnement $ du praséodyme 143, se 
décompose en un fond continu décroissant avec 
l'énergie et un pic correspondant aux électrons 
Auger de la couche Z. La prévision théorique 
tracée sur cette figure tient compte, comme dans le 
cas de l’yttrium, de l’effet d'écran suivant les esti- 
mations de Levinger et correspond pratiquement 
entièrement à une ionisation dans la couche Z, les 
contributions de la couche Æ et des couches #7 et 
suivantes étant négligeables dans ce domaine 
d’énergie. | 

Comme dans le cas de l’yttrium, la comparaison 
entre théorie et expérience est résumée dans le 
tableau suivant : 


ee de 1 à de 3 à de 8 à 
3 keV 8 keV 13 keV 
Spectre continu expé- 
ti source 2 DID DÉMOS 210 
Théorie de l’autoioni- ; 
sation C0 M sl 


Ces résultats ne tiennent pas compte de la rétro- 
diffusion des électrons par la membrane de formvar 
supportant la source. Par contre nous avons vu 
qu’étant donné la qualité insuffisante des sources, 


l’autoabsorption pourrait 
dessous de 4 à 5 keV. 
En plus de ces résultats 


être appréciable en 


relatifs aux électrons 


par minute 
11400 par minute 


Nombre de coincidences 
pour N; 


une expérience en tout point semblable à celle déjà 
publiée [6] pour le prometheum 147, a donné pour 
le praséodyme 143 l'intensité du pic à 5,2 keV des 
raies X-Z du néodyme obtenues par autoionisation. 
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Les compteurs étaient alors remplis d’un mélange 
d’argon sous une pression de 60 em de mercure et 
de propane sous 6 cm. Dansle compteur 1, vu de la 
source sous un angle solide 2x, les photons étaient 
enregistrés à travers un écran de 20 mg/cem? de bery- 
lium. Dans le compteur 2, vu de la source sous un 
angle solide légèrement inférieur à 2x, était enregis- 
tré lerayonnement 8 perdant une énergie supérieure 
à 12 keV (11 400 impulsions par minute). La figure3 
reproduit le spectre des impulsions en coïncidence 
dans les deux compteurs, obtenu dans un sélecteur 
à 50 canaux après soustraction du mouvement 
propre et des coïncidences fortuites connus avec 
précision. Le pic de 5,2 keV dont on veut mesurer 
l'intensité est superposé à un fond continu de frei- 
nage et à une faible raie de 8,0 keV qui est la 
raie X-X du cuivre, correspondant à une ionisation 
des parois du compteur par des photons de plus 
grande énergie ; la remontée vers les basses énergies 
correspond au flanc de la raie X-X de l’argon 
de 3,0 keV, provenant de l’ionisation par le rayon- 
nement 8 des atomes d’argon voisins de la source. 
La séparation du pic de 5,2 keV s’opère facilement 
et permet d’obtenir, compte tenu des angles solides 
et de l’efficacité du compteur, pour les raies X-L 
d’autoionisation, que le pouvoir de résolution du 
compteur ne permet pas de séparer, la valeur 
(3,5 + 0,4) 1074, En admettant pour le rendement 
dé fluorescence de la couche Z du néoydme la 
valeur 0,18 on trouve pour la probabilité d’auto- 
ionisation dans la couche Z la valeur (1,9 + 0,2)10—3 
en bon accord avec celle prévue par la théorie de 
Levinger qui est 8/2? — 2,2.107$,. 


Sodium 22. — Ce radioélément se désintègre par 
émission B* avec un spectre dont l’énergie maxi- 
mum est 0,54 MeV, suivie de l’émission d’un y de 
1,276 MeV. Les sources étaient obtenues comme 
précédemment par évaporation thermique sous vide 
et les compteurs proportionnels étant remplis par 
du propane sous une pression de 57 em de mercure, 
les conditions expérimentales étaient les mêmes que 
dans les cas précédents. La désintégration éven- 
tuelle du sodium 22 par capture électronique ou 
lannihilation des électrons positifs étaient prati- 
quement sans influence dans les conditions expéri- 
mentales choisies ; par contre l’existence du y émis 
à la suite du B* rend plus abondant un des types de 
coïncidences parasites déjà mentionné lors de 
l'étude de l’yttrium 90 ; il s’agit des coïncidences 
obtenues d’une part entre un électron enregistré 
dans le compteur 1 (spectre NW.) et d’autre part 
un y enregistré dans le compteur 2 ; dans ce der- 
nier on enregistre un nombre constant d’impulsions 
correspondant à une perte d’énergie supérieure 
à 25 keV. Le spectre parasite des coïncidences ainsi 
obtenues est proportionnel au spectre N, qui, faible 
en dessous de 10 et 12 keV commence ensuite à 
croître rapidement ; la correction de cet effet para- 
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site est obtenue en lui attribuant au delà de 20 keV 
la totalité du spectre des coïncidences : le coefli- 
cient de proportionnalité étant ainsi déterminé, on 
obtient par soustraction un spectre nul par défi- 
nition à partir de 20 keV ; de 13 à 20 keV l’arbi- 
traire de l’hypothèse admise rend les résultats 
incertains, par contre au-dessous de 13 keV, la 
correction effectuée est faible étant donné le mini- 
mum du spectre NV, dans ce domaine. Le spectre en 
coïncidences ainsi obtenu de 1,5 à 13 keV est repré- 
senté sur la figure 4. 
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Malheureusement la raie des électrons Auger de 
la couche Æ a une énergie trop faible pour appa- 
raître sur le spectre de la figure 4. La qualité du 
radioélément employé permet de penser que les 
sources utilisées contiennent moins de 0,5 ug/em? 
de matière et que l’autoabsorption des électrons 
doit y être faible. Comme dans les cas précédents 
il n’a pas été tenu compte de la rétrodiffusion de 
ces électrons sur la membrane supportant la source. 
La courbe théorique représentéesurla figure 4 a été 
obtenue en utilisant pour la couche X une proba- 
bilité d’autoionisation égale à 0,65/Z2, la correction 
d'écran proposée par Levinger ne semblant pas se 
justifier dans ce cas. La comparaison entre théorie 
et expérience est résumée dans le tableau suivant : 


22Na de 1,5 à 5 keV de 5 à 10 keV 
Spectre expérimental 38.107+ 23 10—4 
Théorie de l’autoionisation 12104 14102 


Ce résultat expérimental est en accord avec celui 
obtenu pour le sodium 22 par Hamilton, Langer et 
Smith [9]. Ces auteurs employant un spectromètre 
magnétique obtiennent au delà de 2 keV un spectre 
d'électrons négatifs décroissant moins vite et 
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d'intensité plus grande que ce que prévoit la théorie 
de l’autoionisation. 

Je remercie M. Jean Mérinis pour l’aide apportée 
dans la préparation des sources utilisées dans ces 
expériences. 


Discussion des résultats. Comparaison avec la 
théorie. — Si l’on considère les résultats expéri- 
mentaux relatifs à l’yttrium 90, au praséodyme 143 
et au sodium 22 et qu’on y ajoute ceux que j'avais 
précédemment publiés en collaboration avec 
G. Charpak [6] obtenus dans les mêmes conditions 
et relatifs au phosphore 32, au soufre 35 et au pro- 
metheum 147 on peut en tirer les conclusions sui- 
vantes : 
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1) Le spectre continu des électrons associés au 
rayonnement $ présente dans tous les cas le même 
aspect de décroissance rapide en fonction de l’éner- 
gle. Au début du spectre à 1 keV, l'intensité des 
différents spectres est pratiquement proportionnelle 
à Z, exception faite pour le praseodyme 143 pour 
lequel nous avons vu qu’une partie des électrons 
de faible énergie devait être absorbée dans la 
matière elle-même de la source. Au delà de 3 à 
4 keV, les Intensités des spectres quoique crois- 
sants avec Z, augmentent moins vite que propor- 
tionnellement à Z ; l’yttrium cependant fait excep- 
tion à cette règle ayant dans cette région une inten- 
sité plus faible. 

Les spectres continus des électrons associés au 


Spectres théoriques 
PS et 535 


Spectres expérimentaux 


rayonnement 8 des 6 radioéléments étudiés sont 
reportés sur la figure 5, ainsi que les spectres prévus 
par la théorie de l’autoionisation. 

Le désaccord important entre ces résultats expé- 
rimentaux et la théorie porte sur plusieurs points. 
D'abord les valeurs expérimentales sont beaucoup 
plus grandes que celles prévues par la théorie ; 
ensuite la théorie prévoit un spectre décroissant 
beaucoup moins vite dans le cas des éléments 
lourds ; enfin la théorie prévoit un spectre dont 
l'intensité décroît avec Z alors que l’expérience 
fournit le résultat inverse. Il faut aussi noter que si 
pour les éléments légers on extrapole à 0 Îes 
courbes expérimentales suivant une loi de crois- 
sance moins rapide que celle prévue par la théorie 
on peut obtenir une valeur totale de l'ionisation à 
peu près en accord avec la valeur théorique. Par 
contre pour les éléments lourds ceci n’est plus exact, 
puisque les intensités des spectres des électrons 
entre 1 et 13 keV seulement sont pour l’yttrium 90 
et le prometheum 147 respectivement 14.107% et 
20.103, les probabilités totales d’autoionisation 
étant respectivement d’après la théorie 6.107 et 
MALO, 


2) Les intensités des raies X, X ou Z que nous 
avons mesurées sont dans tous les cas en accord 
satisfaisant avec la théorie de l’autoionisation. 
(L'écart entre les valeurs théoriques et expéri- 
mentales est toujours inférieur à un facteur 2.) 

3) Les intensités des raies des électrons Auger 
que nous obtenons ne se prêtent pas à des con- 
clusions quantitatives précises mais sont en bon 
accord avec les intensités, connues des raies X et 
les valeurs admises des rendements de fluorescence. 

Il résulte de ces faits que le spectre continu des 
électrons ainsi obtenu ne provient pas en majorité 
des couches internes À et L comme le prévoit la 
théorie de l’autoionisation, mais qu’il faut au con- 
traire envisager une contribution de l’ensemble du 
cortège électronique. Dans la mesure où l’intensité 
de ce spectre est proportionnel à Z, ce qui est vérifié 
entre 1 et 4 keV environ. cela donne une probabilité 
d’ionisation par électron du cortège indépendante 
de Z. En résumé le manque total d’accord entre ces 
expériences et la théorie se manifeste par une dis- 
tribution spectrale des électrons émis simulta- 
nément avec le rayonnement 8 plus intense et 
d’une forme différente que celle prévue, et aussi 
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par une loi en fonction de Z inverse de celle atten- 
due. 

Parmi toutes les autres recherches expérimen- 
tales effectuées sur l’autoionisation en radio- 
activité B, la plupart ont porté sur la mesure de 
l'intensité des raies X des couches À ou Z. Ces 
mesures sont comme les nôtres (voir ci-dessus Pr 
et [6]) en bon accord avec la théorie de l’autoioni- 
sation et ont porté sur les éléments suivants : 
S2P:MO0) 5S 1117; Set 20 -M21et 13) APE 
et 15], Ra£ [12, 13 et 14]. Snell et Pleasonton ont 
abordé l’étude expérimentale de l’autoionisation en 
mesurant la charge du noyau de recul après désin- 
tégration par émission $ du %Kr ; la probabilité 
totale d’autoionisation dans l’ensemble des 
3 couches X, Let M a été trouvée égale à 9,9.1072 
et dans la couche NV égale à 10,9.107?. La théorie de 
Migdal fournit, en tenant compte de l’effet d’écran, 
la valeur 3,17.10 2, pour les 3 couches K, Let M 
et la théorie de Levinger une valeur plus faible ; par 
contre un calcul plus récent de Green [17] fournit 
pour les couches X, Let M la valeur 4,4.1072 et 
pour la couche A la valeur 9.102. Si cet accord 
peut être considéré comme satisfaisant, 1l faut 
cependant noter que la probabilité expérimentale 
de trouver une charge supérieure à 2 est nettement 
plus grande que celle prévue par la théorie. Pour 
l’étude directe du spectre des électrons émis simul- 
tanément avec le rayonnement Gil existe seulement 
l’expérience citée plus haut sur le ?2Na [9] et des 
expériences de Duquesne [18] sur les radioélé- 
ments %2P, %5$ et RaË£, qui par une méthode de 
coïncidences trouvait un spectre d’une intensité 
plus grande que celle prévue par la théorie. 

Les électrons négatifs signalés par Bruner [19] 
dans le rayonnement du #Sc n’ont pas été con- 
firmés par Blue et Bleuler [20]. 

Le phénomène d’autoionisation se présente aussi 
lors de la désintégration par capture électronique. 
Nous avons déjà publié [6] les résultats obtenus 
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dans les mêmes conditions expérimentales avec 
le 5Fe se désintégrant par capture Æ ; nous 
n’avons trouvé aucun électron associé dans ce cas 
à la raie X-X de 5,9 keV avec comme limites supé- 
rieures les valeurs 3 1074 dans la bande de 1 à 
7 keV et 0,5.1074 dans la bande de 10 à 14 keV 
et 0,5.10—4 dans la bande de 15 à 20 keV ; ces 
limites supérieures sont beaucoup plus grandes que 
les valeurs prévues par la théorie [22] (de 15 à 
1.106 par keV dans le domaine de 1 à 20 keV) mais 
beaucoup plus faibles que les intensités anomales 
que nous avons trouvées dans le cas de la radioacti- 
vité B dans le même domaine d’énergie. Miskel et 
Perlman [21] étudiant à l’aide d’un compteur pro- 
portionnel les phénomènes consécutifs au réarran- 
gement du cortège électronique ‘après désinté- 
gration par capture électronique de l’argon 37 ont 
obtenu un résultat en bon accord avec celui prévu 
par la théorie [22] (le spectre des électrons obtenus 
de 6 à 10 keV décroît moins vite que ne le prévoit 
la théorie, le désaccord n’excédant cependant 
jamais le facteur 2). Mesurant au moyen d’une 
chambre Wilson l’intensité relative des él2ctrons 
éjectés par autoionisation lors de la désintégration 
par capture À du fer 55, Porter et Hotz [23] ont 
obtenu dans le domaine d’énergie entre 30 et 
205 keV un résultat en bon accord avec la théorie. 
Cet accord a aussi été obtenu par Langevin [24] 
étudiant au compteur proportionnel les électrons 
Auger et le fond continu d’électrons émis lors de la 
désintégration par capture électronique du ger- 
manium 71. Par contre Soltysik (Thèse non publiée, 
voir référence [9]) étudiant par spectrométrie ma- 
gnétique les électrons éjectés dans le domaine de 10 
à 400 keV lors de la désintégration par capture X 
du berylium 7 a obtenu un résultat de 10 à 40 fois 
plus grand que celui prévu par la théorie et un 
spectre décroissant moins vite en fonction de l’éner- 
gie que celui théoriquement prévu. 
Manuscrit reçu le 20 janvier 1 960. 
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OSCILLATIONS COLLECTIVES D'UN SYSTÈME DE TROUS EN PRÉSENCE DE PHONONS 


Par A. COUMES, 


École d'Ingénieurs Électroniciens, Grenoble. 


Résumé. 


L’effet d'écran électrostatique dû à des particules légères entourant des particules 


lourdes de charge opposée peut, dans certaines conditions, être caractérisé indépendamment de 
la nature des particules écrantées. Les oscillations collectives des trous sont traitées par analogie 
avec les phonons ; leur amortissement semble très faible. 


Abstract. — The screening effect of light particles surrounding heavy particles (the charge of 
which is of the opposite sign) can be characterized independently of the screened particles. The 


collective oscillations of holes are derived by analogy with phonons. 


oscillations is shown to be very small. 


L'étude, suivant les méthodes de traitement col- 
lectif, du comportement des milieux contenant des 
particules électrisées positivement et négativement 
conduit à chercher à généraliser ces méthodes au 
cas où le milieu étudié contient deux systèmes de 
particules positives lourdes (ions et trous) et un 
système de particules négatives légères (électrons). 
En effet, un système de particules lourdes entou- 
rées d’un écran de particules légères et de charge 
électrique opposée peut être Le siège d’oscillations 
collectives : par exemple les oscillations collectives 
des ions écrantés par les électrons constituent les 
phonons [1],[2], [3]. De même Nozières et Pines [4] 
ont montré la possibilité de calculer les fréquences 
collectives d’un système de trous écrantés par les 
électrons, ces oscillations constituent les demons 
(Dependant Electronic Motion). Il semble donc 
intéressant d’étudier les causes d’amortissement 
de ce système d’oscillations. Nozières et Pines ont 
montré que l’amortissement d’un demon par inter- 
action avec les électrons, avec les trous ou avec 
d’autres demons était parfaitement négligeable ; 
mais l’on peut considérer également une autre 
cause d’amortissement : les oscillations du réseau 
cristallin dans lequel évolue le système de 
demons [5]. Dans cet articlenous présenterons uu 
procédé d’étude du système de demons en présence 
du système de phonons. 

Nous considérons un système constitué par des 
électrons (A), des trous (B) et des ions (C) ; donc 
deux sortes de particules lourdes : ions et trous 
dont les actions électrostatiques sont écrantées par 
les électrons. 

Exposons brièvement les principes qui nous con- 
duisent à l’écriture d’un Hamiltonien où apparais- 
sent les opérateurs correspondant aux deux types 
d’oscillations collectives dont nous voulons étudier 
les interactions. 


1. Usage des fluetuations de densité. — Le mou- 
vement de chaque type de particules est décrit en 
faisant usage des fluctuations de densité. 


FRE 


The damping of these 


Par exemple pour les particules A : 


pA(x) — 2 S(X — Lx) — 2 eik.(x—x y) 
d’où é 


A : 
CR — > e—k.x\ 


quireprésente Ja kiène composante de Fourier de LA 
d’où l’écriture de | Hamiltonien du système : 


2 
PE 


PÉ 
Je » LACS 
B 2mp 


A 27MA C 20 


a 97 B 2 
+ 2 7 (ea ee +epek + opt). (1) 
£Z0 k° 
La précision # =£ 0 signifiant qu’une particule A 
par exemple est en mouvement dans un milieu dont 
le potentiel est déterminé par les autres particules. 


2. Introduction de l’effet d’écran et limitation 
aux interactions à grande distance. — L'effet 
d’écran électrostatique que les charges négatives A 
exercent sur les particules B et C sera caractérisé 
par l’introduction de grandeurs Ae et Ae’, On rem- 
place les charges e$ par des charges eg — Ae répar- 
ties avec la même densité 65; de même aux 
charges eç nous substituons des charges e — Ae' 
réparties avec la densité &°, D’où l’Hamiltonien : 


Âe') 2 


LutE A B 
2 (%) Lea se + (eg — Ac) p& + (ec 


x 4T A *B C -*B AN C 
= à 7e [ex Aepk pe + ecAepi Pr + ea Ae’ px px 
EE, ë 


B «+0 Âe B *B 
+ ep À px pk + (a 5) Ae pk Ph 


Âe’ * . 
de (eo = ] Âe’ ck ei | + HKr. (2) 


) 


8 et p*£ désignant les imaginaires conjugués 
de PF et p£. 
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En remarquant que seules les actions à grande 
distance entrent en ligne de compte pour régir le 
comportement collectif de notre système, nous 
pourrons limiter les sommations sur k# à k# < kG 
(ke représentant l’inverse de la longueur d’onde de 
Debye), quitte à ajouter un terme Y = Hp, ren- 


KG 


dant compte des actions à courte distance. 


Remarques : 


1. Les grandeurs Ae et Ae' caractéristiques de 
l'effet d’écran ne sont pas arbitraires : dans la som- 
mation sur # ce sont des fonctions de k dont nous 
verrons plus loin l’évaluation. 

2. Dans nos sommations ce sont les groupes de 
lettres ÆA, ÆB, AC qui constituent en réalité les 
indices de sommation. Lorsque la nécessité appa- 
raitra nous désignerons par k, L ete. les rangs de 
différentes composantes de Fourier. 


3. Introduction de variables collectives. — Dans 
l’'Hamiltonien précédent (2) nous pouvons séparer 
les termes correspondants à plusieurs sortes d’oscil- 
lations collectives en utilisant la méthode indiquée 
par Bohm et Pines [6]. 


a) PLASMA GLOBAL. — Dans le terme quadra- 
tique de (2) introduisons la variable collective P; 
en même temps que nous introduisons un groupe 
de conditions subsidiaires assurant la conservation 
des états propres : 


Pr d—=0 pour tout k k < ko. 


En appelant Q, la coordonnée collective associée 
à la quantité de mouvement collective P, et en 
effectuant la transformation canonique : 


d — eiS/% 


avec pour fonction génératrice : 


1/2 
en (5) 
k%g \k 


RMC RAC CCI ex 


L’Hamiltonien (2) prend la forme : 


CL CNARES 


NPA SP | 
H=D—=+X ES: 2 Si PR 


A 2MmA B 2MB C2 k<k = 
F TVAN CA) (ep — Àe) n8 (ea — Âe’) nc " 
ir) [22 LE Qe Où 
MA mp Inc 
&T\ 1/2 à ue Aile 
+ à al Ok [E CAOPRT RES (Ps a =) 
kckg \ K A TA 2 


S ; uk ak 
+ Zen e-ikxg — (P —) 
B mp 2 


x \ ile 1 
EL ec e—K.XG — (Pr — —)| 
C 


mc 
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4T * C »B RATE 
Je > 7e [ex Ae oi Ces + eu Aer px + ea Ac p% px 


* Ae B *B 
+ €B Ae’ ch a dE (en — ue Pk Ce 


Ae’ } 40 
+ (eo —") Aë çk a] 
POUCES. (3) 


La deuxième ligne de cette expression corres- 
pond à un système d’oscillateurs indépendants 
dont les fréquences sont données par : 


- 


À b AN 
ot GT À , nplen — Ae)?, nolec — Ace | 
N PATES mp É mo 


ha, A», Aç représentent les nombres de parti- 
cules À, B, C par unité de volume. 

Ce sont les oscillations d’ensemble du « système 
global » de trois espèces de particules. Ces oscil- 
lations de fréquence très élevée sont couplées aux 
particules À, B, C par les termes qui constituent les 
troisième, quatrième et cinquième lignes del Hamil- 
tonien (3). Ce couplage est extrêmement faible [6] ; 
nous ne considérerons plus les termes P, Q; dans 
notre étude car ils n’ont que peu d'influence dans 
le domaine qui nous occupe. 

U représente le terme de « phase aléatoire » de la 
forme : 

FUUR 
24 Te a 9e @ 
nous le négligeons par la suite [6]. 

Finalement nous restons avec l’'Hamiltonien sui- 

vant : 


4T A *B C *B / + 
FL Ta [es Acpk pe —+eoAee®px + ea Ac’ pR pr 
k<RG 


/ B * Ae * Âe’ ñ %* 
+ eg Ae’ ph pR° + (en —7) Ae pk pe + (+) Ac! | 


(4) 
et l’ensemble de conditions subsidiaires, qui par la 
transformation canonique est devenu : 


A B } E 
lea 6% + (es —Ac) px + (ec - - Ae') ok] d = 0 pour k < ke. 


b) PHONONS ET DEMONS. — Séparons en deux 
parties chaque composante LE et & des fluctuations 


de densité : un terme individuel et un terme collec- 
FIRE D 


B ; 
CRE AOIE 
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Les nouvelles variables collectives ainsi intro- 
duites sont reliées aux fluctuations de densité et 
aux coordonnées collectives X, et g, par la trans- 
formation canonique d = ei avec 


: [ar ep Âe RUE 
LEEDS B es AL TecAe 
E<k,, k2 Pr Â% Ke °x x. 


La partie intéressante de l’'Hamiltonien devient 
alors : 


P3 Des 2 
TRE >» RTE Æ B de 4 C 
A 2A B 2m © 2mc 


1 + ET np ep Àe , 
Du [re Re 2 à, x] 
ke 2 mp 


4r no ec Âe’ 


1 + ; 
<ke: 2 [re PRE qe ai 


&T e* Ae 172 A+ Ae ec 1/2 * 
+, 2 [() gi ri + (30) TR Pk 


— - 
kg 2 eB 


ame Ae\1? À . fAc'em\"? + 
TE Etre + Fe PE + ane PR TR 


US) 


O 


ee LT np €B Âe 


4Tr ne ec Âe’ 
(> et D  — 


mp mo 


Et les conditions subsidiaires résultant de l’in- 
troduction des variables collectives p} x de- 
viennent : 


pr b=0—/f ie d — 0 


V= eg Àe’ 
k2 


TX 


kr ep Ae 
k2? 


Les oscillateurs cemons et phonons apparaissent 
dans les deuxième et troisième lignes de (5). Pour 
les identifier nous cherchons à calculer les fré- 
quences correspondantes Q, et w,. Pour calculer la 
fréquence des demons Q, Nozières et Pines ont 
suggéré une transformation canonique qui, élimi- 
nant au premier ordre les interactions demons- 
particules A dans l’ Hamiltonien, fasse également 
disparaitre les variables demons des conditions 
subsidiaires. On obtient ainsi une relation entre Ae 
et k ; et cette relation exprime que l’introduction 
de l’effet d’écran Ae et celle des variables collec- 
tives constituent deux opérations compatibles cor- 
respondant aux demons. La même technique appli- 
quée aux phonons (Pz, 4, x) conduit à des valeurs 
correctes pour la vitesse du son V, = «z/k. 

Nous proposerons ici une méthode plus rapide 
pour obtenir la valeur de l’effet d’écran relatif aux 
ions et pour étendre cette valeur aux trous. 

La vitesse de propagation des ondes longitu- 


mx V=0—/f FR y = 0. 
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dinales dans un cristal, calculée à partir des inter- 
actions 1ons-électrons est donnée par : 


OÙ Mn 22 
V2 2 E ebno mA 


9 LA, eà NA MC 


FL] La ” . D r 
E représentant l’énergie de Fermi des électrons (*). 
En comparant cette expression avec 
hr no ec Ae’ 
12 v2 ds oÈ __ AT Ang ec ÀÂé 
; mc 
Il vient 


12 (6) 


«A représentant la pulsation des oscillations de 
plasma des électrons : 
2 AT na eÀ 
OA 
LUN 
L'expression (6) montre que l’effet d’écran relatif 
à la Aïème composante des oscillations d'ensemble 
des ions s’exprime uniquement à partir de termes 
relatifs aux particules À constituant l’écran. Autre- 
ment dit : l’effet d’écran caractérisé par le rap- 
port AeJe est indépendant de la nature des parti- 
cules écrantées. Pour les trous nous pouvons 
prendre également : 


ne ie 


—* ee 
€B 3 © mA 


Cette expression est en accord avec la formule 
donnée par Nozières et Pines. Les vitesses de pro- 
pagation des demons sont de 100 à 1 000 fois plus 
élevées que celles des phonons. 

DIE ue Eè 
Vh = 08%. 
3 MB na eÀ 

c) REMARQUE. — Dans le traitement précédent 
nous avons introduit les variables collectives rela- 
tives aux trous écrantés ainsi que les variables 
collectives relatives aux ions écrantés, mais nous 
n’avons pas considéré de variables analogues pour 
les électrons. La raison en est que les interactions 
électron-électron ne modifient pratiquement pas 
l'interaction entre électrons et phonons de grande 
longueur d’onde ainsi que Pines l’a montré [3]. 
L’approximation de Hartree est donc valable pour 
le calcul des éléments de matrice correspondant 
aux interactions électron-phonon ou électron- 
demon ; nous considérerons les électrons comme 
indépendants, quitte à faire intervenir leur masse 
effective. 


{*) Compte tenu de la neutralité du milieu, qui pour deux 


espèces de particules s’écrit noleg| — naleal et en intro- 
duisant le numéro atomique des ions Z — legfea| la rela- 
tion devient : 
Pere 
z 3 °mo 


Ici la neutralité du milieu s'exprime par 


na ea + ngeg + ncec = 0. 
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4. Étude de l’amortissement des demons par les 
phonons. — La partie intéressante de l'Hamil- 
tonien (5) se borne à : 


D2 6} * - 7* 
Ti FES HS ÿ Le Th + OÙ Xr XK] 
C 


AMATIDA k<e, 


dl -_* * 4x ea Âe A _* 
+ D = [px pe + ox qu gx] + > (52) 0% Th 
E<k 2 L<ke B 


4 Âe’ * 
TEA <) À Fe (7) 
k<kG 


k2 ea 


a) INTERACTION DIRECTE DEMONS-PHONONS.. — 
La collision demon-phonon ne pouvant conserver 
à la fois la quantité de mouvement et l’énergie 1l 
est nécessaire de considérer des états virtuels inter- 
médiaires et l’on calcule la probabilité d’amortis- 
sement par unité d’angle solide autour de la direc- 
tion X par la formule classique : 


U1,i Ui,2 
à E; +: EF; 


OT 


D — 
ñ 


D, 


U;s, Ui2z : éléments de matrice pour la tran- 
sition de l’état initial à l’état intermédiaire et de 
l’état intermédiaire à l’état final. 

E, — E; différence d’énergie entre état initial 
et état final. 

D, densité de l’état final. 

Finalement on trouve une probabilité d’amor- 
tissement de l’ordre de : 

Hi OS CT 
avec À < ko. 

Cet amortissement est minime tout au moins 
aux basses températures. 


b) AMORTISSEMENT AVEC INTERVENTION DES 
ÉLECTRONS. — Nous avons cherché si la partici- 
pation des électrons pouvait conduire à un amortis- 
sement plus important des demons : en effet dans 
le cristal nous disposons d’un « grand nombre » 
d'électrons pouvant participer aux échanges d’éner- 
gie et de quantité de mouvement. Les processus 
envisagés sont du type suivant : 


Électron À + demon — Électron A’ (virtuel) 
— Électron A” + phonon. 
Le schéma du calcul est le suivant : soit un 
demon k, nous cherchons sa probabilité de collision 
avec un électron, le résultant étant l’émission d’un 


phonon /, puis l’on somme cette probabilité sur 
tous les états finaux possibles (valeurs possibles de). 
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Le calcul est assez simple tant quel’on ne considère 
que l’émission de phonons / dont les directions de 
propagation sont voisines de celle du demon initial : 
la seule anomalie pourrait se manifester pour ! —#, 


Ie ile 


mais une discussion montre qu’alors nous n’aurions 
pas d'électrons correspondant à un tel processus. 
Au total nous trouvons une probabilité d’amor- 
tissement minime. 3 

Si maintenant nous faisons intervenir toutes les 
directions possibles pour le phonon émis /, la condi- 
tion /  k ne suffit plus pour éliminer les cas dou- 
teux, cependant une discussion nous à convaincu 
que l’amortissement des demons demeurait négli- 
geable. 


5. Conelusion : recherche des oscillations collec- 
tives des trous. — Les calculs précédents impli- 
quent que le comportement électrique du milieu 
étudié soit essentiellement déterminé par les parti- 
cules légères A (électrons). Ce qui se traduit par les 
inégalités suivantes. 


Gi D &B oÀ > wË 
Va> Vs Va Ve 
kV 10 O0 EVA Oo: 


& : fréquences de plasma. 

V : vitesses moyennes. 

Nous devons donc chercher les demons dans des 
corps possédant des rapports de masse effective 
électronique faible, en regard des rapports de masse 
effective des trous. Les semi-métaux du genre 
bismuth semblent convenir puisqu'on rencontre 
des rapports mÂ/m, allant de 0,5 à 0,005 tandis 
que les rapports m$/m, varient de 1 à 5. Le sys- 
tème d’oscillations collectives des trous excités par 
un champ magnétique semble pouvoir se manifester 
par l’absorption d'énergie électromagnétique. 


Manuscrit reçu le 21 novembre 1959. 
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SUR LE CALCUL DU COEFFICIENT DIAMAGNÉTIQUE DES ATOMES LIBRES 


par CLÉMENT COURT Y, 
Chaire de Chimie physique de la Faculté des Sciences de Lyon. 


Résumé. — À partir des énergies de première ionisation des atomes on détermine un effet d’écran 
que l’on corrige ensuite de «sauts d'écran » expérimentaux. Cet effet d'écran corrigé permet, sans 
modification des nombres quantiques principaux, de calculer les coefficients d’aimantation dia- 


magnétiques des atomes libres. 


La méthode laisse entrevoir le même calcul pour les atomes combinés et les molécules. 


Abstract. — Starting from first ionization energies of atoms, a screening effect is determined, 


which is then corrected for the experimental ‘screen jumps ”. 


enables, without modification of the 
tibility of free atoms to be computed. 


cules. 


On sait, en se basant sur la théorie de Lange- 
vin [1], que le coefficient d’aimantation diama- 
gnétique, x, d’une orbite électronique rapporté à 
Patome-gramme, donc multiplié par le nombre 
d’Avogadro, N’, peut s’écrire : 


D ÉR S r (1) 


(e et m sont la charge et la masse de l’électron, 
r le rayon de l’orbite). 

En unités e. s. cgs, x vaut — 2,462.10731 r? et 
en unités e. m. cgs, plus généralement utilisées 
dans ce cas, y — — 2,862.1010 r2. 

Or, si l’on prend r — 0,529 À pour l’atome 
d'hydrogène on trouve pour x en unités e. m. 
x = — 2,862.2,798.1077 = — 8,009.1077 alors 
que la valeur expérimentale est de l’ordre pour 
l’atome combiné de 20.107. 

Lorsqu’on aborde des atomes plus compliqués et, 
en particulier, ceux des gaz inertes, le désaccord 
est complet. 

Van Vleck [2], d’une part, et Pauling [3] de 
l’autre, ont calculé par la mécanique ondulatoire le 
rayon r des orbites qui devaient intervenir dans le 
calcul des coefficients diamagnétiques et posé : 


Re nr. = 
à = nb" 2 


Z est le numéro atomique de l’atome auquel appar- 
tient l’orbite, s est une constante dite « d’écran », 
n et Lles nombres quantiques connus, à, le rayon 
de l’atome d'hydrogène, non excité. 

Précisément dans le cas de l’atome d’hydrogène, 
nd 7 l}s —0\i— Oetr” — a; (ff) 
(5/2 + 1/2) = 3aë, et le calcul précédent recon- 
sidéré donne y — — 24,027.107 acceptable. 
Lorsque, maintenant, la constante d'écran, s, n’est 
pas nulle, il faut connaître sa valeur et le problème 
devient, de ce fait, difficile. 


(2) 


This corrected screen effect 


principal quantum numbers, the diamagnetic mass suscep- 


The method suggests the possibility of the same calculation for combined atoms and for mole- 


Le problème de la constante d’éeran. — Nous ne 
citerons, des travaux concernant cette constante, 
que le plus connu : celui de Slater [4]. 

Slater a établi, en fonction du nombre d’élec- 
trons extra-nucléaires, une suite de valeurs d’effet 
d’écran très simple et commode, choisie pour que 
ces valeurs combinées avec un nombre quantique 
effectif, n', différent de n, on trouve des nombres 
corrects pour les coefficients y. 

Ces nombres quantiques effectifs, n', obéissent 
aux rapports : 


DD AM EN EN ER MEN) EME ES 
et finalement on aurait pour chaque orbite : 
Le (w' ie ) (n° + 1) 


07 É 3 
106 y — 0,79326 VAN E (3) 


avec les valeurs actuelles des constantes physiques 
figurant dans l’expression de 4. Comme (nous le 
verrons plus loin) le calcul avec le rayon r ci- 
dessus conduit à : 


— 106 y = 0,79326 


n CU USE Te k 
fr | 


cn étant un autre effet d’écran, et que les résultats 
obtenus sont indiscutablement bons pour les 
atomes libres, il va de soi qu’une différence très 
nette doit se faire jour entre s et on. 

Dans les cas où n/ — n, c’est-à-dire pour n < #4, 
s et on ne seraient égaux que pour  — n— 1 (°) 
donc pour les électrons s quand n = 1, p quand 
n = 9, d quand n — 3, mais il n’y a pas d’élec- 
trons d dans les atomes de la troisième période. 

Or, même pour ces valeurs de /, l’effet d’écran de 


(1)7 C'était dans les vues de Sommerfeld l’orbite à excen- 
tricité unité | 
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Slater diffère beaucoup du 6, dont il sera question 
plus loin. 

L'idée de l’effet d'écran de Slater découle de 
l'interaction supposée des électrons en admettant, 
toutefois, que les plus externes n’ont aucun effet 
sur ceux situés en dessous : seuls produiraient un 
effet d’écran sur un électron atomique donné les 
électrons sous-jacents d’une part et ceux de même 
niveau énergétique d’autre part, c’est-à-dire les 
électrons dits équivalents (n et l'identiques respec- 
tivement). 

Il y a alors, pour lui, trois valeurs pour les effets 
d'écran : la valeur 1 qui caractérise un électron 
voisin du noyau pour un électron qui en est très 
éloigné, l’effet 0,85 lorsque l’écart interélectronique 
est moindre et 0,35 quand les électrons sont très 
voisins (même case quantique par exemple). 

La valeur 1, à retrancher de Z, si l’on imagine un 
électron comme faisant réellement écran entre le 
noyau et un autre électron éloigné, pourrait s’ad- 
mettre, à la rigueur, comme un maximum excep- 
tionnellement atteint. Il n’apparaît que lorsque les 
électrons les plus externes ont atteint, puis dépassé 
le stade à nombre quantique principal 3 et c’est 
sans doute une des raisons pour lesquelles, s étant 
trop fort, au-dessus de la valeur » — 3, Slater a 
diminué la valeur de ces nombres quantiques pour 
utiliser les nombres quantiques effectifs plus faibles. 

Mais il y a une autre raison : le fait que le diama:- 
gnétisme n’est pas un phénomène quantifié, en sa 
base. Les nombres quantiques apparaissent dans 
les calculs de 7 en raison du fait que les rayons 
orbitaux, qui eux sont quantifiés, interviennent. 

On peut imaginer un autre système d’eftets 
d’écran. 

Considérons l’énergie de première ionisation, Wi, 
d’un atome de numéro atomique Z. En mécanique 
ondulatoire non relativiste elle peut s’écrire, si 
on néglige l’interaction des électrons : 

’ 72 TE nes 
Wi = (5) 
h étant la constante de Planck. 

Si l’on veut, comme Slater, tenir compte de cette 
interaction, on peut introduire l’effet d’écran de 
sorte que, remplaçant Z? par (Z — 6)? on ait 
toujours l’énergie associée à l’électron déplaçable. 
Mais alors, dès que o est introduit, on peut ne plus 
admettre que cette énergie représente une énergie 
d’ionisation, car le départ d’un électron peut modi- 
fier les effets d’écran dans l’ion résiduel et intro- 
duire une réorganisation des énergies associées aux 
électrons de l'ion. 

_ Prenons, par exemple, les énergies de première 
ionisation de trois éléments qui se suivent : Al, Si, 
P. Elles valent 5,985 ; 8,148 ; 10,977. Les nombres 
atomiques respectifs sont 13 ; 14 ; 15. Si l’on divise 
ces énergies par Z? on trouve les valeurs qui 
devraient être une constante et qui ne le sont pas : 


No & 


0,03541 ; 0,04157 ; 0,04878. Si l’on introduit leffet 
d'écran de Slater il est incontestable qu'on se 
rapproche d’une constante. Ces effets d'écran 
seraient 9,50 : 9,85 ; 10,20. En divisant cette fois 
par (Z—s)? on trouve 0,4885 puis 0,4731 et 
enfin 0,4764. Mais si l’on fait l’opération inverse, 


c’est-à-dire si l’on suppose o inconnu et qu’on le 


calcule on trouve respectivement : 11,0124 ; 11,680 ; 
12,316, qui, naturellement conviennent mieux et 
sont des valeurs expérimentales, en admettant, b'en 
entendu, que l’expression : 


W;i = (Z— o)227n° met [n° h° (6) 


convienne bien à l’énergie de première ionisation 
des atomes. 
La valeur de o qui en résulterait s’écrirait : 


CO —= 77 D nk VW; (7) 
en posant : 3 
k = Lee Von. (8) 


Remarquons que si W; était une constante, la fonc- 
tion o — f(Z) serait représentée par des segments 
de droite parallèles, s’écartant de la droite 6 = Z, 
de quantités croissantes à chaque accroissement 
d’une unité du nombre quantique n. 

Nous avons calculé les constantes d’écran, ©, à 
appliquer à l’énergie associée à l’électron dépla- 
çable de chaque atome et qu’il faut fournir pour 
l’arracher, ce qui revient à admettre un modèle 
d’atome à électrons indépendants. 

Le tableau I donne ces constantes, o, pour la 
plupart des atomes, en utilisant les énergies d’ioni- 
sation données dans l’ouvrage de G. Herzberg [5] 
au tableau 18 des pages 200-201. Nous y avons fait 
figurer les symboles des électrons déplacés en nous 
basant sur les « états normaux » des atomes et des 
ions, représentés, comme chacun sait, par ses 
nombres quantiques atomiques résultants et la 
multiplicité spectrale, u, sous la forme “L;. Pour 
les états p des électrons nous avons marqué la sub- 
division en p, et p, suivant l’ordre de remplissage 
des cases quantiques. Quant aux électrons s recou- 
vrant des sous-couches d ou f, nous les avons dési- 
gnés par sd et sf lorsque ces sous-couches n'étaient 
pas complètes et par sD et sF lorsqu'elles l’étaient. 

Enfin, certains symboles sont accentués s’d (Cr 
et Mn) sd’ (Nb, Mo, Te, Ru) lorsqu'une irrégularité 
se manifeste dans la distribution électronique, ce 
qui modifie l’effet d’écran de manière vraisem- 
blablement irrégulière. 

. Nous avons poussé les déterminations de 6 
jusqu’à la décimale nécessaire pour que, faisant le 
calcul au rebours on puisse retrouver la valeur 
expérimentale de (Z — 6)?. Il va de soi que l’effet 
d'écran dépendant surtout du nombre d’électrons 
extra-nucléaires, o est une fonction de Z comme 
nous l’a montré l’équation (7). Or, si l’on construit 
la courbe 6 — f(Z), on obtient bien le graphique 


N°4 


COEFFICIENT DIAMAGNÉTIQUE DES ATOMES 
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TABLEAU I 
Z ATOME  ÉLECTRON oc Z ATOME  ÉLECTRON o | Z ATOME  ÉLECTRON 6 
2 He ds 0,6558 31 Ga ap, 28,3476 60 Nd Gsdf 55,912 
3 Li 2 1,783 32 Ge Ap; 28,9124 | 61 Pm 6sdf ? 
le Be 2s 2,345 33 As Api 24,482 62 Sim Gsdf 97,817 
ÿ B 2p1 3,440 34 Se hPo 30,169 63 Eu Gsdf 59132 
6 (6 2p: 41798 | 35 Br 4pa 31,271 64  Gd 6sdf 59,785 
7 N 2p, 4,9346 | 36 Kr 4po 31,949 65 EETD Gsd} 60,7851 
ce 2pa 6,002 37 … Rb 5s 34,2397 | 66  Dy 6sdf 51,754 
9 F 2p» 6,7398 | 38 Sr 5s 34707 67 Ho Gsdj ? 
10 Ne 2p 7,485 39 Y 5sd 35,505 | 68 Er Gsdf ? 
11 Na 38 9,1578 40 Zr 5sd 36,4304 69 Tm Gsdf ? 
12 Mg 3s 9,754 41 Nb 5sd’ ? 70 Yb Gsdf 65,662 
13 Al an 11,0124 | 42 Mo 5sd’ 38,3212 | 71 Lu Gsdf ? 
14 Si 3p, 11,680 43 Te 5sd’ ? 72 Hf 6sdF ; 
15 12 3P1 12,316 44 Ru sd’ 40,2346 73 Ta 6sdF ? 
16 S 3Po 13,387 45 Rh 5sd' 11 23200) TV 6sdF 69,366 
17 CI SP» 14,0699 | 46 Pd kd 42,8726 | 75 Re 6sdF ? 
18 A 3P2 14,776 47 Ag 5sD 43,2735 | 76 Os 6sdF 71,197 
19 K 4s 16,7443 48 Cd osD 43,9382 77 Ir 5d ? 
20 Ca 4s 17,3211 | 49 In Sp, :5,7833 | 78 Pt 6sdF 73,148 
21 Sc 4sd 18,19 50 Sn Sp, 6:3337 | 39 * Au 6sDF 74,066 
7202 HE ksd 19,167 51 Sb 5P: 4k7,0438 80 Hg 6sDF 74,751 
23 V 4sd 20,188 59 Te 5Po 17,9382 | 81 TI GP 76,988 
24 Cr &s'd 211817 | 53 I 5Po 48,582 82 Pb 6p: 77,576 
25 Mn &s'd 22 045 54 Xe Spa L9 67e UP 6p 78,395 
26 Fe 4sd 22,9624 55 Cs 6s 0157967 84 Po GP» “4 
97 Co sd 23,9663 | 56 Ba 6s 52,29 85 At 6P) ? 
28 Ni ksd 25,005 56 La 6sd 593,177 86 Rn 6P» 80,674 
29 Cu 4sD 25,992 57 Ce 6sd/ 53,836 87 Fr 7s ? 
30 Zn 4sD 26,679 59 Pr 6sdf 55,078 88 Ra 7s 83,647 


prévu à quelques anomalies près particulièrement 
intéressantes (#ig.). Il est constitué par des portions 
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de droites ayant toutes grossièrement le même 
coefficient angulaire, mais des ordonnées à l’origine 
régulièrement croissantes d’une portion à l’autre. 


Si l’on considère ces accroissements successifs des 
ordonnées à l’origine, que nous appellerons « sauts 
d’écran » on constate qu’ils ont lieu quand on passe 
des électrons enlevés, 1s à 2s, 2s à 2p,, 2p; à 2Pa, 
2p;à 35,35 à 8D1) 2D1 à DPar DD» à 48... ÉbC.. C'est-dr 
dire quand le nombre quantique n s’accroit d’une 
unité [prévu par l’expression (7)] ce qui se traduit 
aussi par le passage d’une fonction d’onde à symé- 
trie non sphérique (électron p) à une fonction 
d'onde à symétrie sphérique (électron s). Il y à 
exception pour les sauts d’écran dus au passage 
d'un électron p, à p,, done quand les cases quan- 
tiques p se remplissent. 

Dans le système de Slater les sauts d’écran 
n’apparaîtraient que dans les passages p - s mais 
n’existeraient pas aux passages S-p3 OÙ Pr - Po. 

Il nous a semblé que ces anomalies étaient du 
plus haut intérêt et ces sauts d'écran vont inter- 
venir indirectement dans le calcul des coefficients 
diamagnétiques des orbites. 

Du tableau [ qui donne les valeurs de c on peut 
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TABLEAU II 
ÊLECTRON ParT SAUT ÉLECTRON Panr SAUT ÉLECTRON ParT SAUT 
15 0,3279 4p: 0,5672 Gsdf 0,8962 | 
0,4852 0,472 0,3458 
25 0,602 BP» 0,665 Gsdf 0,8962 
0,3477 | 1,7494 __ 0,0622 
2p, 0,7473 5s 0,5343 Gsdf 0,8962 
0,53215 Fe — 0,1164 
2p, 0,67275 5sd 0,8899 Gsdf 0,8962 
1,0766 0452 0,1126 
3s 0,5962 5sd’ 0,9378 Gsdf 0,8962 
0,6066 0,4188 
3p: 0,6518 Ad 1,6380 Gsdf 0,8962 
0,3765 — 0,2638 ; — 0,2432 
3Do 0,6945 5s1) 0,6647 Gsdf 0,8962 
1,3915 1,15485 0,1039 
As 0,5768 5pP, 0,65025 Gsdf 0,8962 
—— 0,1301 0,22165 0,0727 
Asd 0,9990 5D» 0,67275 6sdf 0,8962 
0,1304 2,0197 
4s’d 0,8633 6s 0,4933 Gsdf 0,8962 
= 0,1039 
Asd 1,0213 Gsd 0,8570 Gdsf 0,8962 
0,300 02072 
AsD 0,687 Gsdf 0,8962 
1,014 0,3232 
Gsdf 0,8962 
6sdf 0,8962 
6sdF 0,9414 
0,233 
6GSDF 0,685 
1,953365 
GP, 0,70345 
0,1686 
6P» 0,70347 
2,063 
7s 0,455 


extraire les parts contributives à l’effet d’écran, 
apportées par chaque électron, étant bien entendu, 
comme l’entend aussi Slater, que l’ensemble des 
effets d’écran préexistants demeure, et on cons- 
tate alors, comme on le voit au tableau IT, que 
chaque électron de même niveau a un apport iden- 
tique et que cet apport est presque toujours infé- 
rieur à l'unité. Dans quelques cas particuliers, élec- 
trons 4sd et aussi dans celui du palladium (4d), la 
part contributive dépasse l’unité. Ou bien ces 
valeurs sont fausses ou bien, alors, il faut admettre 
que l’apparition de ces électrons provoque un bou- 
leversement de l’ensemble électronique sous- 


jacent, hypothèse qui n’est peut-être pas à exclure. 

Il va de soi que ces apports d effet d’écran dé- 
coulent des valeurs numériques qui nous ont servi à 
décrire les portions de droites. Parfois l’effet expéri- 
mental s’écarte un peu de la valeur que donne- 
raient ces apports. Le tableau III montre que pour 
15 éléments sur 88 considérés il y à un très léger 
écart que nous désignerons par Gcue. — Gexp. dont 
on tiendra compte. 


L'adaptation des effets d’écran aux caleuls de 
diamagnétisme. — Si l’on calcule les coefficients 
d’aimantation avec les effets d'écran o ci-dessus on 
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TABLEAU III 


COEFFICIENT DIAMAGNÉTIQUE DES ATOMES LIBRES 


o 
ÉCART ENTRE QUELQUES VALEURS DE L'EFFET D'ÉCRAN CALCULÉES ET LES VALEURS EXPÉRIMENTALES 


ÉLÉMENT ÉLECTRON Cale. Goaner. (Scale. — Oexper.) 
G 2p, 4,1873 4,1798 0,0075 
F 2p, 6,67475.  6,7398 —_ 0,0065 
Si 3p. 11,6642 11,6800 —_ 0,0158 
CI 3Pe 14,0815 14,0699 0,0116 
Ge &p, 28,9148  28.9124 0,0024 
Br 4 po 31,284 31,271 0,013 
Sn Sp, 16,39355  46,3337 0,05985 
I 5Po 18,61095 48,582 0,02895 
Pb 6p, 77,69155 77,576 0,14555 
Ti &sd 19,189 19,167 0,022 
Co &sd 23,9837 23,9663 0,0174 
Mo 5sd’ 38,422 38,3249 0,1000 
Ru 5sd’ 40,2968 40,2346 0,0622 
Os 6sdF 71,2652 71,197 0,0682 
W 6sdF 69,3824 69,366 0,0164 


trouve des valeurs manifestement trop fortes parce 
que les effets d’écran sont trop forts. 

Remarquons que le calcul qui donne © basé sur la 
mécanique ondulatoire non relativiste admet impli- 
citement un rayon orbital probable : 


n? h? n 


7 —>5 = Mo — 
LT? m(Z — 6) e? VASE) 


(9) 


& étant le rayon de l’atome d'hydrogène non excité. 
On peut recalculer o que nous nommerons o’ en 
remplaçant r’ par r de Van Vleck et Pauling. 
L'énergie dite d’ionisation W; devient W; : 


(Z — 0’)?27°? me 


BU +1) — 1 


(10) 


TV; — 


PA 
pa 


15 
n D2 ne) 
nh V 9 ñn 


Reprenant les énergies de première ionisation 
expérimentales égalées à W; on peut déterminer 6”, 
construire la courbe 6’ = f(Z) ; on retrouve les 
sauts d'écran, atténués toutefois, aux passages des 
électrons p, à p,. La constante 6’, plus faible 
que 6, ne convient pas non plus pour les calculs de 
diamagnétisme. 

Il semble donc logique de corriger les effets 
d'écran 6 pour les adapter à ces calculs. Pour cela 
trois facteurs vont intervenir au premier chef : la 
considération des effets d'écran qui accordent fina- 
lement à l’électron une influence supérieure à 
l'unité ; le respect de la structure de l’atome appa- 
raissant dans les sauts d’écran ; le fait qu’en son 
essence le diamagnétisme n’est pas quantifié. 

Revenons aux sauts d’écran en tant que gran- 
deur. Ils figurent au tableau IT. Sur 33, il n’y en à 

ue 8 qui soient nettement supérieurs à l'unité. 
Is correspondent aux transitions suivantes : 


2p3-3s : 3p3-4s ; 4sD-4p, ; 4p,-5s ; 5sD-5p, ; 5p,-6s: 
6SDF-Gp, ; Gpa-Ts. 

Toutes ces transitions correspondent au passage 
de fonctions d’onde à symétrie sphérique à des 
fonctions d’onde à symétrie différente, et quatre, 
en plus, correspondent à des accroissements 
brusques du nombre quantique d’une unité. Quant 
aux autres,certains présentent uninconvénient équi- 
valent. Prenons, par exemple, l’effet qui apparaît 
entre le dernier électron 3s et le premier élec- 
tron 3p.. Il ne vaut que 0,6066, mais comme il faut 
l’ajouter à chaque valeur apportée par les élec- 
trons 3p,, soit 0,6518, chacun de ces électrons 3p, 
contribuerait à un effet d’écran pour 1,2584, nette- 
ment supérieur à l’unité, donc, sans doute, à dimi- 
nuer ou même à repousser. 

Nous avons pris le parti de les repousser tous, à 
l’exception des sauts relatifs aux passages p; à pa 
pour lesquels le remplissage de trois cases quan- 
tiques doit augmenter singulièrement l'effet 
d’écran. Cette détermination qui amène à retran- 
cher les sauts d'écran correspondant à une augmen- 
tation du nombre quantique d’une unité ou à un 
changement de symétrie de la fonction d’onde fait 
qu’on passe des effets d'écran © aux effets cn plus 
faibles qui conviennent pour le calcul des coefli- 
cients diamagnétiques. Nous donnons leurs valeurs 
au tableau IV. Ce sont des effets d'écran semi-expé- 
rimentaux utilisables sans doute pour le diama- 
enétisme seulement. 

En effet, nous l’avons vu, en raison de la non 
quantification de cette propriété générale de la 
matière il faut diminuer l'influence de x dans le 
calcul de y. Slater prend les nombres quantiques 
effectifs n', nous nous contentons de ne pas tenir 
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TABLEAU IV 


NOTA 


V/ ATOME  ÜLECTRON Om Z ATOME  ÉLECTRON Om Z ATOME ÉLECTRON os 
2 He 15 0,6558 | 31 Ga pi 23,2296 | 60 Nd 6sd/ 48,05865 
3 ii 25 1,2578 | 32 Ge 4p: 23,7968 | 61 Pm 6sdf 46,95485 
ARBRE 2s 1,8598 | 33 As 4p; 24,3640 | 62 sm Gsdf 47,85105 
OUR 2p; 2,6071 | 34 Se 4P: 25,5010 3 Eu Gsdf 48,74795 
CC 2p\ 3,3544 | 35 Br hPa 26,166 64 Gd Gsdf 49,64345 
RE N 2p; 11017. | 36 Kr hPa 26,8310 | 65 Th Gsdf 50,53965 
8e) 2p» 51691 | 37 Rb 5s 27,3653 | 66 Dy. Gsdf 51,43585 
9 F 2P» 5,8419 | 38 Sr 5s 27,8996 | 67 Ho Gsdf (52,33205) 
10 Ne 2m 66021 so EN 5sd 28,7895 | 68 Er 6sdf (53,22825) 
11 Na 3? 7,2483 | 40 Lx 5sd 29,6794 | 69 Tm Gsdf (54,12445) 
12 Me 3s 7,8445 4 Nb 5sd’ 30,6172 70 Yp 6sdf 55,02065 
RAT 3p: 8.963 |42 Mo 5sd' 31,5550 | 721 Lu 6sdf (55,91685) 
14 Si 3P1 9,1481 | 43 Te 5sd’ 32,4928 | 72 JHf 6sdF (56,85825) 
15 P 3P: 9,7999 | 44 Ru 5sd’ 33,306 | 73 Ta 6sdF (57,79965) 
TO RTS 3P» 10,8709 | 45 Rh 5sd’ 34,8684 | 74 - W 6sdF 58,74105 
17 CI 3P: 11,5654 | 46 Pa d 36,0064 | 75 Re 6sdF (59,68245) 
TE 3P) 12,2599 | 47 Ag 5sD 36,6711 | 76 Os 6sdF 60,62385 
19 K 4s 12,8367 | 48 Cd 5sD 37,3358 4 [077 Ir 5d (61,56525) 
20 NCA 4s 13,135 | 49 In 5p; 37,98605 | 78 Pt 6sdF 62,50665 
21 Sc &sd 144125 | 50 Sn 5p: 386363 21-79 AU 6SDF 63,19165 
DAT 4sd 15,4115 | 51 Sb 5p; 39,28655 | 80 Hg 6sDF 63,87665 
GANT &sd 16,105 | 52 Te 5Po 40,18095 | 81 TI 6P; 64,5801 
AMC &s'd 17,2738 | 53 I 5P2 40,8537 | 82 Pb 6p: 65,28355 
25 Mn 4s'd 18,13710 | 54 Xe 5P2 41,52645 | 83 Bi 6P1 65,98700 
26 Fe &sd 19,1584 | 55 Cs 6s &2,01975 | 84 Po 6P» 66,85907 
27 Co 4sd 20,797 | 56 Ba 6s &2,51305 | 85 At 6P2 67,56254 
SRE Asd 21,:010 | 57 La 6sd 43,37005 | 86 Rn 6P» 68,26600 
29 Cu 4sD 21,8880 | 58 Ce 6sdf 44,26625 | 87 Fr 75 68,72100 
20070 4sD 225750 | 59 Pr 6sdf 45,16245 | 88 Ra 75 69,17600 


compte des accroissements de © en fonction de », ce 
qui revient à augmenter le dénominateur de y sans 
toucher au numérateur où figure n. 

D'autre part, en analysant les calculs de Slater 
on constate que, ne tenant aucun compte du nombre 
quantique {, tout se passe comme s’il admettait les 
orbites circulaires, et en utilisant le rayon de Van 
Vleck et Pauling, on utilise aussi des orbites cireu- 
laires. Il paraît difficile, avec les constantes © ou 6’ 
calculées à partir d’une orbite circulaire d’obtenir 
un résultat correct lorsque la fonction d’onde n’est 
pas de symétrie sphérique. Le caleul ne devrait 
convenir que pour les électrons s. Cependant dans 
le travail de Slater on constate que les effets d’écran 
des électrons s et p obéissent aux mêmes règles de 
calcul. C’est pourquoi les sauts d’écran relatifs à la 
transition s-p peuvent être supprimés comme ceux 
qui concernent la transition p-s pour une raison, 
cette fois, différente. Quant aux transitions p:-p2 qui 
ne présentent aucun saut d’écran dans le système 


de Slater, mais en présentent un dans nos déter- 
minations, soit avec le rayon de Bohr, soit avec 
celui de Pauling et Van Vleck, elles sont plus diffi- 
ciles à interpréter. Il n’y a aucune variation des 
nombres quantiques n et !, la symétrie des fonc- 
tions d’onde demeure, nous convenons de conserver 
ce saut d’écran de transition donnant ainsi aux 
cases quantiques un rôle qui nous échappe mais qui 
paraît certain. 
Nous utiliserons donc les valeurs oc. 


Caleul des coefficients d’aimantation diama- 
gnétiques des orbites électroniques. — Ce calcul 
découle de l’équation de Van Vleck [2] applicable à 
l’atome A : 

PU REN | 
TT. He 
C'est la vitesse de la lumière. 
L'énergie d’un électron sur l'orbite de Bohr ou de 
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TABLEAU V 


VALEURS DE À 


n l n? Ë n° en K ÉLECTRON 
1 0 3 2,379781 15 
2 0 42 33,31692 25 
il 30 23,79780 2p 
3 0 207 164,204 82 3s 
1 180 142,78680 3pP 
2 126 99,95076 3d 
F 0 648 514,03248 &s 
1 600 475,95600 4p 
2 504 399,80304 &d 
3 360 285 57360 kf 
5 0 4 575 1 249,38450 5s 
dl 1 500 1 189,89000 5p 
2 1 350 1 070,90100 5d 
6 0 3 258 2 584,44108 6s 
1 3 150 2 498,76900 Gp 
) 2 934 2 327,42%84 6d 
7 0 6 027 & 780,97802 75 
TABLEAU VI 
Z A K(1s) K(2s) K(2p) K(3s) K(3p) 
1 0,79469 2,38407 
2 358 074 
3 345 035 33,32490 
L 340 020 2280 
5 338 014 2196 
6 338 014 2196 23,8014 
7 339 005 2070 801% 
8 334 002 2028 8005 
9 328 2,37984 1776 8002 
10 328 84 1776 7984 
11 328 84 1776 7984 164,20896 
12 328 84 1776 7984 896 
13 327 81 1734 7981 689 142,7886 
mêmes valeurs jusqu’à Z == 18. 
Sommerfeld — la différence entre les deux est Mais le rayon de Bohr ne convient pas, le rayon 
minime et n'apparaît que dans les structures fines de Van Vleck et Pauling donne : 
— s'écrit — Ze?/2r et permet de poser avec NA: A nues 
1 Zion m Foie 96m C? e? r4(Z — s)? Gr 5 | 
tte RES E A TEE 
# mr 10 | 
Nef 2? n+h4 Nh ni 


X — — onmC? Ziebknim 96m C?e? ni 7? et À, avec les valeurs actuelles des constantes, 
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vaut 7,9326.1077. En remplacant s par notre 6m 
il vient alors l'équation (4) déjà donnée qui, plus 
simplement, peut s’écrire : 


K 
ZA tn 12 
fi (Z Fe Gm) ? ( ) 
avec : 
Leg 
K — 0,79326 Ë ni ets | (13) 


Le tableau V donne les valeurs de X pour chaque 
catégorie d’électron. Précisons que si l’on veut 
tenir compte de l’entraînement du noyau il faut 
utiliser la masse réduite de l’électron. Nous en 
avons tenu compte pour les 18 premiers éléments 
chez lesquels À et X ont les valeurs corrigées du 
tableau VI. 


Passage aux atomes libres. — On calculera la 
part contributive au diamagnétisme de chaque 
orbite et on fera la somme. Deux difficultés sur- 
gissent. 

19 L’expérience donne l’énergie de première 
ionisation des atomes. On peut connaître aisé- 
ment on et faire le calcul. Mais comment détermi- 
ner le 6, des autres électrons plus profonds ? 

On pourrait s’adresser aux énergies de seconde, 
troisième … ionisation, mais, si une première 10ni- 
sation ne bouscule pas trop l’édifice électronique, 
rien ne nous dit que la seconde, puis la troi- 
sième … etc. ne conduiraient pas à de grossières 
erreurs ! On prendra alors pour 6, celui qui résulte 
de la première ionisation de l’atome précédent, 
donc de numéro atomique Z — 1, puis Z — 2, etc... 
jusqu’à épuisement du nombre d'électrons conve- 
nable. Ce choix se justifie en prenant au hasard un 
atome, l’étain par exemple, pour lequel on connaît 
les ionisations successives jusqu’à concurrence de 
9 électrons enlevés. Si on compare l’effet d’écran 6 
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de ces ionisations successives calculées à ceux obte- 
nus sur les atomes immédiatement antérieurs In, 
Cd, Ag, Pd, on trouve : 

6 pour les cinq ionisations successives : 


46,3337 44,8329 42,506 41,4825 40,248 
o pour les quatre atomes précédant Sn : 
45,7433 43,9382 43,2735 48,8726. 


Ces dernières valeurs sont supérieures et l’écart 
va croissant avec le nombre d’électrons enlevés : 


0,9104 41,4322 1,7910 2,6246 


ce que l’on peut Htenprétes comme des pertur- 
bations croissantes dans l’édifice électronique pro- 
voquées par l’ionisation. 

20 En toute rigueur, il faudrait faire la somme de 
tous ces résultats, concernant tous les électrons. 
Mais on voit assez vite qu’après avoir considéré les 
huit électrons les plus externes des atomes de gaz 
inertes, les seuls libres, la neuvième orbite donne- 
rait un y largement négligeable, a fortiori les orbites 
suivantes, à partir de Z = 18. 

En général les y à considérer sont ceux des élec- 
trons dits de valence et quelques autres jusqu’à 
concurrence de huit semble-t-il, Pour le mercure, 
par exemple, sept électrons suffisent pour aboutir 
à — 33,558.1076 et Cohen [9] par une méthode 
expérimentale récente très précise a obtenu 
— 33,214.1076. 

Le choix du nombre d’électrons à considérer pré- 
sente une certaine difficulté. 

Le tableau VII donne quelques résultats déjà 
publiés ailleurs [6]. On trouve dans l’excellent 
ouvrage de Selwood [15] quelques valeurs calculées 
pour les gaz inertes. Le tableau VII montre les 
progrès réalisés dans ces calculs. Cependant les 
résultats obtenus avec les atomes qui ne sont sûre- 
ment pas libres ne concordent pas avec les valeurs 


TABLEAU VII 

ATOMES VALEURS EXPÉRIMENTALES CALCULÉES ANIONS Re ee 
He 1,88 [10] 1,91 HA) 1,913 Fe 41 11,919 
Ne 6,660) 65 UI11006,75[1 ie 7,167 Cle 26 25,943 
A 18,13 [10] 19,23 [12] 19,57 [2] 19,2 [13] 20,218 Br— 36 36,261 
Kr 28,02 [12] 29,2 [13] 30,005 1] 52 51,495 
Xe 42,20 [121004410119] 44,600 (*) (d’après Klemm [14]) 
Rn ? 7 49,872 

ATOMES AUTRES VALEURS CALCULÉES (D'APRÈS SELWOOD [15]) 
He 1,853 1,54 1,90 3,51 (Les premières valeurs sont celles de Slater.) 
Ne 527 D 8,6 13,78 
À 18,9 21,5 24,8 21,04 
Kr 31,7 42 33,38 
Xe 48 66 42,90 


NO 


expérimentales (molécules de corps simples poly- 
atomiques par exemple). Une étude en cours, assez 
avancée pour révéler sa fécondité, nous a déjà 
donné en très bon accord avec l’expérience les 
coefficients d’aimantation de H,, C combiné, F,, 
CB CIS BTE IH FH; H.0, HS NH. 
CICH3, BrCH3, ICH3:, CH;3-CH3,.. etc. mais il faut 
faire intervenir les énergies d’ionisation des molé- 
cules ou les énergies de formation et, dans le cas 
du carbone l’hybridation des orbitales, de sorte 
que le calcul s’écarte singulièrement des relations 
simples que l’on a pu admettre jusqu'ici. 

Le cas des atomes paramagnétiques est parti- 
culièrement intéressant. Ils possèdent, on le sait, 
un diamagnétisme de base beaucoup plus faible 
que le para, par suite masqué par ce dernier. Si l’on 
suppose que le paramagnétisme ne perturbe pas 
trop la disposition spatiale des électrons les calculs 
précédents gardent toute leur valeur. 
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C’est ainsi que la considération des électrons 4d 
et 4p du palladium conduit à un diamagnétisme 
de — 2820.10 6. Wucher [7] admet — 30.106. 
Pour le platine Kopp [8] admet — 27,33.10—6. Le 
calcul nous donne — 31,57.10—$6, Toutes restric- 
tions faites ces atomes n’étant pas libres mais liés 
dans le métal, les écarts sont très acceptables ! 

On pourrait finalement se demander si, en corri- 
geant les effets d’écran que nous avons nommés 0’ 
(déduits de W;i par l'intermédiaire du rayon de 
Van Vleck et Pauling) des sauts d’écran qui existent 
aussi dans ce cas, on n’aboutirait pas également à 
des valeurs acceptables pour le diamagnétisme des 
gaz inertes. La réponse est négative car, avec 6’, on 
aboutit pour l’argon à — 105 y — 34 et avec o’ 
(corrigé) à 14. La moyenne même, 24, n’est pas 
valable puisque l'expérience conduit à environ 19 ! 


Manuscrit reçu le 11 janvier 1960. 
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Un Colloque international sur « La physique des forces 
électrostatiques et leurs applications », se tiendra à Grenoble 
du 27 septembre au 127 octobre 1960. Le programme 
comprendra, en particulier, l'étude du triage électrostatique, 
des précipitations électrostatiques, des projections électrosta- 
tiques et des générateurs et moteurs électrostatiques. Pour les 
renseignements complémentaires etl’inscription, s’adresser 
au Professeur N.J. Felici, Secrétaire du Colloque, Institut 
Fourier, place du Doyen Gosse, Grenoble (Isère). 
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ÉNERGIE INTERNE ET CONSTANTES ÉLASTIQUES DU SILICIUM IRRADIÉ 
AUX NEUTRONS RAPIDES 


Par Guy MAYER et Marcez LECOMTE, 


Centre d'Études Nucléaires de Saclay. 


Résumé. — Des monocristaux de silicium ont été irradiés en pile jusqu’à la dose de 3,7.1019 neu- 
trons rapides par cm?. L'énergie interne et les propriétés élastiques ont été mesurées à diverses 
étapes d'irradiation. L'influence de traitements thermiques sur les mêmes propriétés du silicium 


irradié a ensuite été étudiée. 


À l’aide des valeurs numériques ainsi obtenues, nous essayons de préciser la nature et le nombre 
des défauts créés dans ces cristaux par les neutrons rapides. 


Abstract. — Silicon monocrystals were pile irradiated with doses of up to 3.7 X 101? fast neutrons 


per cm?. 
irradiation. 


The internal energy and the elastic properties were measured at various stages of the 
The influence of heat treatment on these properties was then investigated. 


Using the numerical values thus obtained an attempt is made to determine the nature and 
number of defects induced in these crystals by fast neutrons.- 


Dans ce travail nous avons utilisé des échan- 
tüllons monocristallins de silicium contenant moins 
de 10 ppm d’impuretés. 

Au cours des irradiations, effectuées dans la pile 
à eau lourde E.L.2 de Saclay, les échantillons, 
placés dans un tube d’uranium naturel où la tempé- 
rature était de 70 0C, ont été soumis à un flux 
d'environ 101? neutrons.s-!.cm ? pendant des 
durées variables. 

Dans ce chiffre de flux et dans les doses de neu- 
trons « rapides » que nous indiquerons, nous ne 
comptons que les neutrons d’énergie supérieure à 


1 MeV. 


A) Constantes élastiques. -— Nous avons obtenu 
les trois constantes élastiques indépendantes du 
silicium à l’aide de mesures de fréquences de réso- 
nance en modes longitudinaux ou transversaux, 
d'échantillons cylindriques diversement  orien- 
tés [1], [2]. Par exemple, sur un cylindre d’axe [100] 
nous calculons s,, et s,, à l’aide de deux mesures de 
fréquences de modes longitudinaux. 

De la même façon, nous pouvons obtenir les 
quantités 


(1/3) (S11 + 2512 + S44) 


(1/6) (2511 + 4S19 — Sy) 
à l’aide d’un cylindre d’axe [111]. 

L'étude des modes transversaux donne des véri- 
fications supplémentaires. Le tableau [ montre 
l’évolution de ces trois constantes en fonction de 
la dose © de neutrons rapides. 

La figure 1 montre l’évolution de la fréquence 
fondamentale de résonance longitudinale d’un 
cylindre d’axe [111](1— 35,68 mm 9 — 9,56 mm) 
en fonction des doses d'irradiation. A la dose 


et 


< 


de 3,5.101? neutrons.cm ?, la pente est au moins 
dix fois inférieure à la pente initiale. 


TABLEAU I 


EFFET DES NEUTRONS RAPIDES 
SUR LES CONSTANTES ÉLASTIQUES 


® S11 S12 S44 
0 76,70 — 21,24 125,38 
4,75 77,94 — 21,75 127,38 
3,7 78,12 — 21,95 128,20 


Les doses sont exprimées en unités de10#neutrons.cm? 
et les siÿ en 10—4 0, g.s. 


Les courbes décrivant l’évolution des fréquences 
propres de modes longitudinaux ou transversaux 
relatifs à des directions quelques peuvent se cal- 


_— 


[v. 10 7s: 


; : se 10 n.cm°? 


= Re es se 


Fic. 1. —- Effet des neutrons rapides sur la fréquence 
fondamentale longitudinale d’un cylindre de silicium 
d’axe [111]. 


N° 4 


culer à l’aide des chiffres du tableau I ; elles ont 
toutes même allure que la courbe de la figure 1. Il 
en est de même pour la courbe relative à l’énergie 
interne. 

Il est donc probable que la courbure de cette 
courbe ne provient pas d’un défaut de propor- 
tionnalité entre le nombre des défauts produits par 
les neutrons et les effets observés, mais plutôt du 
a que ce nombre n’est pas fonction linéaire de la 

ose. 
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Fic. 2. — Influence de la température sur une fréquence 
longitudinale d’un cylindre [100] (2 — 33,95 mm 
D — 10,40 mm) à diverses doses d'irradiation. 


Keesom, Lark-Horovitz et Pearlman [3] ont 
mesuré la chaleur spécifique du s'licium irradié et 
ont conclu à une diminution de la température de 
Debye qui correspond à une évolution moyenne du 
spectre de vibration vers les basses fréquences. 
C’est bien ce que nous observons à la limite infé- 
rieure de ce spectre. 

L'influence d’irradiations non isotropes de neu- 
trons rapides sur les vitesses d’ondes transversales 
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diversement orientées et polarisées a été observée 
par R. Truell, L. J. Teutonico et P. W. Lévy [4]. 

L’allongement linéaire dû aux irradiations est 
très faible ; il vaut 1,7.10—4 d’après les mesures 
aux rayons X de M. C. Wittels [5] après une dose 
de 4.10% neutrons.cm-?. D’après nos mesures 
macroscopiques cet allongement est de 2.104 
pour ® = 3,7.101° neutrons.cm?. 

Remarquons qu’à cette dose la variation relative 


SAS 


pe ——— 


EPS Re TEA 


L 1,260 


n282 
1,250 (g=1,7.10") 4 
Sd CS ce 
(p=3,7.10"°) 
RÉ NU ne 
Rte 
Et 
78 125 178 223 DESSIN 
Fic. 3. — Influence de la température sur une fréquence 


longitudinale d’un cylindre [111] à diverses doses d’irra- 
diation. 


des vitesses d’ondes acoustiques est de l’ordre 
de 100.104 et nous verrons plus loin (fig. 2 et 3) 
que la variation relative de la pente en température 
à 20 0C des vitesses d’ondes longitudinales est de 
l’ordre de 5 000.1074, 

Le frottement interne des modes longitudinaux 
se trouve, à 20 °C multiplié par un facteur 5 pour 
la dose de 3,7.101% neutrons.cm ?; mais à 
— 180 ©, le silicium irradié a le même frottement 
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interne très faible que le non irradié. Les modifi- 
cations provoquées par une irradiation à 70 °C sont 
stables indéfiniment à 20 °C et au-dessous. 

Nous avons mesuré en fonction de la température 
entre — {80 °C et 20 0C l’évolution des fréquences 
fondamentales . longitudinales de cylindres 
d’axe [100] et [111] (fig. 2 et 3). Dans les échan- 
tillons irradiés (même faiblement), il apparaît un 
maximum. 

Nos mesures permettent de prouver que ce maxi- 
mum existe pour tous les modes longitudinaux, 
mais ne suffisent pas à prouver qu'il en est de 
même pour les modes transversaux. 

On peut montrer de façon très générale, à l’aide 
du principe de Nernst, qu’au zéro absolu les déri- 
vées des constantes élastiques ds:;/dT doivent être 
nulles, donc aussi les pentes des courbes des 
figures 2 et 3. 

Ainsi, il doit exister à une température comprise 
entre 0 °K et la température du maximum de fré- 
quence un point d’inflexion dans les courbes rela- 
tives au silicium irradié. 

En l’absence de données expérimentales à des 
températures très inférieures à celles de l’azote 
liquide (McSkimin [6] s’est arrêté à — 195 0C), 
nous ne pouvons être sûrs que ce phénomène ne se 
produise pas aussi dans le silicium non irradié, 
l’effet de l’irradiation étant alors seulement de dé- 
placer vers les hautes température le maximum. 

Physiquement, il est intéressant d'écrire la quan- 
tité mesurée dv/d7T sous la forme suivante [7] : 


dv dv dv 

d TR 6 log 5% ), (1) 
où y désigne une fréquence propre, v le volume et & 
le coefficient de dilatation. L'existence du terme 
(dv/2 log v), est liée à celle de dérivées d’ordre 3 de 
l’énergie potentielle exprimée en fonction des défor- 
mations. (dv/d log v)r peut être mesurée direc- 
tement [8]. 

(dv/dT'), correspond aux dérivées d’ordre 4 de la 
même fonction. Ce sont les termes du type 
(dv/à log v)? qui sont responsables de la dilatation 
thermique «&. En effet : 


d log Vi 

K est la compressibilité isotherme, C(v) est la 
chaleur spécifique associée à l’oscillation de fré- 
quence v; et la somme est à effectuer sur tous les 
oscillateurs du solide. 

Dans le silicium non irradié, c’est le deuxième 
terme de l’équation (1) qui est en moyenne le plus 
important numériquement. D’autre part, ces deux 
termes peuvent n’être pas de même signe. 

Dans le silicium naturel, « petit à 20 °C 
(æ 6.1076.0C1), devient négatif vers — 160 0C, 
ce qui indique que les quantités (9 log v;:/à log v)» 
ne sont pas toutes de même signe. 


N° #% 


Nos mesures de coefficient de dilatation 
moyen œm entre —— 180 0C et O °C nous montrent 
que 4m varie peu par irradiation : avant irradiation 
am = 2,31.1076,0C—1 ; après une dose de 3,7.101? 
neutrons .cm£, %m — 2,46.107$ 0074, 

Il est donc probable que les quantités 
(à log v;/log v), sont peu modifiées par l’irradia- 
tion et que les effets observés sont dus à une forte 
diminution des termes (av: /dT),. 

+ Des mesures directes des quantités (dv;/0 log v)r 
dans des échantillons naturels et irradiés per- 
mettront de préciser ce point. 


B) Énergie interne. - Quand on chauffe un 
échantillon irradié à des températures supérieures 
à la température d'irradiation, il se produit une 
évolution irréversible de ses propriétés. Si on le 
chauffe suffisamment fort et longtemps, l’échan- 
tillon se retrouve, au terme de son évolution, dans 
le même état thermodynamiquement stable 
qu'avant irradiation. 

La méthode expérimentale de l’analyse ther- 
mique différentielle [2] montre que certaines phases 
de cette évolution sont exothermiques et permet de 
préciser la valeur numérique de la différence d’éner- 
gle interne entre silicium non irradié et silicium 
irradié plus ou moins « guéri » thermiquement. La 
figure 4 montre le résultat de telles expériences sur 


NT 
FO,S 
(B) (A) 
NS 9 
550 200 
Fic. 4. — Analyses thermiques différentielles du silicium 
irradié. 


Courbe À : o — 1,4.101° n .em2. 
Courbe B"9—="3;7.101%n "cm. 


deux échantillons irradiés respectivement aux 
doses de 1,4.10% et 3,7.102% neutrons .em 2. Ces 
expériences, menées jusqu’à 600 0C au rythme 
de 2,7 0G par minute, révèlent des phénomènes 
exothermiques entre 100 9C et 410 0€. 

Pour l’échantillon le moins irradié la première 
bosse correspond à 1,2 cal.g1 et la seconde à 
0,2cal.g=t 

Pour le plus irradié ces chiffres deviennent res- 
pectivement 2,1 cal.g7t et 0,3 cal.g71. Rappelons 
à titre de comparaison que d’après Berthelot [9] la 
différence d’énergie interne AL entre silicium 
amorphe et silicium cristallin vaut 289 eal.g1 : 
d’après Olette [10] entre silicium liquide et solide, 
au point de fusion AL = 432 cal.g71, L’irradia- 
tion abaisse fortement la conductivité thermique : 


_ Nou 


ce phénomène nous gêne dans l'interprétation 
numérique de ces analyse thermiques. 


(®) Description de la guérison thermique. -— 
L'évolution vers son état initial du silicium irradié 
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peut aussi être suivie par des mesures de fréquences 
plus précises que les opérations calorimétriques. 
La figure 5 est relative à la fréquence fondamen- 
tale longitudinale d’échantillons orientés sui- 
vant [100]. La courbe inférieure donne, en fonction 


ER 


T T T T 


400 500 600 


ÉCRIN 


—_—_— 


[1,51 100 200 300 


F16. 5. — Mesures d’une fréquence longitudinale [100] à température croissante. 
1) —@— Silicium non irradié. 2) —O— Silicium irradié à 3,7.104 n.cm?. 


de la température, l’évolution de la fréquence d’un 
échantillon irradié à la dose de 3,7.101° neutrons. 
em? quand on le chauffe au rythme de 2,1 0C par 
minute ; elle ne peut être obtenue qu’une fois sur 
chaque échantillon irradié puisqu'elle révèle des 
phénomènes irréversibles. La courbe supérieure est 
relative à un échantillon non irradié ou complè- 
tement guéri, chauffé de la même façon. 

Une figure tout à fait analogue représenterait les 
résultats relatifs à un cylindre d’axe [111]. 

Entre 120 °C et 300 0C on assiste ainsi à une 
évolution rapide de la fréquence en direction de sa 


valeur d’avant irradiation. Pendant cette évolution 
une augmentation de frottement interne se mani- 
feste, qui disparaît vers 300 0C. 

Si on se reporte à la figure 4 on voit que ce phé- 
nomène coïncide avec la principale « bosse » exo- 
thermique de la courbe d’analyse thermique obte- 
nue avec le rythme de chauffe un peu différent de 
2,7 0C.mn 1. Par contre, la petite bosse exother- 
mique dont le maximum est à 345 0C ne correspond 
à aucun accident de la courbe de fréquences. 
Au delà de 300 9C, la courbe de fréquences du sili- 
cium irradié est encore distincte de celle du silicium 
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naturel ; c’est vers 600 °C que se produit l’évolution 
du silicium irradié qui les amènera à se confondre. 

En résumé, nos expériences montrent qu'il se 
produit entre 120 °C et 300 °C une évolution qui 
affecte fortement et simultanément l’énergie interne 
et les constantes élastiques du silicium irradié. Elle 
traduit la transformation ou la disparition pro- 
gressives de défauts dus à l’irradiation. 

Dans l’hypothèse où il s’agit d’une disparition, 
nous allons essayer de définir à l’aide de nos me- 
sures un paramètre Ë, proportionnel au nombre de 
ces défauts, qui nous permettra d’en étudier la 
cinétique. 

Le rapport BA/CA défini sur la figure 5 décroit 
rapidement entre 120 0C et 300 0C ; les points C 
sont portés par la courbe en pointillé obtenue par 
extrapolation de la courbe de fréquenec du silicium 
irradié, depuis le domaine où il est stable. BA /CA 
encore égal à l’unité vers 100 0C, décroît jusqu’à 
0,25, valeur qu’il garde au-dessus de 300 0C Le 
point A’ défini à chaque température par la rela- 
tion A'A/CA = 0,25 décrit la courbe réelle de fré- 
quences d’un échantillon irradié, puis guéri à une 
température comprise entre 300 0C et 500 °C 


Nous poserons & — BA’/CA'’ Aïnsi Ë décroiît de 1 
à O entre 120 0C et 300 0C. La figure 6 permet de 
1,00 ; ] 
0,75 | 
r 0,50 ; 
+0,25 | 
) 
—— 2 
0 1 
NES Re 0 sec 
F1G, 6, — Fonction£ (7) décrivant la guérison dusilicium, 


1) & obtenu sur un cylindre [100]. 
2) 6 obtenu sur un cylindre [111?. 


comparer les valeurs E(T) obtenues à l’aide des 
fréquences de deux échantillons différemment 
orientés. 

Les mesures d'énergie interne conduisent aussi, 
avec une précision moindre, à la même fonc- 
tion E(T). 

Nous supposons donc désormais Ë proportion- 


N°4 


nelle au nombre » des défauts dont la disparition 
modifie les propriétés du silicium dans le donaine 
120 °C-300 0C. Cherchons si une loi cinétique simple 
du type 

dnfdt — — K(T) n* (3) 
où X(1') est de la forme A e— "7 et À, &, w sont 
des constantes, peut expliquer la forme de E(T). 


1) Essar D’UNE Lor DU 17 ORDRE. — On aurait 


alors 
dnfdt = —Xxe-"#Tn (4) 


Dans nos expériences, ‘la température T est 
fonction linéaire du temps f : 

TRE EU (5) 

dn]de — v(dnfdaT') (6) 


v représente le rythme de montée en température. 
L’équation (4) peut alors s’écrire sans la variable f : 


An [AT — — (fo) e—%ÀT.n (7) 
d?n À s dn nw 
in NID fee PRE 

EVE v° Gr F mr (8) 


Aïnsi dans l’hypothèse d’une loi du 1° ordre on 
aura au point d’inflexion de E(T), situé vers 200 °C : 


1 dé w 
Hide Êl 


Ce qui conduit à æ = 0,55 eV. 

On trouve alors pour À la valeur 5,6.102.s71, 
absolument inacceptable pour un « facteur de fré- 
quence » en phase solide [26]. 

Un tel facteur est en effet généralement [11] de 
l’ordre de 1012:s 74; 


2) Essar D’UNE Lor pu 2° ORDRE. — On aurait 
alors 
dnfdt = — xe "TT n? (10) 
dn [AT = — (x]v) e— "AT n2 (11) 
d?n À dn 02 
— — — pg—Ww/kT en Le 
ar? p° de (2 dr ne ms (2 
On aura au point d’inflexion de la courbe E(T) : 
2 dé sw : 
É AO: Ve 


On trouve 4 — 1,09 e.V. 

On peut, en désignant par n, le nombre initial 
de défauts, calculer An, à l’aide des grandeurs Ë et 
défdT prises à une température quelconque ; en 
effet l’équation (10) permet d’écrire : 


1 dé 
An = D CUET =. 1 
0 G ( EaT (1 1) 
On trouve : 
à 150 0C A6 29,0 108 
à 200 °C NT -RO LUS 
à 250 °C No == JE 
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An, peut être considérée comme constante et l’équa- 
tion (10) est correctement vérifiée. 


D) Signification d’une loi cinétique du 2€ ordre. 
— Bien que réellement ignorant de la nature des 
défauts dont nous suivons la disparition, nous 
allons essayer de justifier l’équation (10) et les 
valeurs trouvées pour An, et w. 

Nous avons longuement et vainement cherché 
dans des solides irradiés très divers comme les 
verres [12], le graphite [13], le quartz, l’alumine et 
la glucine [14], des lois cinétiques décrivant correc- 
tement les processus de guérison ; ceux-ci ont 
toujours refusé de suivre une loi simple et les 
quantités &, À, w qui figurent dans l'équation (3) 
refusaient de rester constantes. Aussi est-ce avec 
émotion et prudence que nous voyons dans le sili- 
cium se manifester un phénomène se conformant 
correctement à une loi du 2€ ordre. 

Si l’acte essentiel de la destruction ou de la 
transformation d’un défaut est sa rencontre avec 
un autre défaut, et si les deux défauts sont présents 
en égale quantité appelée n, la quantité dn/dt sera 
proportionnelle à #?. On peut imaginer qu'il 
s’agit d’atomes interstitiels et des lacunes qu’ils ont 
laissées derrière eux, ou bien uniquement de 
lacunes, etc. 

S'il s’agit de deux types de défauts différents, il 
est probable qu’une des deux espèces sera plus 
mobile que l’autre. 

D’après diverses théories [11] relatives au mou- 
vement des défauts dans un réseau solide, chaque 
défaut effectue en moyenne chaque seconde un 
nombre de sauts égal à 


SN er PR (15) 


où y est une fréquence de l’ordre de 102.871 
comme celle des oscillations atomiques, et w est 
une « énergie d’activation » qui mesure la hauteur 
de la barrière de potentiel que le défaut doit fran- 
chir pour passer d’une position de potentiel mini- 
mum à une autre. 

Supposons désormais pour fixer un langage qu’il 
s’agisse de lacunes mobiles diffusant à la recherche 
d’interstitiels avec lesquels elles s’anmhileront. 
Chaque seconde les lacunes effectuent vn e "#1 
sauts. La probabilité que chacun de ces sauts amène 
une lacune au contact d’un interstitiel est de l’ordre 
de la concentration n/N de ceux-ci. Ainsi : 


(dnfdé) = — ve-wKkTin?]N). (16) 


En identifiant avec nos résultats numériques, 
on voit que vnÿ{N = 7.108s71. En assignant à v 
la valeur 1012 $s—1, on trouve une concentration 
initiale de défauts égale à 7.107. 

Le raisonnement à l’aide duquel nous avons 
établi l’équation (16) est trop schématique. De 
nombreux perfectionnements lui ont été apportés, 
notamment par T. R. Waite [15], H. Reiss [16], 
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Fletcher et Brown [17] ; mais l’ordre de grandeur 
du facteur qui multiplie n?e—"#T dans l’expres- 
sion de dn/dit reste le même dans toutes ces théories. 

Par exemple, Fletcher et Brown arrivent à 
propos du germanium, en tenant compte de sa 
structure cristalline qui est celle du silicium et du 
diamant, à l’expression : 


dnfdt — — ve—"ËT(n2[N) (&rré bNa) (17) 


où r désigne le « rayon de capture » de l’inters- 
titiel que ces auteurs estiment voisin de la maille & 
du réseau ; b vaut a /3/16 et représente la valeur 
moyenne de la distance que franchit à chaque saut 
la lacune en direction de l’interstitiel qui l’anni- 
hilera. N, est le nombre d’atomes par em°. Le 
facteur 4rré bN, qui distingue l’équation (17) de 
l'équation (16) est bien de l’ordre de l’unité. 

Nous pensons que la vitesse de la réaction que 
nous observons dans le silicium irradié est contrôlée 
par un phénomène de diffusion. 

Fletcher et Brown [17] ont cherché à interpréter 
de cette façon la guérison du germanium irradié à 
la température ordinaire. La comparaison des 
points de fusion (respectivement 1209 0K et 
1 685 °K) et des températures de Debye (400 0K 
et 650 0K) dans Ge et Si laisse prévoir que des 
effets analogues à ceux qu’on observe dans Si se 
retrouveront dans Ge à des températures réduites 
d’un facteur voisin de 2/3. 

Ainsi il faudrait irradier Ge à la température de 
— 50 0C, et le phénomène du 2€ ordre que nous 
observons à partir de 120 °C dans Si commencerait 
à — 10 °C dans Ge. En irradiant à 60 °C du ger- 
manium à la dose de 1,2.101? neutrons .cm? nous 
avons observé une baisse relative de fréquence 
longitudinale, dans la direction [100], égale à 
0,75.10 3. Suivant la règle de correspondance des 
températures, il faudrait comparer ce chiffre à 
l'effet d’une irradiation du silicium à 220 °C. 

Il est difficile de savoir quel type de défauts 
diffuse dans le silicium avec l’énergie d’activation 
æ — 1,09 e.V. que nous avons trouvée. 

Bemski [18] et Ross [19] ont observé à l’aide de 
mesures électriques la création de défauts par 
trempe du silicium, puis leur guérison thermique. 
I1 semble que ces défauts liés à la présence d’impu- 
retés [19] ne soient pas de même nature que les 
défauts dûs aux irradiations. 


E) Évaluation de la concentration initiale des 
défauts. — En fixant à 1012s—1 1e facteur v, l'étude 
cinétique donne la valeur 7.10% à la concen- 
tration initiale des défauts qui disparaissent avant 
300 0C. Nous allons essayer de contrôler ce chiffre 
par d’autres méthodes. 


1) Mesures ÉLecrriQuEs. — Wertheim [20] et 
Sonder [21] ont conclu à partir de mesures de con- 
centrations et de durée de vie des porteurs de 
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charge, qu’au début d’une irradiation faite à 30 °C, 
une dose de ® neutrons.cm? créait en moyenne 
dans le silicium un nombre n de défauts égal à 5 ®. 
Les sections efficaces du noyau Si pour les neu- 
trons rapides étant connues, on peut calculer qu’un 
neutron sur 7 cogne un noyau Si dans un cm de 
silicium. Chaque choc crée donc, au début de l’irra- 
diation, 35 défauts. Ces résultats ne nous four- 
nissent que la pente à l’origine de la courbe n(). 
Nous n’avons pas fait de mesures analogues à 
celles de ces auteurs ; nous avons seulement vérifié 
le fait connu qu’au delà de © 101 neutrons .cm ? 
la résistivité de Si est de l’ordre de 106 ohm.cm |; 
aussi n’est-il pas certain que nous comptons nos 
neutrons rapides de la même façon que ces auteurs. 
En supposant que la pente de la courbe n(œ) 
décroît comme la pente de la courbe de constante 
élastique (fig. 1), on trouve qu’à notre dose finale 
de 3,7.10% neutrons.em?, 1l y a 6.10% défauts 
par cm, ce qui correspond à une concentration 


del2=10e, 


2) ÉNERGIE INTERNE. — Appelons DU l’énergie 
de recombinaison d’une paire de défauts et AL la 
valeur mesurée de la variation d’énergie interne au 
cours du processus de guérison, le nombre de 
défauts de chaque type sera égal au quo- 
tient AUL/SUL. En attribuant à SU (qui est sans 
relation avec l’énergie d’activation w — 1,09 e.V. 
qui préside à la mobilité) la valeur arbitraire 
de 2 e.V, on trouve une concentration de défauts 
de 13.104, 

Nous ne savons pas calculer SL, mais nous pou- 
vons lui assigner une limite inférieure ; en effet, 
si dUL était de l’ordre de 0,5 e.V par exemple, les 
défauts que nous étudions se produiraient sponta- 
nément [22] en quantité non négligeable dans le 
silicium solide ordinaire suffisamment chauffé, et 
ceci se traduirait par des anomalies dans la chaleur 
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spécifique. Or les mesures d’Olette [10] menées 
jusqu’à la température de fusion ne révèlent rien de 
tel. 


3) CONSTANTES ÉLASTIQUES. — Nos travaux sur 
la chaleur spécifique du graphite irradié, sur les 
constantes élastiques de BeO, AI?20%, nous ont 
montré qu’en général l’irradiation d’un cristal dé- 
place vers les basses fréquences son spectre de 
vibration. 

Si on suppose que les défauts produits par l’irra- 
diation sont surtout des paires vacances-inters- 
titiels, U faut conclure que l’effet « amollissant » des 
lacunes l’emporte sur l'effet « durcissant » des 
interstitiels. Ces deux effets ont été calculés par 
Dienes [23] dans le cas d’autres structures que celle 
du silicium. 

Dans l’alumine où le nombre absolu de défauts a 
été déterminé par une excellente méthode de dif- 
fraction de neutrons [24], on trouve une baisse 
relative moyenne de fréquences [14] de 5.107 pour 
une concentration de défauts de 1073(9 —53,3.10% 
neutrons .em?). Dans le silicium on crée d’après les 
raisonnements indiqués plus haut une concen- 
tration de défauts voisine de 107%. 

La grande sensibilité des constantes élastiques à 
la présence de défauts est un effet très curieux en 
lui-même [25] et le glissement vers les basses fré- 
quences du spectre de vibration qui est un autre 
aspect du même phénomène a une conséquence 
importante : l’entropie vibrationnelle du solide s’en 
trouve accrue et cela contribue à la stabilité des 
défauts. 


Nous tenons à remercier M. J. Guéron sous la 
direction de qui cette étude a été entreprise et 
M. H. Curien qui nous a constamment aidé de ses 
conseils. 


Manuscrit reçu le 6 février 4960. 
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I. EFFETS MAGNÉTOÉLECTRIQUES ET THERMOMAGNÉTOÉLECTRIQUES 
DANS LES SEMI-CONDUCTEURS. 


Par L. GODEFROY et J. TAVERNIER, 


Laboratoire Central des Industries Électriques à Fontenay-aux-Roses. 


. Résumé. On étudie la conductivité électrique d’un cristal en présence d’un champ magné- 
tique. La méthode utilisée est celle du gain moyen d'énergie. Elle permet d’exprimer sous forme 
tensorielle le vecteur densité de courant en fonction du champ électrique et du champ magnétique 


dans le cas d’une vallée unique. 


. Ce calcul permet ensuite d’obtenir, dans la limite des champs magnétiques faibles, les tenseurs de 
résistivité et de pouvoir thermoélectrique (ce dernier englobe le pouvoir thermoélectrique pro- 
prement dit, l'effet Nernst et la thermomagnétorésistance). Ces résultats sont utilisés dans les 
modèles particuliers généralement admis pour le germanium et le silicium. 


Abstract. — In this paper, the electrical conductivity of a crystal in the presence of a magnetic 


field is investigated, by the method of the average energy gain. 


It has been possible to express 


the current density in tensor from, as a function of the electric and magnetic fields, in the case 


of a single valley. 


Also, in the limit of low magnetic fields, the resistivity and thermoelectric power tensors 
have been obtained (the latter includes the thermoelectric power proper, the Nernst and thermo- 


magnetoresistance effects). 
silicon. 


L'étude systématique des coefficients de magné- 
torésistance électrique a souvent été utilisée comme 
moyen d'approche pour déterminer la structure de 
bandes des semi-conducteurs [1 à 6]. 

Le principe de cette détermination est de se fixer 
a priori un modèle de structure de bandes, de cal- 
culer l’expression du courant électrique en pré- 
sence d’un champ magnétique puis de la comparer 
aux résultats expérimentaux. 

Les modèles utilisés sont constitués par des sur- 
faces d’énergie constante de forme ellipsoïdale ou 
sphérique réparties suivant des axes privilégiés de 
la zone de Brillouin. Le calcul du courant est fait à 
partir de l’équation de Boltzmann en supposant que 
les champs électrique et magnétique sont faibles [7]. 

Bien que la symétrie cristalline permette de 
prévoir [8] la forme des relations entre les vec- 
teurs J (densité de courant), E (champ électrique) 
et H (champ magnétique), les calculs restent longs 
et ne permettent pas de relier simplement les diffé- 
rents coefficients de magnétorésistance aux valeurs 
du temps de relaxation et des masses effectives. 
Pour pallier ces inconvénients et formuler les résul- 
tats dans des cas plus généraux, nous avons appli- 
qué à ce problème la méthode du « gain moyen 
d'énergie » décrite par Aigrain et Englert [9]. Ceci 
nous permet d'obtenir, sous forme vectorielle, 
l'expression de la densité de courant pour une vallée 
ellipsoïdale unique orientée de façon arbitraire. 
Aucune hypothèse n’est faite sur l’intensité du 
champ magnétique. 

Il suffit ensuite pour traiter un modèle parti- 
culier de faire la somme des contributions des diffé- 
rentes vallées. On obtient ainsi l'expression du cou- 


The results are applied to the particular models of germanium and 


rant en fonction des champs électrique et magné- 
tique pour l’ensemble du cristal. Les coefficients 
des différents termes sont reliés explicitement à cer- 
taines valeurs moyennes du temps de relaxation. 
Il est possible ensuite, dans la limite des champs 
magnétiques faibles (2° ordre en H), d’inverser la 
relation précédente et d'obtenir l’expression vecto- 
rielle du champ électrique E correspondant au cou- 
rant «J. 

Une modification simple des résultats précédents 
permet d'obtenir sans nouveau calcul, Pexpression 
du fiux d’énergie isotherme en fonction du champ 
électrique et en présence d’un champ magnétique. 
En utilisant les relations d’Onsager, on en déduit le 
tenseur de pouvoir thermoélectrique, c’est-à-dire, 
les coefficients du pouvoir thermoélectrique clas- 
sique, de l’effet Nernst et de la magnétothermo- 
conductance. De plus, cette méthode permet de 
calculer pour toutes les valeurs de Æ la conduc- 
tivité en courant alternatif d’où l’on déduit les 
résultats relatifs à la résonance cyclotron, et à 
l'effet Faraday. Ces dernières applications seront 
publiées ultérieurement. 


A. Effets magnétoélectriques, 


I. Introduction. — 1. Hyporaèses. — Pour 
effectuer le calcul de la densité de courant J pro- 
duite par un champ électrique Æ en présence d’un 
champ magnétique H, nous nous placerons dans les 
conditions suivantes : 

a) Les surfaces d'énergie constante sont formées 
d'un ensemble de « vallées » elliptiques réparties 
suivant des axes simples du réseau réciproque. 
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Nous ferons d’abord le calcul pour une vallée 
unique (un ellipsoïde d’axes inégaux) orientée de 
façon arbitraire par rapport aux champs Æ£ et H, 
puis nous sommerons les contributions des diffé- 
rentes vallées au vecteur densité de courant pour le 
modèle retenu. Cette méthode implique lPhypo- 
thèse que les transitions subies par les porteurs ne 
se font pas entre vallées différentes, c’est-à-dire que 
les extrémités des vecteurs d'ondes k ne quittent 
pas la vallée qu’ils occupent initialement. 

b) Les transitions subies par les porteurs qui 
comme nous venons de l’admettre sont intravallées, 
conservent l’énergie. 

<) Il existe un temps de relaxation + ne dépen- 
dant que de l’énergie du porteur considéré, pour 
décrire l’évolution du vecteur d’onde moyen de ce 
porteur en fonction du temps et en l’absence de 
champs extérieurs. Le temps de relaxation sera 
supposé inchangé en présence d’un champ magné- 
tique (valable pour &, 7 € 1 61 w est la pulsation 
cyclotron). 

d) Le champ électrique est supposé infiniment 
petit, mais aucune hypothèse n’est faite sur le 
champ magnétique. 


2. PRINCIPE DU CALCUL. — Pour calculer la den- 
sité du courant produit par un ensemble d’électrons 
(ou de trous), nous utilisons la méthode « du gain 
moyen d'énergie » telle qu’elle est décrite par 
P. Aigrain et F. Englert [9]. 

Supposons que les champs E et H sont appliqués 
depuis le temps ? = — co, il en résulte que le 
régime permanent est établi quand { — 0, et consi- 
dérons un électron dont le vecteur d’onde est k, à 
l'instant { — 0. Nous voulons calculer la contri- 
bution que cet électron apporte à la densité de 
courant (évaluée au même instant { — 0). 

Depuis que les champs sont appliqués, le champ 
électrique a fourni à cet électron une énergie que 
l’on peut calculer par la formule : 


0 
a8——ef E<V(E,H,t:k) > dt. 
—Cco 


Où —e représente . la charge de l’électron, 
< V(E, H,t;k,) > représente la vitesse moyenne 
de l’électron considéré à l’instant #. Cette vitesse 
moyenne dépend de k, parce que nous considérons 
un électron qui aura le vecteur d’onde k, à l’ins- 
tant t = 0. 

Dans l’approximation du champ électrique infi- 
niment petit, au premier ordre en Æ, le gain d’éner- 
gie s’écrira : 

0 
AË =—ef" E< VIE = 0,H,1: ko) > dt. 
Les équations de la mécanique étant inchangées 


pour un renversement du temps si l’on remplace H 
par — H, nous aurons donc : 


+00 
JV = — ce] BI 
0 


VIE = 0, — H,1;k,) > de. 
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Il nous faudra alors intégrer les équations du 

mouvement pour obtenir 
<V(E =D; HR 
et pouvoir calculer A6. 

Connaissant le gain moyen d'énergie A6 des 
électrons, nous pourrons calculer leur nouvelle 
répartition dans les différents états d'énergie. En 
effet, soit (6) la fonction de répartition en l’ab- 
sence de champ extérieur (f(6) est la fonction de 
Boltzmann si le gaz d’électrons n’est pas dégénéré 
ou Ja fonction de Fermi-Dirac dans le cas où le gaz 
est dégénéré). ‘ 

En présence des champs extérieurs, le thermostat 
n'aura fourni aux électrons d’énergie & que l’éner- 
gie 6 — A6. Pour obtenir la nouvelle répartition 
dans les états d’énergie, nous remplacerons donc la 
fonction f(8) par f(6 — A6), que nous remplaçons 
par le début de son développement : 


1(6) — A6 f'(6) 
Nous pouvons maintenant calculer la densité de 


courant correspondant à une concentration de x 
électrons par unité de volume, soit : 


— ne | Vf(6 — A6) dk 
a) 


fre _—_ Ag) d'Æ 
où Vest la vitesse et A6 le gain d’énergie de l’élec- 
tron de vecteur d’onde k ; les intégrales sont éten- 


dues à tout l’espace des k. 
Au premier ordre en Æ, la densité de courant est : 


Fa ne [V ABf/(E) d® x 


J — 
fre) ax 


car l'intégrale se Vf(é) dk est nulle, parce qu’en 


l’absence de champ appliqué (A6 = 0) on doit 
avoir J = (. 


II. Caleul pour une vallée unique. — 1. Déri- 
NITIONS. — Nous allons appliquer au cas particulier 
de la vallée unique les calculs généraux exposés au 
paragraphe précédent. 

Nous admettons que pour la vallée considérée, 
l’énergie est donnée en fonction du vecteur d’onde k 
par la forme quadratique définie positive : 


De = D oi; ki k; 
> D 
à . ee r , . . 
où la matrice à des a; est réelle, symétrique et indé- 
pendante de k. 
Cette énergie peut s’écrire sous forme tensorielle : 


— 
ue 
DEN AC: 
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La vitesse d’un électron de vecteur d’onde k est 
donnée par : 
— 


<= 
VAE (fa. 


Dans la suite de ce paragraphe, nous nous pla- 
çons dans un système d'unités dans lequel % = 1, 
e — 1 et où l'unité de champ magnétique est telle 
que l’équation du mouvement d’un électron en pré- 
sence du champ magnétique H s'écrit : 


— 
— 
dkfd= MA H—=— Gk\ H. 


Nous reviendrons à la fin du paragraphe aux for- 
mules valables dans tous les systèmes d'unités. 


2. DÉFINITIONS ET PROPRIÉTÉS DE L'OPÉRA- 
TEURÊ.— Poursimplifier l’exposé,nousintroduisons 
> ce => 
as A Mt > 
l’opérateur linéaire 8 défini par 8 k — « k /\ H. 
En se plaçant dans le système d’axes de symé- 
trie de l’ellipsoide défini par 8 — Cte nous mon- 
trerons en annexe que : 


>> > 
> + 1 + 
x B k — — Fe ee 
Ta 
et 
> > > 
+ _ 
GE — F F OPA ei] \ H 
> > 
( > | 
= (1/2) | a —! 7 (I) Aa _ 
> > 
> à > ; 
où & est l’opérateur inverse de & et P représente 
= 
—> 


le déterminant de la matrice & 


Par la suite nous poserons : 
= 

A > ] 

Oe = — Es Gomd TI 


qui est une quantité positive (& réel) puisque la 
(a, 3-1) asc 
forme quadratique (H. «71 H} associée à la ma- 
trice «1 est définie positive. 
Du résultat précédent, on déduit immédiatement 
que: 


> > Éz 
> - Es ) 9 
B2k — 6 \B?k} = — ob k. 


Cette équation valable quel que soit k nous per- 
met d’écrire l’équation opératorielle vérifiée par £ : 
> > 
+ >> 
B°+ocp — 0, 
d’où nous déduisons que les valeurs propres de 
l’opérateur sont 0, roc et — 1e. 
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> > 


+ > 
Dans la base où 6 est diagonale l’opérateur e É! 

or : 

s'écrit : 


| 1 0 0 
= 
(\ EN 0 eioct 0 
0 0 e—t®ct 
ES 


donc l’opérateur e Ét appliqué à k, fait tourner ce 
dernier d’un angle «. { autour de l’axe défini par 
le vecteur propre correspondant à la valeur 
propre 1. Il est facile de montrer que cet axe est 
porté par H ce qui définit bien le mouvement d’un 
électron dans un champ magnétique. 


EYY 


Ceci nous montre que l'opérateur e F! 
bornée lorsque { — + co. 

Nous verrons que pour évaluer explicitement le 
gain moyen d'énergie A6, il nous faut calculer 

> > 

Rares ARS 
l'opérateur [| 1 — + 8 [. 

Nous montrons en annexe (2) que : 


aune norme 


se Et, RS ee > 

É +) # z Le 2 a 

1 ==710 = Ô - 

de Deer Die 

3. ÉQUATION DU MOUVEMENT. — L’équation du 
mouvement de l’électron s'écrit en représentation 

Ms die 
d’Heisenberg : re il, SC] 

* C 


1e = de . JE Le Joint 


est lhamiltonien total pour £ — 0 (approximation 

énoncée au paragraphel) 

avec J, — hamiltonien de l’électron dans le 
cristal parfait. 

— interaction avec Je champ magné- 
tique. 

JCint — interactions avec les perturbations 
produisant les transitions de 
l’électron représentées par le 
temps de relaxation +(6). 

Sachant que [k, X,] = 0 : mouvement d’un 
électron dans un cristal parfait et que — 1[k, Ju] 
> ue 
= — ak /\ H= — $ k: mouvement de l’élec- 
tron dans un champ magnétique HE, équation du 
mouvement s'écrit : 


JC 


ES 
die 4 


RE b EE uk, JE 


La valeur moyenne du vecteur d’onde k de 
l’électron vérifie : 


rs Ee : 
ee 6 CL ere [k, JCint] À 
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le terme — : <[k, Ji] > représente l’évolution 
de < k > pour les interactions de l’électron avec 
les perturbations du réseau cristallin. L'hypothèse 
du temps de relaxation pour le vecteur d’onde 
moyen consiste à écrire 
Ke 4 NUE 
US (Re, JCint] AE 


T(6) 


D'où l’équation du mouvement déterminant 
l’évolution de < k > : 


ST Pos ce ï L 
d< k> LR DNS RE 
RE Re 
de T 
4, CALCUL DU GAIN D'ÉNERGIE MOYEN. — En 
— 


—> 
tenant compte de la relation V — «x k et de l’hypo- 
_=S 


me 
thèse que « est indépendant de k, on écrit le gain 
moyen d'énergie A6 (dans le système d’unités 
choisi) sous la forme : 
= 
+00 > - 
Ag = — Eu IR 0 De Ed 
Comme le temps de relaxation ne dépend que de 
l’énergie de l’électron et que cette énergie est sup- 


posée rester constante, l’intégration de l’équation 
du mouvement (dans le champ — H) conduit à : 


SS 
_ _ 
T2 1 
<< R(E — 0, — H, t À k,) > — exp B — -Ë) he. 


Ce résultat nous permet de calculer le gain 
d’énergie sous la forme : 


+7 —1 
>> — 
+|+ 1 
uen exp. 


> 
car nous avons montré que exp 8 £.k, reste borné 
en module pour £ — + co. 

D'où : 


AG = — BC, Re ae à) 


< 

En utilisant le fait que « est un opérateur symé- 
trique, nous pouvons écrire l’expression de A& 
sous la forme 


Ë ) FAR ES 

Ie = alé .\a E - Ke he ) 

< to Irene RENE CE ATE 
T° 1 


= 
> 
RE RE 2 p bi (E.H) — (OP ie 2 IR 
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Dans la suite, nous appellerons k au lieu de k, 
le vecteur d’onde de l’électron considéré. 


5. CALCUL DU VECTEUR DENSITÉ DE COURANT. — 
A fin de simplifier les calculs, nous démontrerons en 


annexe (3) la formule suivante : 


JE ‘4 (Re. U] gl6) f'(8) dk nn 6? 8(6) f'(6) dé 
rer ax F fete) a 


dans laquelle U est un vecteur indépendant de k 
et £(8) /(6), deux fonctions quelconque de 6. 

Si de plus, nous remarquons que les fonctions de 
distribution (8) sont telles que 61/2 f(6) — 0 pour 
6 — oo, et que nous définissons la « valeur 
moyenne » < g > de g(6) par la formule 


1 8/2 g(6) f'(6) 46 
de 7 EU j/(8) dé 


0 


CRE 


o 


nous pouvons, en intégrant par parties, transformer 
la relation précédente en 


he x) [e. U] g(6) (6) dk 


rte) ax 


Tenant compte de ces résultats, le vecteur den- 
sité de courant s’écrit 


UE 


> es 
> 
MR CD - ( 


ne L (425) 
Dans le système d'unités C. G. S. l'expression 
précédente du vecteur densité de courant s’écrira : 


< 2 
2 3 ES 
J= nr «EE r - S (22 | 
ne œ n : 1 > pe \% IH \E 


3 
TP) QE 
1 EE o2 7? 


ce? P 


> 
e Ë (H.E) H — 2 2e 


où e est la charge de l’électron, 
c la vitesse de la lumière. 


TI. Calcul pour un cristal quelconque. — 1. Cas 
GÉNÉRAL. — La formule précédente a été établie 
dans l’hypothèse d’une vallée unique, nous allons 
maintenant considérer le cas où la surface d'énergie 
constante est constituée de plusieurs ellipsoïdes 
dans l’espace des k. Ces ellipsoïdes se déduisent de 
l’un d’entre eux en le transformant par toutes les 
opérations qui conservent la symétrie du cristal. 

Nous calculerons le courant total créé par le 
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champ E en présence de H en admettant que les 
différents vecteurs densité de courant J obtenus 
pour chaque vallée s’ajoutent purement et sim- 
plement. Ceci résulte de l'hypothèse faite initia- 
lement de l’indépendance des vallées et de l’inexis- 
tence des transitions intervallées. 

Connaissant la position et la forme d’une vallée 
par rapport aux axes du cristal, la formule vecto- 
rielle ci-dessus permet de calculer le vecteur Jet de 
le rapporter aux axes du cristal, la sommation sur 
les différentes vallées est alors immédiate. 

Soient u,u,u; les vecteurs orthonormés qui 
servent à repérer le cristal, et +, &, ©, ceux qui 
servent à repérer la vallée définie en forme par 

ne 
l'opérateur «. 

On passe des axes # aux axes par la relation 
matricielle 


D; — Ai ui. 


Comme la matrice À ainsi définie est unitaire, 
on à 


Uri — Bit avec B}, Æ Aÿ, 
> 


> 
Les composantes de la matrice « sont données 
dans les axes & par les relations 


2j = Fe 

Mi 
où M1 M, M; sont les valeurs des masses effectives 
prises dans les directions principales de la vallée 
considérée. On obtient les composantes de la ma- 
trice « dans le système des axes w en transformant 
par la matrice unitaire À ce qui donne : 


k yi 
[A—1 o AM = X Ai Ai, 


HET: 


On peut traduire dans les axes w une relation 
—> 


— 
de la forme J — & E par la formule 


E = > AŸ Ai Ex 


ki mi 


En transformant ainsi les différents termes de la 
> 
formule J — ë (H)E puis en faisant la somme sur 
toutes les vallées, on obtientle courant total produit 
par un champ £ donné. Nous allons appliquer cette 
méthode à quelques modèles simples quant à la 
forme et la disposition des vallées. 


A) Modèle « S ». — Dans ce modèle, on admet 
que les différentes vallées sont représentées par des 
ellipsoides de révolution dont les axes coincident 
avec les axes 100, 010, 001. 


On posera m, — Matt my — Em, de cette façon 


en faisant Ë — 1 on doit retrouver le cas de la 
symétrie sphérique à un facteur 6 près. 

Le calcul fait en annexe montre que le courant 
total peut se mettre sous la forme 
J=06E+o, H/\E 


+ [AE.H) H + LE? + VT(E, H)] 


avec 

AR 
EP (1 + 26) MERS 

mi 
6, —= LE DAC AI TA 1 ed 
Î m2 ( UAE RÈTTS +E+2 
Re ee VE pie 
é mn ë 


On remarque que dans ce modèle il existe une 
relation simple entre les coefficients des termes de 
magnétorésistance, savoir : À +u + v = 0. De 
plus, si les ellipsoïdes deviennent des sphères £ — 1 
le coefficient vs’annule (à +u — 0) et l’on retrouve 
la formule connue, valable dans le cas des cristaux à 
symétrie sphérique pour lesquels les termes de 
magnétorésistance peuvent se mettre sous la forme 


Ju = o{E À H) A H. 
En appliquant les formules générales démontrées 
en annexe 6 qui ne sont d’ailleurs valables qu’au 


second ordre en H, on peut écrire le vecteur E en 
fonction des vecteurs J et H suivant la formule 


E= pol + pJ AH 
+ e[L(J.H) H + MJH? + NT(J, H)] 


où 7 représente le même tenseur que précé- 
demment. Les coefficients sont définis par : 


/ m\ 1 1 


{ 


M 


MT LT 5 > EE 
rie 2 €? il 
PRE 27 (2) 
[7 9 SR Die 
A EU ii 26 
a A 
M =— STE 
1 DUMP AN 0 
< 72 2 (2 + E)2 
| Fr. a n 
ET AL RUE 
RTE 
\w =y=E+,—2 
\ <) 
B) Modèle « G ».— Nous étudierons maintenant 


un modèle dans lequel les surfaces d’énergie cons- 
tante sont constituées par huit vallées ellipsoïdales 
de révolution dont les axes de révolution sont 
orientés suivant les directions (111) et dont les 
centres sont situés aux limites de la zone de 
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Brillouin, Le calcul fait en annexe permet d’écrire 
comme précédemment 
J=0E +6, H\E 

FE. H) H + LEH? + YT(EH)] 


avec 
Eu" nil 
kon . 1Ë +2 
le RE ee 
#0 ET gr DURE 
An x 186 +29 
= 2- et > EU + 26) bee 
3m É 
2 (E— 1} 


RNA > ea 2erOos \ tee . 

: A Le FL ERNeSl J 3 E(1 + 2€) 

Il existe, ici encore, une relation simple entre les 
coefficients soit : 


À + u = 0, 


B. Effets thermomagnétoélectriques. 


I. Principe du ealeul. — Suivant la méthode uti- 
lisée précédemment pour calculer la densité de eou- 
rant, nous calculerons le flux d’énergie transporté 
par les porteurs en mouvement. 

Avec les mêmes notations qu’au paragraphe A, 
le flux d'énergie transporté par un courant corres- 
pondant à une concentration de n électrons par 
unité de volume et pour un processus de conduc- 
tion isotherme s’écrit : 


n | EVE — AË) dk 
W — al : 


fre se Aë) dk 


Ne conservant que le premier ordre en E nous 
aurons : 


n | 6V AGf/(6) d’k 
fre dix 


relation qui, tous calculs faits, pourra s’écrire : 


SE. 

ES 
W — K(H)E 
ES 


me 
où K (H) est un tenseur du second ordre dépen- 
dant du champ magnétique H. 

Si la structure de bandes possède plusieurs vallées, 
la 1% vallée contribuera au flux d'énergie par la 
quantité Wi et le flux total d'énergie sera la somme 
des flux partiels W : 


W—Ï>w, 
à 


la somme étant étendue à toutes les vallées. 
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Le flux d’énergie W; sera de la forme : 
a 
> 
W; — K (HE. 


— 


Donc le tenseur. K (H) sera égal à la somme des 
ie 
tenseurs K,(H) : 
> > 
K(H) = D K(H). 


a YŸ 


Or, nous savons que le’tenseur de Peltier x (H) 


est défini par la relation 
> 
+ 
W = x (H)J — 


Eee 
Connaissant le tenseur K — 
calculer le tenseur de Peltier 


> > 
— 


> 
se 
Tr (H) = K(H) o—1 


(H). 

Puis, utilisant les relations d’Onsager, nous éta- 
blirons que le tenseur « pouvoir thermoélectrique » 
peut se calculer par l'équation tensorielle suivante : 


52 ee ie 
z AH) COR CC HS 
Gun = à en. 


où Test la température pour laquelle, nous avons 
calculé les tenseurs de conductivité et de Peltier. 


IT. Calcul du flux d’énergie pour une vallée. — 
Utilisant le résultat établi dans l’annexe 3, nous 
allons donner une relation qui nous permettra de 
simplifier le calcul de W. 

Nous avions montré que l’on pouvait écrire 
l’équation suivante : 


Le À fe. U] g(8) f(8) dx 


Î 18) dk 


où U est un vecteur indépendant de k. 
Le vecteur flux d’énergie s'écrit alors sim- 
plement : 


U < (6) > 


ie 
= TS se) 
T6 1 (E ) 
Se es > Elan HAE 
T° 6 1 =e | 
De FT oË nie [5 HEH — w? à E]. 


— 


— 
Le tenseur X(H) défini précédemment est done 
le tenseur faisant passer de E à H en utilisant 
l’équation vectorielle ci-dessus, 


N°4 
TT: Calcul du tenseur « pouvoir thermoélec- 


+ 
trique » : Q.-— Nous savons [10] que les «flux »: 

J : densité de courant électrique, 

J, : densité de courant d’entropie 
et les «forces » : 

Vu : gradient du potentiel électrochimique. 

VT : gradient de température 
sont les flux et forces conjugués des relations 
d’Onsager, c’est-à-dire que si nous écrivons un sys- 
tème d’équation de la forme : 


> > 
_ LS 
JET n(H)VH dr Le(H)VT 
= > 
5 BUS 
J = Lo(H)Ve + Lo(H)9T 
> > 
oe ° <a ZA . 
il existe entre les tenseurs L,, et L,, la relation 
d’Onsager 
> & 
e + 
Lo{H) = Lo(— Hi) 


Nous pouvons écrire pour J et W des relations 
semblables à celles que nous venons d’écrire pour J 
etre 


ee 
1(H)X; A M, (H)X; 
= 
+ 
- M 


2o(H)X2. 


Les relations d’Onsager s’appliqueront alors aux 
> > | 


+ + 
tenseurs M ,,(H) et M ,,(17) sous la forme : 


= 5- 
_ > 
Mot) = Ma(— H) 


si X, et X, sont les « forces » conjuguées aux 
«flux » Jet W. Pour les déterminer nous devons 
écrire que la source d’entropie est inchangée : 


ice 


= VE + NT = IX: + WXe 


Le flux d'énergie W étant relié au flux d’entropie 
et au courant électrique J par la relation : 


J 
W = Ts —p° 


l'identité ci-dessus nous permet d'écrire : 


VTT 
X vu Ve 


VT 
i 


X, — 
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Utilisant ces résultats, les relations pour J et W 
se transforment de la façon suivante : 


Le . Hj> © 
J = Ma 4)ve + (Mu +, ] 7 


avec 


> > 
M 12(H) = Mat— Hi). 


Interprétons à la lumière des résultats des para- 
graphes précédents les différents tenseurs inter- 
venant dans ces relations. 

Si la température est maintenue constante 
(VT = 0)les courants J et W sont donnés par : 


+ > 

La Le 
J =M;,, Ve =eM,E 
> 

> ‘ > 
W=—M,Ve= eM,E 


Ÿ 


donc 

> > 

> 

M 1H) = (Je) © (H) 
et 

£> se 

Ee 1+ 

Ma(H) = = K(H). 


Sinous imposons maintenant que le courant élec- 
trique soit nul (J = 0), il apparaîtra un gradient 
du potentiel électrochimique qui sera relié au gra- 
dient de température par l’équation : 

_ > 
Èe É na vT ik 


> 
> 
E 1] LATE 


M::(H)ve + |M,.(H) T 


e 
Le tenseur « pouvoir thermoélectrique » défini 
par 


> 
> 
ve — Q(H)UT pour 


HN) 
peut s'exprimer sous la forme : 
> > 
< Lee Na (H) M UT 
se 11 12UÆ) 
C — = = — . 
à (HT) 27 | T 


En tenant compte des relations précédemment 
établies, nous avons : 


_ ie 
z U ce —UH) À (— H) 
QU) Hi - 
re _ & 
# E Le o—1(H) r (— H) « (— H) 
eT T 


256 JOURNAL DE PHYSIQUE 


Si la structure de bande présente plusieurs 
vallées, le tenseur « pouvoir thermoélectrique » se 
calculera par la relation suivante : 


> se Ke — 
> 6—UH)DKi—H) v > 
Q(H) — T re 


Ja somme Ÿ étant étendue à toutes les vallées. 


IV. Applications. — A. Nous appliquons Ies5xe 
= 


Ÿ 


> 
sultats précédents à un modèle isotrope (4.= 1 /m) 
en l’absence de champ magnétique. 
Dans ce cas, nous avons (dans le système 
CAGRS A 
JE 


me 


W = (ne]m) < 7567 E 


donc : 


Considérons maintenant le cas particulier d’un 
semi-conducteur non dégénéré dans lequel le libre 
parcours moyen des électrons serait proportionnel 
à 6". Le temps de relaxation serait alors propor- 
tionnel à 81/2, On calcule facilement que : 


ÉREE 


a 


et 
& — ÀT Log (n[M) 


où M, est le nombre d’états de la bande de conduc- 
tion. Le pouvoir thermoélectrique s’écrit ainsi : 


he = (2 + r — Log (n[Mc)]. 


Ce résultat bien connu, nous permet de vérifier la 
validité de notre théorie. 

B. Nous allons maintenant appliquer les mêmes 
résultats au cas du modèle «S » pour lequel nous 
avons donné précédemment les tenseurs de résis- 
tivité et de conductivité. 

ONE - > 
> > La 
Il faut calculer le tenseur À —6c —1(H) K (— H) 
— 
à partir des relations de définition W — KE et 
E — 67 (H)J. En combinant ces deux relations, 
on montre qu’au deuxième ordre en # on obtient 
> 
AV= aV+aHAV ; 
+ alp(V.H)H+ q VE? + rT(V,H)] 


N° 4 


où V est un vecteur arbitraire et où les coefficients 
do ds An P, q, r sont donnés par les relations 
suivantes : 


rot 
Dpt 2 Zee 
: NE MN Ar 
ne My 1 HIIE VOLRT> ae 

6 2 LISE > DEMI 
SE 7 À + 26 a TT 


NA LEE 26) 7 
RITES 
f<mE > TES 
cr €? 
ma +orfi+E+e) Se 
Aro FE) 
RTE Ne 
2 2 FERRER Te 
d miUE Eee, 
LT LT > 
| HER 
2 1 EE 
air = (EE d) den 
RE ET TC 
Ne tee à 
On vérifie encore ici la relation simple 


p +g+r —0. De À on déduit facilement le 
tenseur de « pouvoir thermoélectrique » par la 
relation : 


> > 
+ + 
VA 


Les coefficients à, 4, a, pqr sont donc simple- 
ment liés aux effets thermomagnétiques, tels que 
l'effet Nernst. 


V. Conclusion. — La méthode de calcul des phé- 
nomènes de transport que nous avons développée 
permet de trouver les relations entre les différents 
coefficients qui décrivent les effets magnétoélec- 
triques et thermomagnétoélectriques. Les calculs 
faits sont transposables directement au cas de tous 
les modèles comportant des vallées elliptiques. 

Comme les valeurs moyennes de + qui inter- 
viennent dans les différents coefficients ne sont pas 
identiques, on peut espérer que les résultats expéri- 
mentaux complets permettront de choisir entre 
différents modèles de structure des surfaces d’éner- 
gie constante et entre les différents processus de 
diffusion possibles. 
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ANNEXE 1 


>> > 
++ _ 
Caleul de x 6 k et B? x 


Pour effectuer ces calculs, nous choisirons notre 


système d’axes tel que la matrice « soit diagonale. 
>> 


+ F7 es * 
La 2% composante du vecteur & 8 k s’écrit : 


Ê> } A 
a Bk “mile — à D €ù à; À Ht 
jt 


en utilisant l'indicateur de Kronecker : 

en — 0 : si 2 des indices sont égaux ; 

ej — +1: si(z7l) est une permutation paire 
des indices (1 2 3) correspondant au sens direct des 
axes ; 

eh — — 1 :si(il)est une permutation impaire. 

Nous pouvons écrire : 


6> ) 
«x Bk = Ye oi «5 o4 ki ay + Hi, 
il 


Remarquant que 


Ë 2 
HROROTENUEUNS/— PCT — 
== 


= 


. a , . 
notre expression de « f k s’écrit : 
=) a 
(; 23 RS an = — FA el 
D'où finalement la relation vectorielle : 
>> _ 
++ + 
x Bk — 1(/P) F4 on. 
= 


> 
Cette formule nous permet de calculer 8 ? k 


G%k — 6 Bk — En) n H=U/P)\R\«— H} \H 
> > 
ee en He 4 TL (5 al 
ANNEXE 2 
> >] —1 
+ + 
Caleul de l’opérateur |; — 7 | 
se 
L'opérateur B possédant les valeurs propres 0 
la et —io, nous savons que l'opérateur 


> > 1-1 
> > l 
[5 — T7 | peut se décomposer sous la forme : 


= > 
> Dit > + 
(5 +) des 
ARR NÉS PERMET 
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J 

> > > 

+ + $ 
où E, A,B sont des opérateurs que nous allons 
calculer. Pour ce faire, nous allons multiplier à 
droite, les ur membres de l’équation précédente 


Se 

ds CL ram | 
par \1 — + $/etidentifier en +,ce quinous conduit 
aux relations suivantes : 


> > 
£ À 

CA 
E=— 2 


> >> 
+ >+ 
AB — 


oO 


EYY 


> > 
+ > 
B=0  B 


> + > 

de à ++ + 

En utilisant l’équation 8 (Res «°? == 01e 

deux premières relations nous donnent en élimi- 
> 


PA 
nant B: 
> > 
> + 
A= 6 
et par suite 
ES 
> > > > + 
+ + RE 
B=—$Gja et Ë=1 +2 
D'où le développement cherché : 
> 
> œr 
> mn Noire po 
(F #4 RC 2 
A MS bete 
«e Z ee 
À + + É° 
ce CET À + © T7 
ANNEXE 3 


Calcul d’une fonction de la forme 


[Ex [k.V] elE) €) dk 


fre) ax 


où U est un vecteur indépendant de k 


et 


LS 
6 = É k).r les intégrales sont étendues à tout 
l’espace des k. 

Remarquons tout d’abord que l’élément de 
volume d°k peut s’écrire 


dSé dé 
Av 6| 


où dSg est un élément de surface de l’ellipsoïde 
6 — Cte et 


dk — 


= 
<= 
Ve 6 = & ki 
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L'intégrale du numérateur s’écrira alors : 


fie. 0) ete) F8) àse dë 


Le as 8(6) f @a8 f," n[k.U] dS£ 


> 
ak 
en posant : n — = os — normale extérieure à 
k 
l’ellipsoïde, & — Cte, : surface de l’ellipsoïde 
considéré. + 


Le théorème d’Ostrogradsky permet d’écrire la 
seconde partie de l'intégrale précédente sous la 
forme : 


k.U1 ds = RU Heu) dk 
[" Rose AT ER Q(8) 


Q(6) étant le volume de l’ellipsoïde limité par la 
surface Sg. On caclule alors facilement que 


dk — (4x/3) VP 6/2, 
fe 4) c 


Finalement l'intégrale du numérateur de 7 
s’écrit : 


Se 
> 
JE k[k.U!] 2 


8) f'(8) dx 
= EU VP [07 er g(8) (8) à8. 


Quant à l’intégrale du dénominateur, il est facile 


de montrer en utilisant le système d’axe qui diago- 
ES 


—> 
nalise «, qu’elle peut se transformer sous la forme 


2RVE [JE 


D'où l’expression de la fonction de U cherchée : 
7 GAZ 6) dé 


) 8/2 dé. 


1/2 
1 8"? f(6) d8 
ANNEXE 4 


Caleul du courant dans le modèle « S » 


Dans ce modèle les différentes matrices À qui 
définissent les vallées peuvent se représenter par : 


a b c 
ARR c a 
C a b 


où a, b, c sont trois nombres dont l’un est égal à 
plus ou moins un et les deux autres nuls. On com- 
mence par faire le calcul avec la matrice À ainsi 
déterminée puis on fait la sommation sur les diffé- 
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rents cas possibles : savoir : a= +1 b=c—=0, 
a=c=0 b=+la—b—=0 c—= +. 

On voit facilement que lon pourra par la suite 
négliger les termes qui contiennent le produit de 
deux lettres a, b, c, différentes et ne conserver que 
ceux dans lesquels ‘chaque lettre figure seule à une 


‘puissance paire. 


Les seuls termes non nuls parmi les [AT1& AŸ 
sont 

[A—'aA]" = (a?]m;) + (b?]me) 

LA! A]? = (b?fm:) + (fm) + (a? /ms) 

[AT AIS = (otfm) + (at{ms) + (B ma). 


Évaluons les différents termes qui composent J: 
S 


> 
{Terme : J = na E<T>;la formule générale 
donne après sommation 


+ (c?]ms3) 


| 
m; Ma m3 
=> 
ere 
28 terme n << 72 "on transiorme 


M Mo Ma 
«1 H par la formule de changement d’axes précé- 
dente où il suffit de changer 1/m; en m; et on 
obtient : 


M + Mo + M3 HI | E 


M1 Mo Ma 


Jr = 2n < 7? > 


32 terme : Le premier morceau (E.H).H ne dé- 
pend pas de « la sommation est immédiate. 
Le second morceau contient 


= 12) ne 


calculons-le pour une vallée : 


1 
= [let 


2 12 2 
nn Ma + b?ma + cm3) H? 
+ (b%m: + cms + am) H2 + (cm, + am + b? ma) H2]. 


On peut alors évaluer la somme de termes de la 
= 


ES 
forme ne ” on obtient un vecteur de compo- 
santes : 


MOT Te TTL 


| SE, H3+ (+20) r re 


Fa (= Ces . Ha _ E, I? 


M, M 


MA Te ITS ; 
(+ 7e 2 FE, H? + 3E, H3 
HO Va ie 


mm m 2 M3 
(nes) pere 
Mo M3 m1 


mi Mo M3 
(se) 7, 
Hs DES OR 


+ (5 +742) 5, #3 + 38H 


HS ln ep 


N° % 


en posant m, — m3 et m, = Emil vient : 


ll 
SE, H3 +(1+6+5)Am 41e 


œ 
E {ir + 


c’est-à-dire un vecteur que l’on peut écrire sous la 
forme 


Ai 


(l 
+ :) E, H? + 3Æ, Hè + (: DAvEs :) E, 12 


Ÿ 


HR 


41 
+ j) EH + (1 LEE 1) E, H? + 3E, Hi 


(1 eee à E.IE® + (2 ke : —:) T(E, I) 
€ (es 


où T(E, H) représente un vecteur dont les compo- 
santes rapportées aux axes w sont Æ, H?, E, H3, 
E, H? (E; et H, étant les composantes de E et H 
rapportés aux mêmes axes). 

En regroupant tout ce qui provient du troisième 
terme, on obtient : 


Lors > 


mi-m3 


SES GE.H) 4 — (1 +5+3)AE 
« 


+(+5—2) Ten] 
Le 


Le vecteur courant total peut donc se mettre 
sous la forme : 
J—6, E +6, H\ E 
+ [AE.H) H + LEH? + VT(E, H)] 


avec 
1 2 NZ $ 
ner | +) = 2» (1 + 26) 
He mi 
2M ER EE 
D him <>. 
M1 m3 m2 
Î ÉS 
NEA RE — g2 
M1 M mi 
D REES 
uw =—{1+€+2) 
y =E+5—2 
e 
on remarque que 
À +u + v—= 0. 
ANNEXE 5 


Calcul du courant dans le modèle « @ » 


Nous suivrons un plan analogue à celui de 
l'annexe 4 sans reprendre les différents raison- 
nements,. 
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Les matrices À sont maintenant définies par : 


0 nS2 _Æ5 

2 V3 V3 

€ ET EC € STONE 1 

À — 1e + | | LS | — — 
2 V3 2 V3 V3 

CAN ET | ef 1 Es Es | 1 

2 V3 2 V3 k V3 


OÙ € € Eg, SONt égaux à + 1. La sommation sur 
les vallées revient ici à sommer sur les huit combi- 
naisons possibles pour l’ensemble €, €, €3, puis à 
diviser par deux pour tenir compte du fait que la 
première zone de Brillouin ne contient qu’un demi- 
ellipsoïde de chaque vallée, 

En tenant compte de la symétrie de révolution 
des ellipsoïdes (m2 — m3) la matrice transformée 
de « s’écrit : 


CRT) 
à Fe A SC 
3 Mi Me 3 [mi au os 
ORAN 


£o 
| 
E Î 
SU TT nt: Mi - m | ce 1 | 


3 
monde aa) 
D TC 0 3 [mu 3 
On procède comme dans l’annexe 4 pour PRE 


les différents termes. 
1 terme : 


Jo — 
22 terme : 


I=T< 2 > [ete] HA E. 


M Mo 


meer 2] E 
3 na N ne NE 


3€ terme : 
(JP) Ha H) 


A m1 + 2m 
= - ——— H2- 
3  Mims 3 


[es €o Hi Ho + 


I 


o 
2 M1 — Me 
M: M 


etes Hi He cntEs Ho ds] 


a 
(s Où =) " 


EM CE De 1 
9 mimi \m: 


2 
)£: H? 
Mo 


8 = 1 1 
T9 . 5 É m LS Ha Bent HMS A5] 
J 1 Mo 1 2 
( _ ] ï +2 1 # 
2 mm; me ( Fa ] E, H? 
> [O7 4 12) 2 9 M ma m Î Mo 2 


+ 8 Mi — Ma ( 1 


) [H HEs+ H,H,bi 
Mo 


9 mimi \m: 


24 
(s &2 à a 


 _8mMmi—mMaftl 
ÉT 


EAU 2) : 
2 k Ma + 1 a ; ] E, H? 
mi 


O7 m5 Mo 


1 
—) (Ho HE + H3H3 


mm? \m: 
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km, + 2m 1 2 
L — | E H° 
9 mms Mi M 


8 M1 — Ma 


LA [(H(E.H) — T (E, Hi] 


Or MENT Tite 


ms (E.H) 1), 


3 
M M 


Le courant total peut donc encore se mettre sous 
la forme 


J—0,E +6, H\ E 
+ [XE.H) H + LEH? + VT(E, H)] 


où Ta la même signification que dans l’annexe 4 
et où l’on a : (avec Ë — m,/m) 


m=;<r>n[e ++ 


Mi Mo 


SE 1 o) 
ju <= > (1.+ 24 
! 1 DTA 
nee net |. 
3 Ms Mol Mo 


Fee <#>n[ft+i]s 


N° 4 


On constate qu’il n’y a pas de relation simple 
évidente qui relie À y v indépendamment de E. 


ANNEXE 6 


Caleul du champ électrique E 
en fonction de la densité du courant J 
et du champ magnétique H 


Nous avons trouvé une expression de J sous la 
forme 
J—-0E+6, H\ E 
+ SAE.H) H + VH°E + YT(E, H)]. 


Paritération jusqu’au second ordre en H, on voit 
facilement que E s’exprime en fonction de J et 
de H par une relation de la même forme que celle 
donnant J en fonction de E et de H. Nous posons 
donc : 


E=poJ+p HA 
+ P[L(I.H) H + MH? J + NT(J.H)]. 


En portant cette expression de E dans celle 
donnant J et en identifiant on trouve : 


3 TNT 0 = 1/00 
= <> EN + 20 P1 = — Po 51 = — C1l0$ 
‘ 
à ; . Pa = — P O2 — — 02/06 
Î Mo 4 ?) 
PA = —— [3 + = (m; — Ma) ( =_)] 
Mo + 2m: 3 Mi M ee 00 5? He 6? 
GR C2 Go S2 
orne 2 2 
Me Po NT . 61 
+ Lien MR re ce 60 Se 
Li 3 M; 3 él N = ÿ 
(PE 2 (ms — ma)? 2 LE En We Manuscrit reçu le 11 août 1959. 
\ 3 MmilMo + 2m) 8 E(1 + 2Ë) Modifié le 25 septembre 1959. 
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LETTRES A LA RÉDACTION 


TEMPS DE VIE DES NIVEAUX DU !l$n (*) 


Par Rae STIENING et Martin DEurscu, 


Laboratory for Nuclear Science, 
Massachusetts Institute of Technology, Cambridge, 
Massachusetts. 


Meessen et Grenacs [1] ont signalé des coïncidences 
retardées entre les rayons gamma du Sn produit par 
la désintégration du 11In+(53m). Selon leur mesures, 
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IG 


(*) Ce travail est subventionné par des crédits de 
l’Atomic Energy Commission, l'Office of Naval Research 
et le Air Force Office of Scientific Research dans le cadre 
du contrat AT(30-1)-2098, et du contrat AF 49(638)-181. 


une raie avec une énergie d'environ 2,15 MeV suit la 
raie de 0,41 MeV avec une période de 10 + 1 ns. La 
cascade retardée constituerait à peu près deux tiers de 
l'intensité globale du rayonnement à 2 MeV. 

Nous avons cherché à vérifier l’existence de ces coïn- 
cidences différées en utilisant un dispositif perfectionné 
dans notre laboratoire depuis un an. Les détecteurs 
étaient les suivants : 

Pour la raie de 2 MeV, un crystal de Nal(Tl) de 
3 inch de diamètre et d’épaisseur, couplé à un photo- 
multiplicateur RCA 7046. Pour la raie de 0,4 MeV, un 
crystal de ! inch et un P. M. type RCA 7265. 

0,7 mg d’indium métallique fut irradié dans la 
pile MITNR, puis dissous dans l’acide chlorhydrique. 
Les sources furent préparées à partir de cette solution. 

La figure 1 montre le spectre de délais obtenu pour 
la cascade 2,1 MeV — 0,4 MeV et un spectre étalon 
obtenu avec une source de coïncidences quasi-simul- 
tanées. On peut déduire de la coïncidence des deux 
courbes que l’intensité d’une cascade éventuelle avec 
un temps de vie de 10—8s doit être inférieure à 2 % de 
l’intensité totale. Au moins 75 % des coïncidences ont 
une période inférieure à 0,5 ns. 

Dans une deuxième série de mesures, la bande de 
sélection d’énergie fut élargie de façon à admettre 
tous les gammas avec une énergie supérieure à 1 MeV 
en coïncidence avec la raie de 0,4 MeV. Des résultats 
tout à fait analogues furent obtenus. 


Lettre reçue le 8 février 1960. 
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INFLUENCE DES FLUCTUATIONS 
DU NOMBRE DES MOLÉCULES 
SUR LA FRÉQUENCE D'UN MASER 
OSCILLATEUR A FAISCEAU MOLÉCULAIRE 


Par T. C. Wanc, 
Laboratoire de l’Horloge Atomique du G N. R. S$., 


Paris-Besançon. 


Les fluctuations et déplacements de fréquence les 
plus importants ont été traités assez complètement 
dans des articles publiés précédemment [1], [2]. 

La présente note étudie les fluctuations de la fré- 
quence d’un « maser » oscillateur provoquées par les 
variations du nombre de molécules dans le jet. Ces 
variations modifient en effet la saturation de l'émission 
stimulée. Ce type de fluctuations en fréquence n’a pas 
été sérieusement étudié Jusqu’iei. 
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L’équation caractéristique décrivant le fonction- 
nement de la cavité d’un maser en présence d’un fais- 
ceau de molécules exitées est l’équation (39) de la 
référence [2] que nous noterons SWT. 

Cette équation s'écrit : 


Toutes les notations du présent article sont con- 
formes à celles de SWT. En considérant la distribution 
de la polarisation P et du champ électrique E dans 
une cavité résonnante (mode TM), et en suivant la 
méthode de SWT avec une légère modification pro- 
posée par Shimoda [3] on obtient : 


1 PORC) 
5 (1 + 21Q ES ) 
7 RD 
L sin — 
D ne 
Nr (1 cos ©) — : 2e 2 
KAL v « k 
N U 
(1) 
où 4 = /Jw? + 2° et do = © — y. ; 
La partie réelle de l'équation (1) conduit à 
LT? Qn u? AE) 
PAT eu () 
ou encore 
 ., (a HAS La 
ET (uw) nhQ ; ee 
La partie imaginaire donne : 
Sin RL 
&e— © _hr°nu?, L'or 4%0 ) : (3) 
RAA EEE SN k 


L’équation (2) est la même que léquation (29) 
de SWT qui exprime l’égalité de la puissance émise 
par le jet et de la puissance dissipée dans la cavité. 
L’équation (3) détermine la fréquence. On peut 
traiter (2) et (3) comme des relations entre les va- 
riables ©, £ et n. Naturellement on a x — . ou ü 
est également une fonction de Æ. Pour simplifier on 
peut admettre avec SWT que uw est constant. En diffé- 
renciant (24) et (3) par rapport &w, Æ et n et en sim- 
plifiant on peut obtenir les deux relations simultanées 


à dW KL KL dr 
PP SR) a de te a te 
F5 Le 1 + 2(o om) do É Cote, | : 5 0, 


Oc — O9 dr 
do [r .| Æ 


(o — wc) (© — oo) n 
dW kL kL \L 
[21e — &o) do + x? w| [1 = cos - Sin | 
ke? “ ( À 1) 
2 |[— — sin — 
U 


Les équations et (5) donnent les relations entre 
les variations de paramètres (©, W et n) compatibles 
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avec le fonctionnement du maser oscillateur. Élimi- 


nant Le entre les équations (4) et (5) on obtient 


(&o — «c) (© — 0) 
1o— 
Oc — O9 
se cos Le — in he] 
RL SR NS PAREAU : U 
25 CS jee de) 
D — sin — 
V V dn (6) 
kL kL n 
2 Got — À | 
finalement x 


soit ® le crochet de la formule (6). Quand 0 (aui 
est égal à nn) est égal à 0 ou x, D est égal à 1. Quand 


a à : : 
0—= 7 ce qui correspond approximativement aux con- 


ditions normales de fonctionnement d’un maser 
(cf. SWT) D — 1/2. Il est facile de voir à l’aide de 
l'équation (46) de SWT que l’expression 
(© <a Ge) (co sd @6) 
(©ce — ©) 
est toujours du signe contraire à © — &,. La valeur 
absolue maximum de cette expression est obtenue 


| Oe + O9 
OTSQUE © = —5—. 
L’équation (6) donne alors (en prenant ® — 1/2 
q P 
(do) = _ (Ga) 


Ceci est seulement une approximation, pratiquement 
Oe —- [@) 


on ne peut pas avoir © — G ° ; cette condition ne 


pourrait être réalisée que si la surtension de la cavité 
était égale à celle du jet. L’équation (46) de SWT 
montre alors que l’on a : 


DS LL O 


Oe — ©  ‘Q 
L’équation (6) jointe à cette dernière relation donne : 


dn 


do)m: ique — : 
( }max, pratique Q (wo &c) ;: 


(6b) 


Dans les conditions normales de fonctionnement 
, . x . A 
d’un maser oscillateur à ammoniac (6b) le nombre n 
de molécules dans l’état exité est d'environ 1013, le 


dn RE ; ë 
facteur Es est d'environ 3 X 107. La fluctuation do 


est nulle quand la cavité est exactement accordée à la 
fréquence actuelle de l’ammoniac ©,. La bande 
d'accord maximum d’un maser oscillateur conduit à 
une valeur de ©, — «, au moins égale à 8 MHz. Avec 


une cavité à grande surtension mA peut être de l’ordre 
1 


de 1/100. Par suite on a pratiquement : 


(d ©}maxima 
nd, + 13 PR = 12 
= SNS 0 TO (6c) 


pour la raie 3-3 de l’ammoniac, 
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C’est précisément la limite de détection des fluc- 
tuations de fréquence par les méthodes actuelles. Les 
fluctuations changent de phase quand ©, — &, ; ce 
résultat peut donner une nouvelle méthode pour dé- 
tecter l’accord de la cavité. D’autre part on peut 
traiter 5, £ et n comme variables au lieu de ©, et n. 
On obtient un résultat presque identique à celui donné 
par l’équation (6) mais de signes différents. Physi- 
quement il est clair qu’il s’agit d’un glissement de fré- 
quence dû à la saturation. C’est un effet analogue à 
l'effet d’élargissement dû à la pression de Van Vleck 
et Weisskopf [4], mais valable cette fois pour une raie 
d'émission et dû à la saturation. La distribution des 
vitesses des molécules dans le faisceau peut modifier 
quelque peu les valeurs numériques de l’équation (6c) 
mais non la nature de l’effet. 

La plus grande partie de ce travail a été effectuée 
dans les laboratoires de la Compagnie Arthur D. 
Little Inc. L'auteur voudrait remercier le Dr H. O. 
MeMahon pour l’aide et l'intérêt qu'il a bien voulu lui 
témoigner ainsi que M. S. Shapiro pour ses utiles 
suggestions. Il est également reconnaissant au 
Pr C. H. Townes pour les commentaires qu'il a bien 
voulu lui adresser. Il voudrait enfin exprimer sa grati- 
tude au Centre National de la Recherche Scientifique 
pour l'invitation qui lui a été adressée et aux 
Professeurs Kastler et Uebersfeld pour leur accueil 
au Laboratoire de l’Horloge Atomique. 


Lettre reçue le 15 février 1960. 
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MÉTHODE DE DÉTERMINATION 
DANS L’'ULTRAVIOLET LOINTAIN 
DES CONSTANTES OPTIQUES 
DE COUCHES MINCES ÉVAPORÉES 


Par Mlle Denise FABRE, 
MM. Jacques Romaxp et Boris Vopar, 


Laboratoire des Hautes Pressions, C. N. R. $., Bellevue. 


Les méthodes de détermination des constantes 
optiques des lames minces solides évaporées sont nom- 
breuses [1], [2], [3], [4], mais la plupart de ces mé- 
thodes sont difficilement applicables dans l’ultraviolet 
lointain, soit parce qu’elles nécessitent la connaissance 
des constantes optiques du support (qui sont toujours 
mal connues dans cette région), soit parce que Putili- 
sation d’un support transparent est nécessaire, Ce qui 
limite alors le domaine dans lequel les mesures peuvent 
être effectuées. Les méthodes utilisées dans cette région 
spectrale sont peu nombreuses [5], [6] et font appel à 
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des mesures de pouvoir réflecteur sous plusieurs angles 
d'incidence. l’état de polarisation de la lumière, qui 
intervient dans les équations, peut être éliminé [5], 
mais la méthode devient alors peu précise. On peut 
aussi, comme l’a fait Tousey [6], le déterminer à l’aide 
des mesures mêmes de pouvoir réflecteur. Cette der- 
nière méthode, qui semble donner de bons résultats, 
fait appel à des abaques. 

Toutes ces méthodes s'appliquent à des couches 
d'épaisseur uniforme. Par contre, la méthode que nous 
proposons s'applique à des couches d’épaisseur va- 
riable. Cette méthode est une application particulière 
d’une méthode interférentielle [7] s'appliquant en prin- 
cipe aux matériaux peu absorbants : On réalise une 
couche évaporée d'épaisseur variable en fonction de la 
distance à une des extrémités de la lame, selon une 
loi de variation qui est vérifiée ultérieurement par une 
mesure directe des épaisseurs. On effectue alors un 
enregistrement du pouvoir réflecteur en fonction de la 
distance sur la lame, c’est-à-dire de l'épaisseur, et ce, 
pour chaque longueur d’onde du domaine intéressé ; la 
mesure se fait sous incidence normale, pour laquelle 
l’état de polarisation de la lumière n'intervient pas. 
Les enregistrements obtenus présentent une série de 
maxima et de minima d'intensité dont l’amplitude 
s’amortit avec l’épaisseur du dépôt. Pour simplifier 
l'analyse du phénomène observé, 1l suffit de remarquer 
que, les pouvoirs réflecteurs étant généralement mé- 
diocres dans l’ultraviolet lointain, seules participent à 
l’interférence les deux premières ondes réfléchies ; nous 
pourrons donc négliger en première approximation la 
troisième onde réfléchie et cela a fortiori, dans le cas 
d’un milieu d’absorption non négligeable. Soient Z, 
l'intensité de la lumière incidente, Z, celle du rayon 
réfléchi à la surface vide-couche, 7, celle du rayon 
émergeant après réflexion à l'interface couche-support, 
e le pouvoir réflecteur vide-couche, 9” le pouvoir réflec- 
teur couche-support, v et x les constantes optiques 
(telles que l’indice complexe n — v — 7x), et d l’épais- 
seur : 
8rxd 

À Vis 


Ti = pl T2 = pt — p}°e 


La composition de ces ondes donne pour l’intensité 
totale : 


TTL ET, LONTIT cos o 


D étant le déphasage des ondes qui est égal à : 
D— p + — 2vd ; © étant la somme des déphasages 
à la réflexion qui sont absolument inconnus. Fina- 
lement : 


ATX 


— — &TVd 
+ 2 Vop’(1— p)e À cos (e | us (É (1) 


À 
Le rapport À est justement la quantité mesurée 
expérimentalement ; la figure 1 représente l’aspect des 
variations de + en fonction, de d. On voit que À se 


0 . . 7 
présente sous la forme de la somme de trois quantités, 
l’une constante, la seconde exponentiellement dé- 
croissante et la troisième périodiquement variable et 


264 JOURNAL DE PHYSIQUE NES 


d'amplitude exponentiellement décroissante. Les quan- 
tités inconnues, p, p’, x et v sont caractéristiques de 
la couche pour une épaisseur donnée, mais peuvent 
varier, en même temps que la structure de la couche, 
en fonction de d. 

Nous remarquerons tout d’abord que les seuls points 
repérables avec exactitude sur l'enregistrement sont 
les extrêma de pouvoir réflecteur. Or, les extrêma de 
la fonction Æ(d) ne se produisent évidemment pas pour 
des valeurs de d telles que ® — kr; cependant, si on 
appelle À et B (fig. 1) les courbes passant respecti- 
vement par les maxima et minima de À, et « et 
celles passant par les points correspondant à 


cos D — + 1 et —1, on montre que la quantité 

Arr B te 2 
C=—;5—, ne difière en valeur relative de 
6 — ET E que de 


2x? | CR na 
x? + v2 3 : __8xxd |? 
ASE les 


ce qui, dans notre cas, peut être négligé. Nous assimi- 
lons done la courbe € (fig. 1) obtenue point par point 
à partir des extrêma expérimentaux, à la courbe c, 
qui est celle représentant la fonction : 


1 __8rxd 
Re RE PA (2) 


Des courbes À, B et C on tire les grandeurs incon- 
nues de la manière suivante : 

19 Le procédé utilisé pour tirer x, v et b des données, 
implique la connaissance préalable de 9. La déter- 
mination de la valeur de b qui est l’ordonnée de 
l’asymptote commune aux trois courbes À, B, €, 
peut être faite par le calcul (quand il n’est pas com- 
mode de l’obtenir directement à partir de mesures 
faites sur des couches de forte épaisseur). Considérons, 
en effet, l’ensemble des valeurs : GA — R,, GB — R;, 
GC = À’, relatives aux courbes À, B, C (fig. 1) pour 


une même valeur de d; on montre que GI = bp est 
donné (à l’approximation définie plus haut), par (1) : 
26 = _ 2e VAR 


On détermine ainsi les valeurs de b pour une série 
de valeurs de l’épaisseur d, afin, d’une part, d’avoir 
plus de précision, et d’autre part, de mettre en évidence 
une éventuelle variation de p en fonction de d, par 
suite d’une variation de structure de la couche. 

20 Ayant déterminé op, on porte sur un graphique 
les valeurs de Log (R°— +) en fonction de d. On 
obtient ainsi une droite dont la pente donne la valeur 
de x et l’ordonnée à l’origine la valeur de Log p"(1 — p)? 
d’où l’on tire p”. 

30 On obtient l'indice de réfraction v à partir des 
points C, (fig. 1) intersection de la courbe C avec la 
courbe expérimentale, et qui correspondent, avec 


D — (24 + 1) 5: à des ondes en quadrature ; d, étant 


l'épaisseur correspondant à C,, on a alors 


À 
dy+s = dy + Cru 


-On peut porter sur un graphique les valeurs de d, en 


fonction du nombre entier n pris à partir d’une origine 


Lite Ar 


arbitraire ; si la courbe obtenue n’est pas une droite, 
on met ainsi en évidence une variation de v en fonction 
de l’épaisseur de la couche. 

Les essais d'application de cette méthode à des 
couches de fluorure de magnésium ont été très encou- 
rageants et nous nous proposons de donner prochai- 
nement quelques exemples d’application pratique, 
ainsi que la discussion des erreurs. 


(?) Nous remercions M. Damany, qui a suggéré cette 
solution simple. 


Lettre reçue le 29 février 1960. 
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